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The nature of the penetrating component of the cosmic 
radiation has been studied by making G. M. counter 
measurements in a cavern at a depth below rock of ap- 
proximately 60 m water equivalent. Observations on the 
vertical counting rate were made inside and outside the 
cavern with various thicknesses of lead between the 
counters. In agreement with the results of previous in- 
vestigators the counting rate with no absorber between the 
counters is approximately one-twentieth of the counting 
rate outside. The small absorption of the penetrating 
ionizing radiation at this depth in 550 g/cm? of lead shows 
that the penetrating component is to be associated with an 
ionizing particle. It is therefore not necessary to assume 
the transmission of cosmic-ray energy to these large depths 


by nonionizing particles, e.g., neutrinos. Comparison of 
the absorption of the penetrating component inside and 
outside the cavern shows that the penetrating ionizing 
radiation undergoes a hardening on transmission through 
large thicknesses of material. The percentage of the radi- 
ation of the less penetrating type is very nearly the same 
inside and outside the cavern (25 percent to 30 percent, 
respectively). This suggests that such soft radiation is in 
approximate equilibrium with the penetrating component 
even at the earth's surface. The character of the absorption 
curve for the soft component under a depth of 60 m water 
equivalent shows that it must be possible for the penetrat- 
ing component to generate secondaries of considerable 
energy. 


INTRODUCTION 


T is now quite definitely established that the 
general cosmic radiation is composed of a soft 
and a hard component. The experimental data on 
the soft component, particularly the variation of 
cosmic-ray intensity with altitude at different 
geomagnetic latitudes, can be satisfactorily ex- 
plained if one assumes that it is composed of 
electrons and photons which behave according 
to the theory of radiation.+? It is, however, 
quite impossible on this view to account for all 
of the observed intensity at and below sea level 
* A report on this paper was presented at the Washington 
Meeting of the American Physical Society, April 1938. 
1 Heitler, Proc. Roy. Soc. A161, 261 (1937). 


*L. W. Nordheim, Phys. Rev. 51, 1110 (1937); 53, 694 
(1938). 


so that some of the radiation which reaches low 
altitudes must be associated with a more pene- 
trating type of radiation. The transmission of 
cosmic-ray energy through large thicknesses of 
matter has been extensively studied in a number 
of depth-ionization investigations. Such measure- 
ments have recently been extended by Wilson* 
down to depths of approximately 1500 meters of 
water equivalent. In order to study some of the 
properties of the penetrating component we have 
performed an experiment designed to measure 
the absorption of the ionizing particles present 
in the cosmic radiation when filtered by a con- 
siderable thickness of material. 


§ Wilson, Phys. Rev. 53, 337 (1938). See the literature 
there cited for earlier intensity versus depth measurements. 
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Fic. 1. Counting rate at earth’s surface as a function of 
lead absorber thickness. 


THE EXPERIMENT 


The experiment was carried out in Linville 
Caverns in western North Carolina at an altitude 
of approximately 2500 feet. The depth of rock 
above the experimental site was approximately 
seventy-five feet measured by a _ surveyor’s 
transit. The rock is of a limestone formation. 
We have calculated the depth to be sixty meters 
water equivalent, assuming the specific gravity 
of the rock to be 3.28. This calculation must 
necessarily be only a rough approximation be- 
cause of a possible inhomogeneity in the rock 
and because of the sloping contour of the 
mountain overhead. The solid angle defined by 
the rays which pass through the counters was 
such that the minimum distance to the mountain 
side was approximately sixty-four feet ; the maxi- 
mum distance was, of course, considerably more 
than seventy-five feet. There is apparently a 
transition effect iri the first 20 m water equivalent 
in going from air to rock® or from air to water.‘ 
By making measurements at a depth of 60, m 
water equivalent we are fairly certain of avoiding 
complications due to such transition effects. The 
counters and counter circuits were essentially the 
same as used in previous investigations.® Insofar 
as possible the apparatus was mounted in 
paraffin covered ply-wood boxes in which drying 


4 Ehmert, Zeits. f. Physik 106, 751 (1937). 
a 3) Z. Morgan and W. M. Nielsen, Phys. Rev. 52, 564 
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Fic. 2. Counting rate under an additional 60 m water 
equivalent as a function of lead absorber thickness. 


agents were placed to safeguard against diffi- 
culties due to excessive dampness The flexible 
leads to the sealed counter tubes were wound 
with rubber tape. In this work we have operated 
four counters in pairs as shown in the inset of 
Fig. 1. The two groups of counters were placed 
vertically above each other with their axes 
parallel. The distance between the center of the 
upper and the center of the lower pair of counters 
was 63 cm. 

Measurements were made of the counting rate 
with various thicknesses of lead between the 
counters outside and inside the cavern. The 
efficiencies of the counters outside the cavern 
were measured by the method of Street and 
Woodward.® From such measurements we esti- 
mate the over-all efficiency of the apparatus for 
detection of coincidences to be approximately 
80 percent. 


RESULTS AND DISCUSSION 


The curve of Fig. 1 shows the data observed 
outside the cavern. The data obtained inside are 
shown in Fig. 2. The counting rates at the two 
locations with no lead between the counters are 
235 c/hr. and 14 c/hr., respectively. The ratio 
of the corrected outside counting rate (293 c/hr.) 
to that inside’ is in approximate agreement with 


® Street and Woodward, Phys. Rev. 46, 1024 (1934). 

7No attempt was made to measure the counting effi- 
ciency of the apparatus inside the cavern. The error made 
by assuming that it recorded all rays which passed through 
both groups of counters wili not appreciably affect the 
ratio indicated above. 
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the measurements of Wilson* and Ehmert* who 
observe a reduction of about twenty-fold be- 
tween ground level and a depth of 60 m water 
equivalent. 

It is obvious from the curves of Fig. 1 and 
Fig. 2 that the ionizing radiation both outside 
and inside the cavern consists of a soft and a 
hard component. Observations similar to those 
shown in Fig. 1 have been made by a number of 
investigators’ with essentially similar results. 
There are, however, several important conclu- 
sions which may be drawn from a comparison of 
the curves and we shall therefore discuss them in 
some detail. 

The absorption coefficient of the penetrating 
component calculated from the data of Fig. 1 is 
0.0004 cm?/g, a value in reasonable agreement 
with the results of previous investigators.* 

The straight line of Fig. 2 has been drawn to 
correspond to an absorption coefficient of 0.0002 
cm?/g, a value slightly smaller than that given 
by the data of Wilson at approximately the same 
depth. Obviously, the absorption coefficient is 
not here accurately measured. It is to be empha- 
sized, however, that the value obtained applies 
to the ionizing particles present in the cosmic 
radiation after filtering through a considerable 
layer of material. The large penetrating power of 
such particles as here observed, makes it un- 
necessary to assume that the cosmic radiation is 
to any appreciable extent, transmitted down to 
these large depths by a nonionizing and pene- 
trating type of radiation, e.g. neutrinos.® It 
seems clear that our results indicate that prac- 
tically all this energy is carried down by ionizing 
rays and that the role played by neutral particles 
must be a minor one. 

A comparison of the data contained in the two 
curves shows that the ionizing radiation under- 
goes a hardening as it penetrates to lower depths. 
A decrease in the absorption coefficient of the 
cosmic radiation with increasing depth below 
sea level has generally been observed in ionization 
chamber’? and counter* 4 measurements. Such 
observations in the past have been concerned 
with the change in slope of the depth-intensity 


§ Auger, Leprince-Ringuet and Ehrenfest, J. de phys. 7, 
58 (1936) and others. 

* Heisenberg, Zeits. f. Physik 101, 533 (1936). 

Rossi, La Radiation Cosmique (Hermann, Paris 1935). 


curve obtained by measurements of the cosmic- 
ray intensity for various depths of absorbing 
material above the sensitive portion of the 
apparatus. It is clear that such measurements 
give no direct information on the penetrating 
power of the observed ionizing radiation. It is, 
therefore, of some interest to point out that the 
results of the present investigation show that the 
hardening of the radiation is directly to be 
associated with the penetrating ionizing particles. 
This conclusion also follows immediately if one 
accepts the well established depth-ionization 
curves of others and if, as shown above, the 
cosmic radiation is transmitted through the large 
depths of material by penetrating ionizing 
particles. 

We have made an estimate of the amount of 
soft radiation in the two sets of data by sub- 
tracting the contribution due to the penetrating 
component as shown in Fig. 1 and Fig. 2. We 
thus find the percentage of soft radiation to be 
very nearly the same outside and inside the 
cavern (30 and 25 percent, respectively). This 
suggests that even at approximate sea level 
the penetrating component must be very nearly 
in equilibrium with its secondaries. The parallel 
variation in the intensity of the soft and the hard 
components suggest rather strongly that a large 
percentage of the soft component at both levels 
is generated by the more penetrating type of 
radiation. The fact that a soft component is 
observed at all at large depths implies that it 
must be generated by a more penetrating type of 
radiation which as we have shown is of an 
ionizing character. 

An examination of those portions of the curves 
which correspond to the absorption of the softer 
component shows that the half-value thickness 
corresponds to a thickness of 20 g/cm*. Any 
difference in this respect between the soft com- 
ponent inside and outside the cavern is not 
large enough to be disclosed by the experimental 
data. There can hardly be any doubt but that 
the soft component is closely associated with 
showers. We are therefore forced to conclude 
that showers are generated either directly or 
indirectly by the ionizing penetrating component. 
It would also appear that the complexity at any 
rate of the smaller showers under heavy layers of 
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matter is approximately the same as that of the 
showers observed at sea level. The general 
features of this conclusion are in agreement with 
the results of an earlier investigation® on shower 
production under thick layers of various ma- 
terials. We observed transition effects which 
indicated that the showers were electronic in 
character. We suggested that such showers might 
have their origin in the generation of shower 
producing radiation (possibly secondary elec- 
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trons and photons) by the penetrating com. 
ponent. It also follows from these transition 
curves that some of such secondary particles 
must have considerable energy. 

We wish to thank Mr. J. Q. Gilkey of Marion, 
North Carolina for his kindness in making Lin- 
ville Caverns available for these measurements 
and to acknowledge the interest of Professor L, 
W. Nordheim and his helpful suggestions jn 
our work. 
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The probability of the astrophysically important reaction H+H = D+e? is caiculated. For 
the probability of positron emission, Fermi’s theory is used. The penetration of the protons 
through their mutual potential barrier, and the transition probability to the deuteron state, 
can be calculated exactly, using the known interaction between two protons. The energy 
evolution due to the reaction is about 2 ergs per gram per second under the conditions prevailing 
at the center of the sun (density 80, hydrogen content 35 percent by weight, temperature 2-107 
degrees). This is almost but not quite sufficient to explain the observed average energy evolu- 
tion of the sun (2 ergs/g sec.) because only a small part of the sun has high temperature and 
density. The reaction rate depends on the temperature approximately as 7° for temperatures 


around 2-107 degrees. 


$1. INTRODUCTION 


T seems now generally accepted that the 
energy production in most stars is due to 
nuclear reactions involving light elements. Of all 
the elements, hydrogen is favored by its large 
abundance, by its large internal energy which 
makes a considerable energy evolution possible, 
and by its small charge and mass which enable it 
to penetrate easily through nuclear potential 
barriers. Again, of all reactions involving hydro- 
gen, the most primitive is the combination of two 
protons to form a deuteron, with positron 

emission : 
(1) 


H+H=D-+e. 


In fact, this reaction must stand in the beginning 
of any building up of chemical elements; it has 
already been discussed in this connection by v. 


Weizsicker.' However, there seems to be a 
general belief that reaction (1) is too rare to 
account for any appreciable fraction of the 
energy production in stars and that it can serve 
only to start the evolution of elements in a star 
which will then be carried on by other, more 
probable, processes. It is the purpose of this 
paper to show that this belief is unfounded but 
that reaction (1) gives an energy evolution of the 
correct order of magnitude for the sun. 

On the other hand, we do not want to imply 
that reaction (1) is the only important source of 
energy. An analysis of all possible nuclear 
reactions with light elements? shows that the 
capture of protons by carbon and nitrogen will 
also play an important role. It is more important 


ly. Weizsacker, Physik. Zeits. 38, 176 (1937). 
2 Bethe, to appear shortly in the Physical Review. 
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than (1) for heavy stars, less important for light 
ones and about equally important for the sun. 

In calculating the energy evolution from reac- 
tion (1), it must be considered that (1) is followed 
by a number of other reactions which are all 
“fast” in comparison with (1) because they 
involve the emission of radiation or of heavy 
particles rather than of 8-rays.’ The deuterons 
formed in reaction (1) will first capture another 
proton, with y-ray emission : 


D+H =He'*+y. (2) 


The ratio of the rates of reactions (2) and (1) is 
about 10'* which means that, in equilibrium, 
there will be about one deuteron for 10'* protons. 
This makes collisions between two deuterons 
very improbable and thus excludes any ap- 
preciable production of neutrons by such 
collisions. 

The fate of the He* formed in (2) depends to 
some extent on the question whether this nucleus 
or H® is more stable. The most recent determi- 
nation’ of the energy evolution in the reaction 
H?+H?=He'’+ 2! seems to show that He’ is more 
stable while older determinations> gave the 
opposite result. The processes which occur are : 

If He* is more stable: 


He*+Het= Be’+y, 
(3) 
Li?+H = 2He’*. 


If H* is more stable: 


He®+e- 


H?-+H =He!+y. (4) 


The net effect is in both cases the same, viz. the 
combination of four protons and two electrons to 
form an a-particle. (The emission of a positron is 
equivalent to the consumption of an electron 
since the positron will ultimately annihilate an 
electron.) One a-particle is formed for each 
process (1). The energy produced per process 
(1) is therefore equal to the difference in weight 
between four hydrogen atoms and one helium 
atom, viz. (cf. reference 5, 6) Q=4-1.008 13 


‘Cf. Bethe, Rev. Mod. Phys. 9, 188, § 76C (1937). 
‘Bonner, Phys. Rev. 53, 711 (1938). 

(gan, Livingston and Bethe, Rev. Mod. Phys. 9, 373 
* Bainbridge, Phys. Rev. 53, 922A (1938). 


— 4.003 86=0.028 66 mass units, i.e., 
Q=4.3-10~ erg. (5) 


While the energy evolution in (3) and (4) is 
exactly the same, the processes lead to the 
formation of different intermediate products 
which could be of influence on the building up of 
heavier elements. However, it can be shown? that 
the building up of elements heavier than Het‘ is 
negligible under any circumstances so that the 
question of stability of He* or H® is actually not 
important. 

In order to calculate the probability of the 
proton combination, the following factors have to 
be considered : 

(1) The probability of collision of two protons 
which involves the penetrability of their mutual 
potential barrier. This can be calculated very 
accurately since the force between two protons is 
very well known from scattering experiments.’ 

(2) The probability of emission of a positron 
during the collision. This involves a theory of the 
B-decay. All available experimental evidence 
points to the conclusion that the original Fermi 
theory gives good agreement with the experi- 
mental dependence of lifetime on energy. 

(3) The energy distribution of the protons in 
the star which is given by the Boltzmann law. 


§2. THE PROBABILITY OF POSITRON EMISSION 


According to Fermi’s theory, the probability of 
B-emission (per second) is 


B=gf(W)|G|*. (6) 


Here G is the matrix element of the nuclear 
transition, 


G= S (7) 


y; and yy, the initial and final state of the 
nucleus, except for the substitution of one 
neutron by a proton. f is a function of the energy 
W of the emitted 8-particle, 


3 2 


1 
f0W) = 
30 20 15 


1 
log (W+(W?—1)!], (8) 


7 Tuve, Heydenburg and Hafstad, Phys. Rev. 50, 806 
(1936). 
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where W is expressed in units of mc? and includes 
the rest mass. A lower limit for the Fermi 
constant g can be deduced by assuming that G=1 
for transformations such as e+. This 
reaction has a half-life of 11.0 minutes and an 
upper limit of the 6-spectrum (observed) of about 
1.25 Mev, therefore B=log 2/660=1.05-10-% 
sec.!, W=3.5, f(W) =11.86 and 


g=0.9-10- sec.-". (9) 


The (maximum) energy of the positron emitted 
in the combination of two protons follows from 
the very accurately known masses of proton and 
deuteron,® v7z.: 


Two protons = 2- (1.008 13—0.000 55) =2.015 16 
Deuteron = _ 2.014 73—0.000 55 =2.014 18 


Positron energy, including mc? 0.000 98 
(mass units), equivalent to 1830-0.000 98 mc? 
= 1.80 me’. 


For W=1.8, we have f(W) =0.132 and therefore, 
with (6) and (9), 


B=1.2-10|G]? 


A similar value is obtained by extrapolating 
the empirical Sargent curves to lower energy. 
B would be larger, if the matrix element G for the 
transformation N“—>C® is less than one. On the 
other hand, the Konopinski-Uhlenbeck theory 
would give about 10 times less for the decay 
constant. However, that theory seems to give too 
strong a dependence of decay constant on energy 
throughout, and will therefore not be used. 

In calculating the matrix element G for the 
proton combination (1), it must be considered 
that the wave function of the deuteron is sym- 
metric with respect to interchange of space and 
spin coordinates of proton and neutron while the 
wave function of two protons is antisymmetric.* 
In particular, when the two protons are in an S 
state (which is most favorable for their coming 
close together), their spins will be antiparallel 
(singlet) whereas the ground state of the deuteron 
is a triplet S state. The transition is therefore 
allowed only if the Gamow-Teller form® of the 
B-theory is used. This formulation permits a 


(10) 


8 We are indebted to Professor Oppenheimer for drawing 


our attention to this point. 
* Gamow and Teller, Phys. Rev. 49, 895 (1936). 
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change of direction of the spin of the proton 
which transforms into a neutron, and therefore a 
change of the total spin by one unit. Strong 
evidence for this theory is found’ in the short 
life of He®; the 6-transformation He®—>Li® leads 
very probably from a 'S to a 4S state and is 
closely analogous to our process 2H—>D. We shall 
therefore accept the Gamow-Teller theory in this 
paper. Moreover, we shall assume that the matrix 
element G in (10) can be calculated simply by 
integrating over the spatial coordinates of the 
two particles, i.e., that the summation over spins 
which is also implied in G gives unity. This 
assumption does not seem serious in view of the 
uncertainty in the numerical factor in (10). 

If the original Fermi theory were taken instead 
of the Gamow-Teller theory, the transition 
H+H=D-+e* could occur only by virtue of the 
‘“‘small’’ components of the Fermi interaction, viz. 
@heavy particle * Welectron-neutrino- Since the Dirac 
operator @heavy (velocity) is “odd,” its application 
on the *S wave function of the deuteron ground 
state will give a function of *P character. This can 
combine with the *P part of the wave function of 
the incident protons. Due to the smallness of the 
latter and of the e-operator itself, the transition 
probability will be greatly reduced, viz., ap- 
proximately by a factor 


(W/D)*(E/ Mc’) = 10, 


where W=1.8 mc? is the energy evolution in the 
process, [=4.4 mc® the binding energy of the 
deuteron, E~50 kev the kinetic energy of the 
protons and Mc?=931 Mev. In this case, then, 
the energy evolution in the sun due to proton 
combination would be negligibly small. 


§3. THE COLLISION Cross SECTION 


According to the general principles of quantum 
theory, the cross section o for the combination of 
two protons of relative velocity v is given by (6), 
(7) if we insert for y; the wave function of two 
protons normalized to unit incident current. 
Thus we have 


o=gf(W)u| Sv |’, (11) 


where ya is the wave function of the ground state 
of the deuteron and y, is normalized per unit 


10 Goldhaber, Phys. Rev. (to be published). 
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density at infinity. The integration goes over the 
space of the relative coordinates of the two 
protons (or proton and neutron, in the deuteron). 

If we assume a square potential well of radius 
r,and depth Vo, the normalized wave function of 
the deuteron is" 


_ (B/r) exp (—(x—x0)) for r>ro 


(B/r) sin ux/sinuxo for (12) 
with the abbreviations 

x=r/b, (12a) 

b=h( Me)-!=4.37-10-" cm, (12b) 

u=(Vo-e/e)}, (12c) 

(12d) 


M is the proton mass, ¢ the binding energy of the 
deuteron, i.e., 2.17 Mev. Depth and width of the 
well are related by the condition 


cot —1. (12e) 


The wave function of two protons, normalized 
to unit density, has the form e“ at large dis- 
tances. (For considerations of symmetry, see $4.) 
It can be expanded in spherical harmonics and, 
at the small velocities prevailing in stars (~10 
kilovolts), only the zero order term will be 
important. It has, at large distances, the form 


¥p=sin (p— ¢)/p, (13) 


where p=kr, k=Mv/2h (13a) 


and ¢ is a phase (depending logarithmically on 
p). At small distances, i.e., inside the nuclear 
potential well, we have 

w iA 


sin 
ror 


(MD)! 
(14) 


where D is the depth of the potential well 
between two protons, assumed “‘square”’ and of 
radius ro, A is a-normalization factor. Outside of 
the well, y, is a solution of the Schroedinger 
equation in the Coulomb field which may be 
written” 


y,=e'* cos K(F+G tan K)/kr. (15) 


Here F is the regular and G the irregular solution 


" Bethe and Bacher, Rev. Mod. Phys. 8, 110 (1936). 
2 Yost, Wheeler and Breit, Phys. Rev. 49, 174 (1936); 
J. Terr. Mag. 40, 443 (1935). 


of the Coulomb equation which behave asymp- 
totically as sin (p—go) and cos (p—go), re- 
spectively. As can be seen easily, (15) goes over 
asymptotically into (13). K is the phase shift due 
to the nuclear field, fixed by the condition of 
continuity of (14), (15) at r=ro; we have" 


tan K = F*6/[1— FG6] (16) 


dlog F dlogw 
with ) ‘ (16a) 
dp dp ro 


For the small proton energies concerned, F(ro) 
is very small (see below), both because of the 
Coulomb barrier and the small value of &. 
Therefore K will be very small (K =0.0017 for 
v=e/h, ie. for a relative kinetic energy of 
12.5 kev which gives the largest contribution to 
the cross section at a temperature of 5-107 
degrees) and e'* cos K (cf. (15)) can be replaced 
by unity. Then the factor A in the internal 
wave function (14) is 


F(ro) tan K 
ksin 


According to Yost, Wheeler and Breit,” the 
Coulomb wave functions may be written (for 
orbital momentum L=0) 


F=Cp?, (18) 


(17) 


where 
n=e/hv (19) 


contains the effect of the Coulomb barrier, and 

’ and © are slowly varying functions of r which 

have the value 1 for r=0. Inserting into (16), (15) 

we get 

tan K (20) 
A= (20a) 


(= 
d log p /ro 
dlog@ dlogw 
= (14 ). (20b) 
dr dr ro 


 (r>r0), (21) 


= C_[(r0) (1) ] 
r 


sin vr/b 
(r<ro). (22) 


sin vro/b 
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In the last expression, we have introduced the 
abbreviation 


v= (MD)\h-'b =(D/e)!, (22a) 


which will be useful for the integration. \ repre- 
sents the effect of the resonance due to the 
nuclear potential well. 

For slow protons, the wave functions at small 
r are almost independent of the proton energy 
and have the form® 


| (23) 
2!3! 3!4! 
O(y) = —2y'Ki(2y!) =1+y(log vee —1)¢ 
ae (23a) 
sis+1! t+1 
with 
y=2r/a, a=2h?/Me?=5.75-10-% cm, (23b) 
y=0.577--- (Euler's constant). 


The deviations from (23), (23a) for 7=1 are of 
the order of a percent. 

Collecting our formulae, we may now write for 
the cross section (11) 


W) (4xb*)? 
v 
where 
= Ai +Ao+Asz, (25) 
0 0 0 
A\= sin ux sin vxdx, (25a) 
Sin SiN VX» 
x0 
A3=Axo exp (25c) 


x0 
The integrations are elementary and give 


SIN SIN 


1 = 


sin (u—v)Xo 


(u—v) x0 


sin (u+v)xo 
| (26a) 
(ut+v)xo 
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Ae=1+x0+ (2b/a)(1+x0+ 3x0") 

A3=Axo{1+(2b/a)[log (2b/a)+2y—1]As 
—[(2b/a)?3(1+x0+ 3x0") 
+3(2b/a)*(7/3)(1+x0+ 3x0? + 


(26b) 


with 


f exp (—(x—xo))x" log xdx, (26d) 


x0 
fi=1+(1+2x0) log xo — exp (x0), (26e) 
fo= 3+ (1+x0+ 3x07) log xo 


—Ej(—xo) exp (x0), (26f) 
fa= (11/6) +(5/6)x0+ 6x0" 
+(1+x0+ 2x0? + log x0 
— exp (xo). (26g) 


Table I gives the numerical results for two 
values of the radius ro of the potential well, 
viz. e?/mc? and e?/2 mc. The depth D of the 
proton-proton well was taken from Breit, Condon 
and Present.” All other quantities were calcu- 
lated from the formulae in the text. It is seen 
that the contribution of the inside of the poten- 
tial well, Aj, is rather small even for the larger 7, 
which shows that the result will not depend 
sensitively on the shape of the potential well. 


TABLE I. Numerical results for two values of the radius. 


ro =e2/mc? ro =e2/2 mc? 

Xo 0.645 0.322 
Vo (Mev) 20.9 66.5 
D (Mev) 10.3 47.0 
2.94 5.45 
v 2.18 4.65 
(rd log w/dr) 0.236 0.110 
1.050 1.025 
O(ro) 0.769 0.854 
(rd log &/dr) 0.050 0.025 
0.814 0.915 
nN 2.63 4.80 
Ai 0.689 0.277 
Ag 1.949 1.547 
As 1.205 1.030 
(1+x9)(1+ 1.835 1.367 
A 2.84 2.44 
A? 8.08 5.93 


13 Breit, Condon and Present, Phys. Rev. 50, 825 (1936). 


(26b) 
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Furthermore, the contribution A; which is due to 
resonance is smaller than A, which would be 
present even if there were only the Coulomb 
field between the two protons. The final result 
for A? increases somewhat with increasing radius 
of the well. For »=1, a numerical calculation 
gave A°=8.8 instead of 8.1. 


$4. SYMMETRY, STATISTICS, ETC. 


If N is the number of protons per cm’, the 
probability of finding one proton with velocity 
vy; and one with ve, is 


N*¢(v1) ¢(V2)dvidvo, (27) 


where ¢ is the Boltzmann distribution function. 
In calculating the total probability of our 
process, we must integrate over each pair of 
volume elements dv,dv2 only once; if we want to 
integrate over v; and v2 independently, (27) 
must be divided by 2 (this corresponds to the 
fact that the total number of proton pairs is } N*). 
We can transform to relative velocity v=v;—Vse 
and center-of-gravity velocity V=}(vi+ve) and 
integrate over the latter, then we obtain (in- 
cluding the factor } mentioned) 


2N*¢(v)dv (27a) 


g(v)dv=(M/4rkT)! 
Xexp (— Mv?/4kT)4rv°dv. (28) 


Of the proton pairs with given velocities v,Ve, 
one in four will have opposite spins, three in 
four parallel spins. The latter have an anti- 
symmetrical spatial wave function and therefore 
do not contribute to the cross section of our 
nuclear process. The former have a symmetrical 
wave function which is, normalized to unit 
density : 


[(exp 
+ (exp (ike +7k, Te) ) (29) 


Separating off the motion of the center of 
gravity and expanding into spherical harmonics, 
we obtain 


(exp (tkz)-+exp (—ikz))/v2 
=Vv2 sin kr/kr+---, (29a) 


which differs from (13) by a factor of v2. This 
gives a factor 2 in the cross section and, with the 


TABLE II. Energy evolution due to proton combination, for 
pcu? = 10 g/cm’, 


é(million degrees) | 5 10 |20 |30 |40 100 
19.8 15.70) 13.72) 12.56) 10.89) 9.90) 9.18 7.28 
e(ergs/g sec.) 0.0040 | 0.15) 0.76; 2.2 | 9.1 [36 150 


factor } for the a priori probability of opposite 
spin, there remains a factor }. 

The total number of processes per cm* per sec. 
becomes thus 


o(v)dvve. (30) 
By the insertion of (24), (28), (19), this becomes 


M\} 
p’ = ) [ve 
kT/ 


( Mv" (31) 
ex 
hv 4 hv 


The integrand has a strong maximum very 
close to 
v=(4rekT/hM)!. (31a) 


Approximating it, in the usual way, by a 
Gaussian around the maximum, we find 


p’ =167-3-*?- gf( W) (32) 
with 
(33) 


If ¢ is the temperature in millions of degrees, 
we have 
7=33.8¢-4, (33a) 


§5. RESULT 


We insert into (32) the values of the constants 
b (cf. (12b)), gf(W) from (10), A® from Table I 
(for ro=e?/mc*), and assume that the concentra- 
tion of hydrogen (by weight) is cy and the 
density (in g/cm*) is p. Then the number of 
processes per gram per second is 


p= (6.02 - 10%)? 
(34) 


Each combination of two protons leads (§1) 
ultimately to the formation of an a-particle out 
of four protons and two electrons, with an 
energy liberation of 4.3-10-° erg. Therefore the 
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energy production is 
ergs/gram sec. (35) 


The density at the center of the sun" is about 
80, the hydrogen content cy ~0.35, therefore 
pcx’~10. Table II gives the energy evolution 
due to the proton-proton reaction at various 
temperatures. With a hydrogen content of 35 
percent, the central temperature of the sun is 
about 20-10° degrees. At this temperature, the 
energy production is about 2.2 ergs/g sec. This is 
of the same order as the observed average 
energy production of the sun (2.0 ergs/g sec.). 
Thus we come to the conclusion that the proton- 
proton combination gives an energy evolution of the 
right order of magnitude for the sun. 

For a quantitative comparison, it must be 
remembered that both the temperature and the 
density of the sun decrease fairly rapidly from 
the center outwards, and that the rate of re- 
action decreases with both of these quantities. 
The average energy production per gram will 

4 Stromgren, Ergebn. d. exakt. Naturwiss. 16, (1937). 
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thus be considerably smaller than that at the 
center, perhaps by a factor of 10. Since oyr 
calculations are rather accurate, it seems that 
there must be another process contributing some- 
what more to the energy evolution in the sun. 
This is probably the capture of protons by 
carbon (reference 2). 

For many problems, it is necessary to know the 
temperature dependence of the energy produc- 
tion e. It is convenient to express this dependence 
as a power law, e~7™". Then (cf. (35)) 


n=d log ¢/d log T=}(r—2). (36) 


At 20 million degrees, this gives in our case 
n=3.5. This is large enough to make the point 
source model of stars a rather good approxima- 
tion. On the other hand, it is too slow a de- 
pendence on temperature to explain, with tem- 
peratures of the order 2—4-10’ degrees, the very 
high rates of energy production found in very 
heavy stars. However, we believe that our 
process is the principal source of energy in stars 
lighter than the sun. 
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The barytron theory of nuclear forces shows a close correspondence to the theory of the 
electromagnetic field. Estimates based on this analogy are given for the probabilities of the 
processes leading to actual production of barytrons. A comparison with the experimental 
evidence regarding the occurrence of barytrons in the cosmic radiation shows that the theo- 
retical cross sections are too small to explain the properties of the hard component. The diffi- 
culty is increased by the short lifetime of the barytrons which is estimated to be of the order 
~10- sec. from radiative 8-decay. Therefore no simple picture of barytron production in terms 
similar to radiation theory can be given. This failure, however, indicates only the inapplicability 
of perturbation calculations, but does not constitute an actual disproof of the link between 


nuclear forces and cosmic radiation. 


INTRODUCTION 


T present there are two different lines of 
evidence for the existence of a new par- 
ticle, the barytron,' of electronic charge (positive 


* A preliminary report on this paper has been given by 
L. W. Nordheim and E. Teller at the Washington Meeting 
of the American Physical Society, April 1938. 

1 This name was proposed at the Washington Meeting of 
the American Physical Society, April 1938, 


and negative) and of a mass yp about 200 times 
the electronic mass m. 

From cosmic-ray investigations there can be 
little doubt that a large percentage of the hard 
component, which itself constitutes the greater 
part of all radiation from sea level downwards, 
can be neither electrons nor protons? but must 


2 Anderson and Neddermeyer, Phys. Rev. 51, 884 (1937); 
Street and Stevenson, Phys. Rev. 51, 1005 (1937). 
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consist of rays of intermediate mass. An actual 
determination of the mass* from ionization- 
range-curvature measurements has only been 
possible for a few rays but they all show, within 
the rather wide limits of error, masses about 
200 m. 

The other line is purely theoretical and arises 
from the attempts to build up a field theory of 
nuclear forces similar to the connection of the 
Coulomb law and the electromagnetic field. As 
is well known such a theory was proposed by 
Heisenberg and others‘ on the basis of Fermi's 
theory of B-decay. According to this theory a 
heavy particle (proton or neutron) can emit or 
absorb an electron-neutrino pair much as a 
charged particle can emit or absorb light quanta. 
Such a scheme leads to an interaction between 
heavy particles due to a “playing ball” effect, 
which implies the emission of the ‘“‘field par- 
ticles” (electrons and neutrinos) by one heavy 
particle and subsequent reabsorption by another. 
With the f-field as determined from radioactive 
decay, however, grave difficulties arose.° The 
interaction between heavy particles is far too 
small for distances of the known range of nuclear 
forces (~10-'* cm) but diverges at small separa- 
tions at least as r~® because of the simultaneous 
action of two transmitting particles. This be- 
havior can be avoided only by rather artificial 
modifications which cannot be formulated in a 
consistent and relativistically invariant way.® 

It was noted by Yukawa’ as early as 1935, 
before the experimental evidence on a new 
particle, that a part of these difficulties could 
be avoided by assuming the interaction field to 
be associated with single particles of finite rest 
mass (i.e. barytrons). 

It has been suggested by Yukawa‘ and several 
other authors® that the new cosmic-ray particles 


‘Compare the discussion by Corson and Brode, Phys. 
773 (1938) ; Williams and Pickup, Nature,141, 836 

* Compare Nordsieck, Phys. Rev. 46, 234 (1934); Tamm 
and Iwanenko, Nature 133, 981 (1934). 

5 Compare f.i. the report by Bethe and Bacher, Rev. 
Mod. Phys. 8, 82 (1936). 

®For such an attempt compare Camp, Phys. Rev. 51, 
1046 (1937) and Nordheim, Nordheim, Oppenheimer and 
Serber, Phys. Rev. 51, 1037 (1937). 

7 Yukawa, and others, Proc. Phys. Math. Soc. Japan 17, 
58 (1935); 19, 1034 (1937); 20, 1 (1938). 

§ Yukawa, Proc. Phys. Math. Soc. Japan 19, 712 (1937). 

* Oppenheimer and Serber, Phys. Rev. 51, 1113 (1937); 
Bhabha, Nature 141, 117 (1938). 


are identical with those introduced in the field 
theory of nuclear forces. In order to see, whether 
such an interpretation is possible on the basis of 
the ‘‘nuclear force field’ theory now under dis- 
cussion, it is necessary to investigate conditions 


-under which the barytrons can actually be 


created. 

The only method available for this task at the 
present time is a correspondence treatment on 
lines similar to those used in ordinary radiation 
theory. As a logical preliminary step for the dis- 
cussion of the connection of cosmic-ray phe- 
nomena and nuclear forces we have tried 
therefore: 

Firstly, to work out as far as possible the con- 
sequences of such a correspondence treatment. 

Secondly, to deduce from observational data 
the cross sections for barytron production neces- 
sary to account for their actual occurrence in the 
atmosphere. 

The result so far has been a negative one, i.e., 
the estimated cross sections come out too small 
by a factor larger than 10 even in the most 
favorable case (barytron production by photons). 

We show this in Part I by rather general and 
qualitative arguments which should apply to any 
field theory which can be treated in the same way 
as ordinary radiation theory and in Part II for a 
special form of the nuclear field (the one initially 
proposed by Yukawa) for which the calculations 
can be carried out explicitly. 

This failure of the correspondence treatment 
does not imply necessarily that such a con- 
nection between nuclear forces and cosmic-ray 
phenomena is ruled out altogether. It may mean 
only that it is not legitimate to extrapolate the 
nuclear field theory in its present state as 
developed for the low energy phenomena of 
nuclear forces and §-decay to the higher energy 
range of cosmic radiation. In the formalism of 
the theory such a failure will manifest itself in 
the occurrence of divergences in the application 
of perturbation theory. It is known that such 
divergences already occur in ordinary quantum 
electrodynamics, though the omission of the 
divergent terms leads in that case to reasonable 
results in agreement with experiment. From our 
investigation it seems then necessary to conclude 
that this expedient does not work for the nuclear 
field theory and that it is not possible, therefore, 
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to give a simple correspondence picture for the 
production of barytrons. It is not impossible, 
however, that the barytron production can be 
described on lines indicated by Heisenberg’? in 
his discussion of field theories containing a 
fundamental length. 


Part I. PossiIBLE PROCESSES AND ORDERS 
OF MAGNITUDE 


Section 1. General properties of barytrons 


An interaction between heavy particles due 
to the exchange of one single particle with finite 
rest mass gives rise to a potential, the space 
dependence of which is given by 


I(r) =(G?/r) exp (—1r/%); Xo=he/u, (1) 


where G is a constant of the dimension of a 
charge and X the “Compton wave-length” of 
the barytron."' This interaction corresponds to 
a screened Coulomb field. The proportionality 
to r— is the same as in the Coulomb case because 
only one particle has to be transmitted; the 
screening is due to the finite rest mass in contrast 
to light quanta. To understand this factor one 
might say that it represents essentially the 
probability distribution of finding a barytron at 
distance r from a proton or neutron when there 
is not energy enough actually to create it. 
Equation (1) gives the order of magnitude” 
and range for nuclear forces for h~ M/10 and 


A=G2/hce~p/M~1/10, (2) 


where A can be considered as the analog to the 
fine structure constant a=1/137=e?/fc in 
ordinary radiation theory. The fact that A<1 
gives some hope that a correspondence treatment 
of the nuclear field has some meaning. It signifies 
that higher order processes (i.e. those involving 
a greater number of barytrons) should be of 
lesser importance than those of lowest order. 
We shall therefore study only these. In case the 

10 Heisenberg, Zeits. f. Physik 101, 533 (1936). 

11 We denote the rest energies (i.e., mass Xc) of electrons, 
barytrons and protons by m, u, M, respectively. 


12 For a short range potential, zero binding energy for the 
two particles problem (deuteron) results when approx- 


imately 
I(r)dr2he/M; (1a) 


(comp. Bethe and Bacher, reference 5, p. 109), Inserting 
(1) into (1a) the relation (2) is obtained. The actual value 


for A might be slightly larger than w/M. 
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value (2) should turn out not to be small enough 
to ensure this convergence it would be, of course, 
not possible to derive any conclusions by the 
standard methods of perturbation theory. 

In order to have a simultaneous explanation of 
radiative B-decay it is necessary to assume that 
the barytrons have electronic charge (positive 
and negative), obey Einstein-Bose statistics, and 
that they can disintegrate spontaneously into an 
electron and a neutrino. The lifetime can be 
estimated in the following way: The probability 
of finding a barytron in the neighborhood of a 
neutron-proton pair in a nucleus is of the order 
1/10 (ratio of binding energy per particle and rest 
energy of the barytrons). The lifetime ry of a 
nucleus subject to B-decay should be therefore 
~1/10 of the lifetime of a free barytron rz, if 
just the energy is available to create one. Now in 
Fermi’s theory ry is proportional at least to 
where is the available energy (maxi- 
mum energy of electrons). Therefore rg/7y(«) 
~ (u/M)(€/u)*. The lifetime of light elements of 
€)~ 5m is of the order 100 sec. and therefore the 
lifetime of a barytron should be ~ 10-8 sec. This 
estimate will be verified for the model discussed 
in Part II but it apparently should not depend 
sensitively on special assumptions regarding the 
mechanism. 

The limited lifetime would lead to the conclu- 
sion that the barytrons have to be produced as 
secondaries in the atmosphere itself. As rz is so ex- 
tremely small (a barytron of energy 10" ev =100u 
would only have a free path ]/=c-7,-100 cm 
~ 300 m) the barytrons must be created very near 
to the place where they are observed or their 
connection with the 8-decay has to be given up 
altogether. 

On basis of these conceptions the most prob- 
able mechanisms for the actual production of 
barytrons are discussed in the next two sections. 


Section 2. Production of barytrons by protons or 
neutrons 


This mechanism is quite analogous to the 
emission of light quanta by charged particles. 
A proton carrying with it a ‘‘nuclear force field”’ 
can actually emit barytrons when passing with 
sufficient energy through another nuclear force 
field, e.g. of a nucleus at rest. If the theory con- 
verges in the sense discussed in Section 1 the 
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emission of a single barytron should be the most 
probable process and representative of the whole 
effect. 

Because of the resemblance of J(r) to a 
Coulomb field a very rough estimate of the order 
of magnitude of this effect can be obtained by a 
simple “translation” of the ordinary formula of 
radiative emission. The dictionary is: 


Radiation Nuclear force field theory 
a=e/hc=1/137; A=G*/hce=p/M 
=al37u/M~1/10; 
Ro="G/ 


~0.075(u/M)-ro 
~7.5X10-*r0; 
Ro. 


The cross sections of ordinary radiative processes 
for high energies are of the order of magnitude 
Zar, where Z is the atomic number. From this 
we obtain for example for the ratio of the prob- 
ability ¢g for the production of a barytron by a 
proton over the one ¢,, for the production of a 
light quantum by en electron of the same energy 
as the proton, the order of magnitude 


dp X AR? 


— 


Phy ary 


u\i X 
~0.73(—) —<10-*. (3) 
M/ 2 


The factor X is introduced to take care of the 
composite structure of the nucleus and should be 
approximately equal to its mass'* number. The 
ratio of these probabilities therefore is of the 
order of 10-*. 

To show the significance of this result we 
remark that an electron has about an even chance 
to create a quantum of comparable energy in one 
unit length of the radiation theory of showers," 
which is for air ~ 1/30 of the whole atmosphere. 
Therefore a proton will generate approximately 
only 30X10-*~1/30 barytrons when travers- 
ing the whole atmosphere. As the number of 
protons in the cosmic radiation is certainly much 
smaller than the number of barytrons, the above 

8 As the range of the nuclear forces is of the same order 
as the average distance between its constituents one should 


expect an independent superposition of their action with 


regard to our emission effect. 
“ Bhabha and Heitler, Proc. Roy. Soc. A159, 432 (1937); 
Carlson and Oppenheimer, Phys. Rev. 51, 220 (1937). 


emission process cannot explain the occurrence 
of the latter near sea level so long as the anal- 
ogy formula has any meaning. 


Section 3. Production of barytrons by light 
quanta or electrons” 


A production of barytrons by light quanta or 
electrons can only be effected by a simultaneous 
action of the electromagnetic and nuclear-force 
field. The simplest mechanism of this kind is the 
direct photoeffect, i.e. a photon is absorbed by 
a proton and a barytron is reemitted. The process 
is due to the combined effects of the interactions 
photon-barytron and barytron-proton (com- 
pare Section 6). It follows then from dimen- 
sional arguments that the cross section must be 
of the order po 
where F is a homogeneous function of the 
primary energy E and the rest energies of the 
particles concerned. If this function is assumed 
to be of the order unity for high energies'® the 
ratio of the cross sections for producing a 
barytron by a photon over the cross section for 
the production of a normal electron-positron 
pair would be of the order of magnitude 


én” X poRyo X 137m 


2 
) <10-. (4) 


dei” ar, 


This is larger than (3) by roughly a factor 10 and 
as we can expect the presence of many more 
quanta than protons in the high atmosphere it 
might offer a more promising picture. 

On the assumption that the hard component 
of the cosmic radiation is produced as a secondary 
radiation by photons, one can now estimate from 
observational results how large the ratio (4) has 
to be in order that enough hard rays are pro- 
duced. If, in the production process, the whole 
energy of a quantum is given to the barytron, 
the total number of barytrons of energy E 
produced in the atmosphere will be 


N(E)= f x)¥(x)dx, (5) 


% The hypothesis that the hard component is produced 
by the soft component in the atmosphere itself has been 
first put forward by Bowen, Millikan and Neher, Phys. 
Rev. 53, 217 (1938). 

16 As we shall see in the later discussion this assumption 
is subject to grave doubts since, with the present form of 
the theory, the function F generally decreases with 
increasing energy. 
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where ¢,’"(/) is the cross section as a function of 
energy, n(E, x) the number of quanta of energy 
E as a function of depth x below the top of the 
atmosphere, and Y the number of atoms per unit 
volume. If we measure the cross section in units 
of the cross section of normal pair production 
and the depth in the unit length /) of multiplica- 
tion theory'* (4,~0.35 m water equivalent for 
air) dl= ¢.,""Ydx, and (5) reduces to 


N(E)= f $n(E)n(E, dl. (6) 


The distribution n(£, x) of quanta in the atmos- 
phere can be obtained in sufficient approximation 
from the analysis of one of the authors.'? Accord- 
ing to the multiplication theory the number of 
quanta will be slightly larger (by a factor of 
~1.5) than the number of electrons of the same 
energy and shows approximately the same energy 
distribution, which is 


n(E)dE=f(x)E-dE, for E>Eo, (7) 


where f(x) gives the dependence on depth. The 
exponent s lies between 2 and 3, and £, is the 
energy down to which the multiplication of rays 
extends, i.e. ~ 1.5108 ev for air. The form of 
the distribution is nearly independent of altitude 
and therefore f(x) is approximately proportional 
to the total soft intensity in the atmosphere 
which, in turn, gives the greater part of the 
intensity in the atmosphere itself. The number 
of slow electrons below E, should, furthermore, 
be of the same order as those above EF) (i.e. those 
obeying the distribution (7)). Therefore one can 
say that the number of quanta of energy above 
E, at a certain depth will be roughly equal to the 
total number of soft particles, i.e. to nearly the 
total measured intensity, and that their energies 
will roughly follow the same distribution (7). 
Integrating the intensity distribution (Fig. 3 of 
reference 17) at the equator'* over depth one 
obtains in the units used there 600 particle 
meters or, as 1 m water equivalent ~3 unit 

17 Nordheim, Phys. Rev. 51, 1110 (1937) ; 53, 694 (1938). 

18 Only few of the photons produced by the field sensitive 
part of the soft primaries would be energetic enough to 
roduce barytrons which could travel down to sea level. 

his would be consistent with the smallness of the geo- 
magnetic effect at sea level (15 to 20 percent), and it would 


seem therefore quite possible to account for the features of 
this effect on the photon hypothesis. 
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lengths, about 2000 photon lengths for photons 
above an energy of 1.5X10° ev. The number 
Q(£) of photon lengths for quanta of an energy 
larger than E will be therefore 


Q(E) = 2000( 
E>Ey)=1.5X 108 ev. 


The number of hard particles at sea level’? js 
around 7 (in the units used above) and certainly 
more than half of these, i.e. around 4, have 
energies larger than 10° ev. 

We make now the most favorable assumption 
that no barytrons get lost in the atmosphere 
itself by spontaneous disintegration. (This would 
mean, of course, that the theory of B-decay, 
cutlined in Section 1, could not be correct.) Then 
most of the barytrons will come from the upper 
part of the atmosphere where the intensity of 
quanta is largest, and will have already lost an 
energy of 1 to 210° ev by ionization. Therefore, 
only quanta above about 2 to 3X10° ev should 
be counted. With s~2.5 we obtain from (7a) 
Q~40 for E=2X10° ev. 

This means that such a quantum must have 
a chance of about 1 in 10 to create a barytron 
of comparative energy, before it is absorbed 
through the normal production of an electron 
pair. This estimate should hold for any mecha- 
nism. In case the production cross section depends 
on energy, the ratio 1/10 would mean an average 
over all energies, and if the production occurs in 
multiple processes, the cross section would be 
smaller in proportion. But it has to be noted, 
that the actual energy distribution of the hard 
rays as measured by Blackett”? demands that 


19 Compe Table III, reference 17. These figures were 
estimated from experimental data and not deduced from 
the tentative analysis given there. 

20 Blackett, Proc. Roy. Soc. A159, 1 (1937). It is inter- 
esting to remark that the energy distribution for the hard 
rays would be in agreement with the one observed and 
compatible with the absorption curve of the hard rays 
underground if the cross section for barytron production 
is independent of the quantum energy like the cross section 
of electron pair production. In this case the energy dis- 
tribution would be the same as the energy distribution of 
the producing quanta, i.e., of the form (7) which is at least 
in qualitative agreement with observation as shown in 
reference 17. According to this hypothesis of the secondary 
nature of the hard component it is not necessary to assume 
that the hard rays produce hard secondaries in great 
amounts. The main cause for absorption of the hard rays 
will then be again energy loss through ionization and the 
distribution law (7) will give also in this case an absorption 
law E~x* ! (x=depth) as observed. 
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also barytrons of very high energies are pro- 
duced.” 

The correspondence formula (4) gives for the 
ratio of barytron production to electron pair 
roduction by photons for air a value of about 


2x10 (m/M=1850; Z°=50, X~20), which- 


is thus too small by a factor ~50. It seems, 
therefore, not to be possible to give a simple 
account of the barytron production on lines 
similar to radiation theory. 

The discrepancy would still be larger, if the 
lifetime of the barytron is actually as short as 
estimated in sections 1 and 4. The connection 
between the barytrons and radioactive 8-decay is, 
however, not cogent, and it might be changed by 
a reformulation of the 8-theory. 


Part II. A Mopret THEORY FoR 
BARYTRON PRODUCTION 


Section 4. Interaction Hamiltonian and lifetime 
of barytrons 


The considerations of the preceding sections 
were of a purely qualitative nature and based on 
simple analogies to the theory of radiation. It 
remains to be seen whether the characteristic 
properties of the barytron field (i.e. that it 
describes the motion of charged particles of 
finite rest mass having no negative energy 
states) introduce any features which might in- 
validate this analogy. For this reason we give 
in the following sections an actual model theory 
for the effects in question. For this it is necessary 
to have an interaction which can be applied to 
high energy phenomena. This demands that the 
perturbation developments converge at least in 
the same way as those in quantum electro- 
dynamics. The only nuclear force field so far 
proposed which satisfies this requirement is the 
scalar one,” originally introduced by Yukawa. 


"Similar considerations would apply to the purely 
electromagnetic production of positive and negative pairs 
of barytrons by quanta. The cross section for this process 
is, as well known, for high energies, nearly independent of 
energy the ratio to normal electron pair production being 
of the order (m/n)*. As two barytrons are created in a 
single act this ratio had to be of the order 1/20 or x5 m. 
For the probable value of u«~100 to 200 m this kind of 
production can only contribute less than 1 percent of the 
total. The action of electrons compared to photons will 
in both cases (mixed electromagnetic-nuclear field and 
mits ‘ne production) be smaller by a factor 
~ 
# An alternative form of the barytron field of vectorial 
nature and resembling therefore more closely the electro- 
magnetic field, has recently been suggested independently 
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Though this field admittedly does not give satis- 
factory results for the spin dependence of nuclear 
forces, it seems to be worth while to study the 
mechanism of barytron production with a 
definite model. The result is that the orders of 
magnitude come near to the estimates of Part I 
though there are some differences in details. 

For the scalar field of Yukawa the interaction 
with the heavy particles is described by the 
Hamiltonian 


he-Gf (vn pr) 
(8) 


where yy, wp are the wave functions of a heavy 
particle in the neutron or proton state. They are 
assumed to obey Dirac’s wave equation so that 
(yx*Byp), etc., are the simplest scalars. 8 is the 
Dirac operator and ¢ is the scalar wave function 
of the barytron, which is assumed to obey the 
Klein-Gordon equation, and is subject to second 
quantization according to the procedure of 
Pauli and Weisskopf** for Einstein-Bose par- 
ticles. G is an interaction constant with the 
dimension of a charge; the other factors are 
chosen so that the simple form (1) of the poten- 
tial, without any additional constant, is obtained. 

For the quantization of the barytron field 
following the procedure of Pauli-Weisskopf, a 
Fourier decomposition can be used. 


¢= exp ((i/hce)B-r); 


i 
=— bp*); 9 
adpt+b,*) (9) 
. 


by Yukawa (reference 7), Fréhlich, Heitler, and Kemmer, 
Proc. Roy. Soc. A166, 154 (1938); and Bethe, Phys. Rev. 
Abstract (1938). The possibilities for the interaction 
Hamiltonian of the barytron field with heavy particles have 
been studied in detail by Kemmer, Proc. Roy. Soc. A166, 
127 (1938). The interaction (8) used here is the same as 
Kemmer's with his interaction constant g, #0 and all others 
zero. It seems that a formal perturbation treatment of the 
emission processes with his generalized interaction would 
give rise to terms which contain positive powers of the 
initial energy and which, therefore, would become very 
large at high energies. On the other hand, the divergences 
in this theory are far more serious and such formal de- 
velopments would have no physical meaning. It might be, 
however, that this'lack of convergence gives an indication 
of qualitatively new phenomena (multiple processes of 
higher order in the sense of Heisenberg) which might 
finally bring the solution of the barytron problem, though 
at present no method exists for treating these questions. 

*3 Pauli and Weisskopf, Helv. Phys. Acta 7, 709 (1934). 
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Here Ez and Bare the energy and the momentum 
vector (in energy units i.e. B=momentum Xc) 
and the quantities ag etc. stand for the operators 


az*=creation of a barytron of positive charge, 
b,z* =creation of a barytron of negative charge, 
ag=annihilation of a barytron of positive 
charge, 
bg=annihilation of a barytron of negative 
charge. 


From (8) and (9) we obtain the following matrix 
elements for the change of a proton into a 
neutron with momentum N/c and emission of a 
positively charged barytron 


—iGhc(27r)} 
Kp; 
(Ex)! (10) 
P=N+B*; 2’, 


where up and uy are the four component Dirac 
amplitudes of the corresponding states (all wave 
functions are assumed to be normalized per unit 
volume). The matrix element for the reverse 
process of absorption of a barytron with negative 
charge by a neutron is the conjugate complex to 
(10), and for the process N—P+B-, i.e. the 
emission of a negatively charged barytron by a 
neutron, the subscripts P and N in (10) have to 
be exchanged. 

From the matrix elements (10) the interaction 
potential between a proton and a neutron is 
obtained by the normal second order perturba- 
tion formula. For states of motion for which the 
kinetic energies of the heavy particles are small 
compared to yw the 8-brackets in (8) can be 
replaced by unity, and the energies of the inter- 
mediate states are simply the energies of the 
intermediate barytrons. Positive and negative 
barytrons give the same contribution, and one 
obtains therefore for the interaction potential 


exp (—(i/hc)B-r)dB 
V(r) +2 f f f 


G 
=—— exp (—ru/hc)-P, (11) 
r 
where r is the distance between neutron and 
proton, and the factor (27/c)~* represents the 
density of barytron states in momentum space. 
P is an operator which exchanges the proton- 
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neutron character of the heavy particles. Its 
effect on a function antisymmetrical in the 
character variables (as must be assumed for an 
3S state which is symmetrical in all other varj- 
ables) is a change of sign, while the sign is pre- 
served for a symmetrical character function 
(\S state). Equation (11) gives therefore a 
Heisenberg force of space dependence (1), but 
repulsive for the *S state of the deuteron.” *4 

In order to describe the radioactive 8-decay 
an additional interaction between barytrons and 
electrons and neutrinos has to be introduced. 
The simplest possible assumption is, in analogy 
to (8), that 


Bei) 


+conjugate complex}. (12) 


which couples the emission of an electron y, 
and the absorption of a neutrino y, (that is 
emission of an antineutrino) with the production 
of a positive or disappearance of a negative 
barytron and so on. Here g is a new interaction 
constant also of the dimension of a charge. A 
second-order perturbation theory over the stages 
N—P+B-;B- electron + antineutrino gives 
then the formula for the probability of 8-decay as 


4 m 
(7) 


where ¢ is the energy of an emitted electron and 
€) its maximal possible energy (determined by the 
total available energy). 

This formula is identical with the one derived 
from a Fermi interaction of the usual form 


+conj. compl.}, 


(14) 
where Gr is a dimensionless constant and 

2rgG 
( he we 


Gr= 


As Gr is of the order®® 5X10-", we obtain by 
use of (2) 


1 he 
— Gr?~ 10-4, 
he 


Lamb a and Schiff, Phys. Rev. 53, 651 (1938). 

25 Calculated from the data given by Nordheim and Yost, 
Phys. Rev. 51, 992 (1937). For light positron emitters for 
which the heavy particle matrix element is ~1 one has 
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The decay probability of a free barytron (i.e. its 
reciprocal lifetime) at rest becomes on the other 
hand 


je. the lifetime is of the order 10~° seconds in 
agreement with the estimate of Section I. 

As another simple effect for barytrons we 
mention the scattering by a neutron or proton,?® 
ie. the analogy of the Compton effect for light 
quanta. The total cross section integrated over 
all scattering angles for this effect is, with the 
interaction (8), given by 


21+) Eo 
> 
1+2y (1+2y)? M 


where Fy is the initial energy of the barytron. 
For Ey) > M the cross section decreases as 
as in the ordinary Compton effect. For Ey< 
(analogy to Thomson scattering) one obtains a 
factor (.//E)*? to the geometrical cross section 
+R,. This factor is due to the special form (8) of 
our interaction. In ordinary radiation theory 
one has instead of the 8-brackets in (8) the Dirac 
a-operators which are of the order v/c where v 
is the velocity of the charged particles. Expressed 
in the formalism of perturbation theory the 
ordinary Compton scattering is due to virtual 
transitions into negative energy states while with 
the 8-interaction this is not necessarily so. 

As most of the results in this section have 
already been obtained by Yukawa we have 
omitted the detailed derivations. 


Section 5. The radiative emission of barytrons 
by heavy particles 


According to the conservation of charge, 
barytrons can be emitted in collisions of heavy 
particles according to the schemes given in (17). 


I. 

II. 
Ill. Pi+P:>-P+N+Bt 
IV. Nit+N2-N+P+B-. 


(17) 


The cross sections for these processes will be 


** A positive barytron can only be scattered by a neutron 
and a negative one only by a proton. 


equal so that we need only to consider the reac- 
tion (I). With the matrix elements (10), connect- 
ing the change of character of a heavy particle 
with the emission or absorption of a barytron, the 
reaction has to go over two successive inter- 


mediate states and there exist the six possi- 


bilities given in (18). 


a b c 
Po N'+ BY Py N'+ Py > N'+ B*’ 
Not BYP, = Not BYP, 
N'’>P,\+B- N’>P2+B- Not Bt’ 

(18) 
Po N'+ Bt’ Noe P2+B'O 
Pot Pot+B’>N" 
Not 


So long as the heavy particles can be con- 
sidered as free the formulae (17) and (18) are 
directly the equations for the conservation of 
the respective momenta. The dashes in (18) 
denote intermediate states. We consider the 
process from the system of reference in which the 
neutron is initially at rest, i.e. N,=0. According 
to perturbation theory the differential cross 
section will then be 


2x Ep, 
——| I] sos Eso)d ppd pr, (19) 
h Poc 


with > (19a) 
—Ey)(Ey =e Ey) 


where the indices 59, s’, s’’, s denote the initial, 
two successive intermediate, and the final state 
of the total system and dp, and dpp, the density 
functions for the corresponding particles 


BEgdE pdQy P,Ep.dQp, 
(2rhc)* (2rhc)* 


with dQ, and dQp, the elements of the solid 
angles for the barytron and one of the final 
protons, respectively. The factor 6(Es— Eso) 
denotes the conservation of energy. Only those 
values of the momenta will contribute to the 
integral cross section for which the resonance 
denominators in the sum (19a) are small. We in- 


vestigate separately the case A of high energy 
of the oncoming proton, i.e. Ep. > M and case B 
of small energy i.e. Ep.< M. 

A. Relativistic case: Eps>M.—The smallest 
resonance denominators are obtained for small 
values of the momenta B*’ and B~’ for the inter- 
mediate barytron, i.e. for small momenta trans- 
fers between the two heavy particles. If we 
neglect the symmetrization in the final protons, 
i.e. in P; and Pz, and designate with P, the par- 
ticle in which the original neutron Noy goes over, 
then only the series a and e in (18) contain the 
corresponding small resonance denominator and 
need alone to be taken into account.” For these 
terms we have then 


BY 
N’= N”=P,+B-~P). 


(21) 


Since only very small momentum transfers con- 
tribute (see the formulae (30) and (31)), Ep: 
will still be only ~ A/. Therefore we can neglect 
differences Ep:— En» against Ep,’ even in the 
denominators and we have as an expression for 
the conservation of energy 


Ep.=Ep,+Es-. (22) 


Denoting the energies of the heavy particles 
from now on by Fo, £;, E’, E”’ and those of the 
barytrons by and E,;’= we obtain from 
(10), (18), (19a), (21) and (22) 


| (uy 
Eo 


(une 
R” 


27 This procedure is allowed as long as Ep, and Ep, are 
of different order of magnitude, i.e., as long as Es— <Ep,. 
We have checked the correctness of these and other similar 
neglections later on by more detailed calculations, which 
are elementary but tedious and have, therefore, been 
omitted here. 


Ksos = 
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Summing over positive and negative energy 
states of the intermediate neutron and the spins 
of the initial and final states and transforming 
the resonance denominators 


—Po-B— 
=2(Eoe’ —P-B’), 

— (Eo+ = 2(Eoe’ + B’ + 
=2(Eoe’ +Po-B’), 


(24) 


we obtain for the differential cross section (terms 
containing B’ in the nominator are neglected) 


1 

0 


fo 12 (25) 
| Ey’ — (Po-B’)?}? 


This expression contains the dependence on the 
direction of the final particles through the scalar 
product (Po-P;) in the nominator and through 
the intermediate momentum B’ and energy ¢ 
in the denominator. 

For the integration we introduce the following 


variables: the vector 

Q=P.+B (26) 
through its absolute value Q (which determines 
the opening angle @ between P; and B with 
QdQ= P,B sin 6d6) and the polar angles 3 and ¢ 
(the latter is cyclic) with respect to the axis Py, 
and finally the azimuth y of the plane (PB) 
around Q (y is contained only in Py-P;). The 
integration over the azimuths ¢ and y leads to 


f f (M24 


=4n?{ M?+ +0 cos 8} 


and we obtain for (25) 


EoM*{20(M?+ — Po(P2+@2—B?) cos 8} sin ddddQ 
J | Egte? — Po?(Po—Q cos 8)?}? 


and the new variables 


z=1-—cos? and y=Po—Q (29) 


e energy 
the spins 
'sforming 


(24) 


n (terms 


ected) 


1 


(25) 


> on the 
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cos J} 


(27) 


(28) 


(29) 
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we obtain (only the lowest powers in z and y which give the chief contribution to the cross section are 


retained ) 


vo 2 
dy | dz 
4 0 


(30) 
1 w+ 


The limits for y are 


M «—B 
A=P o—(Pi+B) . 
€ 


(31) 
yo= Po—|Pi:—B| =smaller of or 


(relativistic approximations and (22) are used). 
Only the lower limit A is important which sig- 
nifies the minimum possible momentum trans- 
ferred to the initial neutron. Carrying out the 
integration over y we obtain in sufficient approx- 
imation as final result 

Ey—ede 


M 
do= J 
—. (32) 


x 
1+(uEo/eM)? Eo 


This cross section has a maximum as a function 
of the emitted energy for e~(u/M) Eo. 

It is interesting to compare (32) with the 
result one would obtain in the limit of vanishing 
rest mass of the barytron (which would give a 
pure Coulomb law for the proton-neutron inter- 
action). Putting p—-0 in (30) and integrating 


one obtains 
32 Ey (Eo— de 
do,+0~—A Re— (33) 
3 M ee M 


We see that this would give a much larger effect 
than the analogous formula of radiation theory 
and especially a strong divergence for small e. 
This is apparently due to the used 6-interaction 
instead of the a-operators (see the end of section 
4) which increases the cross section enormously. 
The finite rest mass, however, introduces a very 
effective limitation for the decrease of the 
resonance denominators. 

By use of the correspondence picture of 
Weizsaecker-Williams for the description of 
radiative effects we can explain the factors in (32) 
as follows: The increase by the factor (//,)? 


2M* pw dy 1 
~ f 


w+) Ee 


comes from the choice of the 8-interaction. (For 
this reason it would perhaps be more logical to 
introduce G?/u as the critical “‘barytron radius” 
in place of G*/ M as in Section 2.) The decrease 
by the other energy dependent factors is due to 
the screening expressed by the exponential 
factor in (1). 

Integrating (32) over ¢ we obtain for the total 
cross section for the production of a barytron 
of any energy 


and for the average energy loss of a proton (N 
=number of atoms per unit volume and X =the 
effective number of scattering particles in the 
atom, which is of the order®® 3Z, 


—dE/dx~X (35) 


For energies from 10° to 10" ev (34) is of the 
same order of magnitude as would follow from 
the analogy consideration of Section 3. 

B. Nonrelativistic case: Ep.< M.—It is neces- 
sary to investigate this case separately because 
of the large effect of the 8-interaction at small 
energies. In this case (E)<M) the energies of 
all heavy particles concerned will be smaller 
than M and therefore all the B-brackets in the 
matrix elements can be replaced by unity. 
Furthermore all energy differences of the heavy 
particles in the denominators can be neglected 
compared to the barytron energies ¢ and ¢’, but 
all six series of the intermediate states (18) have 
to be taken into account. The evaluation of the 
cross section is otherwise similar to the relativ- 
istic case and leads to the result 


2 T)'(To— ©)! Bde 
(36) 
(2Ty—«)* 


28 According to the scheme (17) a proton can emit a 
barytron ina single way when encountering another proton 
but in two ways in an encounter with a neutron. 
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and integrated 
lg To/n, 


where 7)=P,?/2M is the kinetic energy of the 
oncoming particle. We see that the cross section 
(37) is a smooth continuation of the relativistic 
case (34) with only a slightly different numerical 
factor. 

To obtain an estimate of the importance of 
the emission effect investigated in this section 
we calculate the probability W that a proton 
creates a barytron when passing through matter 
before it is stopped entirely by normal ionization. 
Since the specific energy loss (35) is small com- 
pared to the specific energy loss y by ionization, 
the energy as a function of depth x is 


(37) 


E(x) = yx, 


where E, is the initial energy. We have then with 
the help of (34) 


Eo dE 
W= XN f o(E)— 
M 7 
M 
~16xNARi(—) —, (38) 
7 


if Ey is sufficiently large. Inserting the numerical 
values for air we obtain 


W~(X M/u) X10? ~2 X10 


which as has been discussed in Section 3 is 
insufficient to explain the occurrence of barytrons 
in the cosmic radiation. 


Section 6. Barytron production by photons 
The simplest reactions of this type are 


I. k+P—-N+Bt, 
Il. k+N—-P+B-, 


where k=/hy and k are energy and momentum 
(in energy units) of the light quantum. Since 
both reactions will have the same cross section 
we consider only (I). The process goes in two dif- 
ferent ways over one intermediate state each. 
The first possibility consists in the conversion of 
the light quantum into a_positive-negative 
barytron pair and subsequent reabsorption of 
the negative barytron by the proton. The second 


(39) 
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consists in the emission of the positive barytron 
by the proton and subsequent absorption of the 
light quantum by the barytron. The momentum 
relations are 


k—B++B- P—B++N 
P+B-—--N B’+k—B* (40) 
B’= —B- 


The relations between the energy of the 
emitted barytron and its direction relative to the 
direction of the initial light quantum (angle 9) 
are very similar to those of the original Compton 
effect, i.e. in sufficient approximation (if we 
assume the initial proton to be at rest, i.e. P=0) 


k 
14(k/M)(1—cos 


e+~k; 


et 


(41) 
kR<M. 


The matrix elements for pair production and 
absorption are according to Pauli and Weisskopf* 


r\! ex: (B+—B-) 
KHpair= 
2k (ere)! 


\'e,-(Bt+B’) 
2k 


(ere’)! 


where «+ and e«~=€’ denote the energies of the 
barytrons and e, the unit vector of polarization 
of the quantum (e, is orthogonal to k). The 
matrix element for the second step is again (10). 
We have then 


e.:(B+—B~-) =e,-(Bt+B’) =2e,-Bt 


and the cross section becomes 


(43) 
2e-(uy*Bup) |? 
= Ge*5( Es — Es») 
ke 


2° Pauli and Weisskopf, reference 23 formulae (49) and 
(53). It is interesting to note that in (42) the difference 
B+—B- enters (in reference 23 the definition of one of the 
momenta is reversed). This clearly has to be so as the 
resultant electric current of the two barytrons is propor- 
tional to this difference. 
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We treat again separately the cases of high and 
low energy. 

A. Relativistic case: k >p.—We sum over the 
spin states of the heavy particle and use the 
following relations coming from the conservation 
of energy and momentum (compare (39) and 


(40)) 
§(Es— Es.)dQ 


det EyB Mdet+ 
=2r——| sin ddd=2r 
dE Mk 


et? 


Then we obtain for the differential cross section 
with Ey = M+k— et 
2M?(k— et) Bidet 


. (46) 
MR? 


The total cross section will be then (putting 
e+ =e) 


[2M+(k—©) |(k—©ede 
o=4rRopo 
+2M(k—0) 


The main contribution to this integral comes 
from the region kR—p<e<k which means that 
practically always the whole energy is given to 
the barytron. Carrying out the integration over 
ewe have with sufficient approximation 


TRopoF(k), 
(48) 


2 


9 k 


uM 
ra) (1+ -+—}. 
k? 4 2M 


This expression contains the factor Ropo as dis- 
cussed in Section 3 but the dimensionless factor 
F(k) = F(hv) decreases strongly for high energies 
similar to the ordinary Compton effect, which 


again means that the effect cannot account for 
the number of observed barytrons. 

B. Nonrelativistic case: k~p.—Here the 
6-brackets in (10) are unity, 


det 
dE 


et ~k, 


and the integral cross section becomes after some 
calculation 


o=2rRopo- — lg ——-2 


B k-B 


k+B 


(49) 


This is of the same order as one would expect 
from the extrapolation of (48) to low energies. 
For k~u, (49) contains a factor M/y in addition 
to the dimensional expression Ropo. This cross 
section would actually approach the magnitude 
required by observation (comp. Section 3) if it 
were independent of energy. 

The authors wish to express their thanks to 
Professor E. Teller for stimulating discussions 
on the subject. 


Note added in proof: A more exact evaluation of the 
fundamental constant G of the theory has recently been 
given by Sachs and Goeppert-Mayer (Phys. Rev. 53, 991 
(1938)). They obtain a value ~0.3 for G*/he instead of 0.1 
as assumed in this paper. The cross sections of the pro- 
duction processes are thereby increased very effectively 
(especially for the process discussed in Section 3) so that 
they approach already those required by observation. 
Furthermore, the convergence of the perturbation cal- 
culations becomes doubtful, even for the scalar barytron 
field. This means that the energy dependence of the cross 
sections cannot be determined any more. Even if the pos- 
sibility for a better estimate of the production probabilities 
seems thus to be still more remote, the prospect for 
establishing the discussed connection between nuclear 
forces and the hard component of the cosmic radiation is 
certainly much better now. 


= 
= 
(41) 
on and 
= | 
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Measurements were made of the transmission of liquid 
para- and normal hydrogen for the slow neutrons obtained 
from a Rn-Be source by scattering in paraffin. In some of 
the experiments the paraffin was at room temperatures 
and in others it was cooled with liquid air. Teller, and 
Schwinger and Teller, have shown how the spin de- 
pendence of the interaction between neutrons and protons 
can be investigated by such measurements. As the neutrons 
captured by the protons were less than 1 percent of the 
scattered neutrons, the scattering by the thin layers of 
liquid hydrogen used was determined from the trans- 
mission measurements. The scattering cross section of 
liquid hydrogen for slow neutrons was found to depend 
upon the velocity of the neutrons and on the ortho-para 
composition of the hydrogen. The following values were 
obtained for the neutron scattering cross sections, o, of 
liquid ortho- and parahydrogen: for ~300°K neutrons 
(paraffin at room temperatures), o (ortho) =56X 10~** cm? 
per molecule, (para) =29X 10°"; for ~120°K neutrons 


(paraffin cooled with liquid air), (ortho) =79 and 
o (para) =18X10-*. We were able to prove conclusively; 
(1) that the interaction between neutrons and protons js 
dependent upon the relative alignment of their spins, and 
(2) that the energy of the singlet state of the deuteron, in’ 
which the spins of the proton and neutron are antiparallel, 
is greater than the energy of these particles when far apart, 
that is, the singlet state of the deuteron is virtual. The 
scattering+capture cross sections of liquid oxygen and 
liquid nitrogen were measured for ~300°K neutrons, and 
of water for ~300° and ~120°K neutrons. The cross sec- 
tions for the 300°K neutrons were: o(O)=4X10-* cm? 
per atom, o(N)=13X10-* cm?* per atom and o(H,O)= 
91.2 10-** cm? per molecule. For the ~120°K neutrons 
o(H,O) = 116 10-** cm? per molecule. The scattering cross 
sections of the protons in water, (1/2)[o(HO) —o(0)], 
were 43.6 X 10% cm? per proton for ~300°K neutrons and 
56 X 10-* cm? per proton for ~120°K neutrons. 


INTRODUCTION 


HE abnormally large scattering cross section 
of hydrogen for slow neutrons led to the 
assumption that the interaction between neutron 
and proton depends on the spins of these par- 
ticles. This hypothesis, attributed to Wigner by 
Feenberg and Knipp,'! had been generally ac- 
cepted although it lacked direct experimental 
verification.” Teller’ pointed out that the exist- 
ence or nonexistence of a spin dependent neutron- 
proton interaction could be determined from a 
study of the scattering of slow neutrons by para- 
and orthohydrogen. He showed that if the inter- 
action proved to be spin dependent it should be 
possible to tell from such experiments whether 
the singlet level of the deuteron is stable or 
virtual. In the detailed analyses subsequently 
* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
** National Bureau of Standards, Washington, D. C. 
¢ Columbia University, New York, N. Y. 


tt Columbia University, New York, N. Y. Now at the 
University of Illinois. 
anne Feenberg and J. K. Knipp, Phys. Rev. 48, 906 

2 See reference 4 for a discussion of the evidence to be 
obtained from Amaldi and Fermi’s measurements of the 
cross section for radiative capture. 

3E. Teller, Phys. Rev. 49, 420 (1936). 


carried out by Schwinger and Teller* it was 
shown that the effects were large enough to be 
readily detected experimentally, and that it 
might also be possible to obtain information con- 
cerning the range of the neutron-proton forces. 
Schwinger® has since pointed out that such 
experiments would also serve to confirm the 
choice of 3 rather than 3 for the spin of the 
neutron, a choice hitherto made entirely because 
of its simplicity. 

In discussing the para- and orthohydrogen 
scattering cross sections for slow neutrons, briefly 
referred to as the para and ortho cross sections, 
three possibilities are to be considered : 


(1) No dependence of forces on spin 

(2) Spin dependent forces with 
(a) a real (stable) singlet level of the deuteron 
(b) a virtual singlet level of the deuteron. 


In case (1) the only difference to be expected 
in the para and ortho scattering cross sections 
arises from the different spacing of their rota- 
tional energy levels. At room temperature only 
the first few levels are occupied, and for neutrons 
of energy comparable to the excitation energies 


4 J. Schwinger and E. Teller, Phys. Rev. 52, 286 (1937). 
5 Julian Schwinger, Phys. Rev. 52, 1250 (1937). 
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of these levels parahydrogen should have the 
larger cross section, since the closer spacing of its 
levels makes inelastic scattering more probable. 
The difference should be small however and 
should be negligible for neutrons either above or 
below this energy range. In the absence of spin 
forces all scattering must be of the para—para or 
ortho—ortho type, that is the rotational quantum 
number of the scattering molecule must change 
by 0, +2, +4, 

In cases (2a) and (2b) the presence of spin 
forces makes possible in addition para—ortho and 
ortho— para scattering (AJ= +1, +3, ---). If, as 
is the case in liquid hydrogen, the scattering 
molecules are in their lowest rotational levels 
(J=0 for para, J=1 for ortho) spin forces 
should cause the para and ortho cross sections to 
differ markedly for neutrons of energy insufficient 
to excite the 0—1 transition. For neutrons of 
this energy no inelastic scattering is possible in 
parahydrogen, while in orthohydrogen the scat- 
tering process involving the 1-0 transition 
remains. 

For slow neutrons whose wave-length is com- 
parable to the distance between the protons in a 
hydrogen molecule interference effects are im- 
portant, and it is chiefly in these that cases 
(2a) and (2b) differ. The phase of a neutron 
wave scattered by a proton depends on the sign 
of the binding energy involved, and changes by 
approximately a when the binding energy 
changes sign. Hence in case (2a) the scattered 
waves corresponding to singlet and triplet inter- 
action would have like phases, while in case (2b) 
they would have opposite phases. Case (2b) 
brings a possibility of obtaining destructive 
interference in parahydrogen with its antiparallel 
proton spins, for here interactions of the singlet 
and triplet type should occur simultaneously. 
The scattered waves should interfere destruc- 
tively, with the result that for low energy neu- 
trons the scattering cross section of parahydrogen 
should be much smaller than that of ortho- 
hydrogen, where this type of interference cannot 
occur. However, simple calculations based on the 
idea that the spin of a neutron impinging on a 
parahydrogen molecule is parallel to one proton 
spin and antiparallel to the other are not legiti- 
mate, although the trends they indicate are in 
the right direction. Schwinger and Teller have 


considered the problem more carefully and have 
given the necessary detailed calculations based 
on the use of appropriate wave functions. 
Experiments to test their predictions, under- 
taken jointly by Columbia University and the 


‘National Bureau of Standards, have been carried 


out in the Bureau’s Low Temperature Labora- 
tory. Preliminary reports®:’ have already been 
published. 

Considerations of intensity and geometrical 
arrangement render it more feasible to perform 
transmission than scattering experiments. Since 
the neutron-proton capture cross section®: ° is of 
the order of 0.3 less than 1 percent 
of the average scattering cross section, the 
difference between the scattering+ capture cross 
section given by transmission experiments and 
the true scattering cross section is inappreciable. 
Good geometrical conditions are exceedingly 
important in transmission experiments. Experi- 
ments under other conditions, employing non- 
parallel neutron beams, and with both source and 
detector close to the scattering sample have been 
shown to be especially unreliable with scattering 
samples containing hydrogen. The uncertain 
angular distribution of the neutrons, together 
with large corrections for multiple scattering, 
slowing down processes, back scattering, etc., 
make such experiments very difficult to interpret. 
The design of the present apparatus was such as 
to make the errors introduced by geometrical 
factors negligible. 

Since the apparatus was well adapted for use 
with other liquefied gases, the transmissions of 
liquid O2. and Ne were also measured. As a useful 
check on the accuracy of the work, several 
measurements were also made with water in the 
apparatus. 


APPARATUS AND PROCEDURE 


The apparatus divides naturally into two 
parts: one the equipment for preparing and 
analyzing liquid hydrogen and maintaining it in 
the path of the neutron beam, and the other the 
source, detector, and amplifier system. The 

6 J. R. Dunning, J. H. Manley, H. J. Hoge, F. G. Brick- 
wedde, Phys. Rev. 52, 1076 (1937). 

7J. R. Dunning, H. J. Hoge, J. H. Manley, and F. G. 
Brickwedde, Phys. Rev. 53, 205 (1938). 


8 E, Amaldi and E. Fermi, Phys. Rev. 50, 899 (1936). 
9G. A. Fink, Phys. Rev. 50, 738 (1936). 


Fic. 1. Arrangement of source and detector relative to 
scattering sample. Scattering material was condensed in 
one or both of the two cells in the lower vacuum chamber 
of the special Dewar flask. 


arrangement of this second part of the apparatus 
is shown in Fig. 1. The sample was contained in 
either or both of two cerium glass cells which 
formed part of a specially designed Dewar flask. 
The Rn-Be source of approximately 500 milli- 
curies was contained in a “howitzer” in which the 
neutrons were slowed down. This howitzer" con- 
sisted of a mass of paraffin in the shape of a deep, 
very thick walled cup contained in a Dewar 
flask. Copper tubing for the circulation of liquid 
air ‘was imbedded in the paraffin, and thermo- 
couples were provided at two different points for 
temperature measurement. Immediately in front 
of the source was a solid copper cylinder to 
shield the detector from gamma-rays. The source 
was placed in a small source Dewar (2 cm 
diameter) containing a thermocouple and a 
heater. These were used when the howitzer was 
cooled with liquid air, to avoid any change in 
efficiency of the source due to localized con- 
densation of the Rn. When using liquid air it 
was also necessary to heat the outer copper baffle 
at the mouth of the howitzer, to prevent the 
collection of moisture or frost. The effectiveness 
of the howitzer in cooling the neutrons will be 
discussed later in the paper. 

The neutron beam was well diaphragmed and 
shielded by cadmium as shown in Fig. 1. The 
opening in the cadmium diaphragm placed 
between the two cells was 3.5 cm in diameter, 
1 cm less than the diameter of the cells. The 
detector was an ionization chamber 3 cm in 
diameter and 13 cm in length, filled to a pressure 
of 3 atmospheres with BF; gas. The average 
solid angle subtended by the chamber at the 


10 P, N. Powers, H. Carroll, J. R. Dunning, Phys. Rev. 
51, 1112 (1937). 
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cells was less than 0.02 steradian. The chamber 
was connected to a linear amplifier whose output 
was recorded by a scale-of-two thyratron counter. 

All glass in the path of the neutron beam was a 
special cerium glass, whose exceptionally low 
boron content results in very small neutron 
absorption. Total glass in the path of the beam 
was about 6 mm, the cell walls being slightly 
less than a millimeter thick and the Dewar walls 
slightly more. Total absorption by empty cells, 
Dewar walls, and silvered surfaces was about 
6 percent for room temperature (~300°K) 
neutrons. Cell volumes were determined by 
weighing mercury and diameters by direct meas- 
urement. These data gave the mean internal cell 
thicknesses, which appeared very uniform over 
the entire area of a cell. In any case only a 
difference in the mean thickness of the exposed 
(60 percent) and shielded (40 percent) area of a 
cell would cause appreciable error. Transmission 
measurements with water in the cells were used 
to check the cell thicknesses. Satisfactory agree- 
ment was found in all cases except one. In this 
case (0.234 cm cell) the transmission data were 
used to determine the thickness of the cell, 
resulting in a small reduction of the value 
originally found. 

Figure 2 is a diagram of the apparatus used in 
the preparation and analysis of the liquid hydro- 
gen samples. Glass tubes leading to the cells 
passed through liquid hydrogen in the upper 
part of the Dewar. Hydrogen of the desired 
composition was condensed in these tubes and 
ran down into the cells. This arrangement per- 
mitted the amount of hydrogen in the cells to 
be kept constant during the course of the 
measurements, for any hydrogen lost by evapora- 
tion was simply recondensed and returned to the 
cells. Since the level of the liquid in the cells 
could not be observed, the volume of gas con- 
densed was measured. This was done either by 
taking the hydrogen from a cylinder or tank of 
known volume and measuring the fall in pressure, 
or by using a Toepler pump of known capacity. 
A considerable margin of safety was allowed, 
sufficient to bring the level of the liquid well up 
into the connecting tubes. 

Hydrogen to be condensed in the cells was 
taken from a supply tank of especially prepared 
electrolytic hydrogen, the purity of which was 
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believed to be in excess of 99.98 mole percent. 
This hydrogen was passed over hot copper 
(600°C) and through a liquid-air trap to remove 
oxygen, which is an active ortho-para conversion 
catalyst in liquid hydrogen." Passage over hot 
copper also insured attainment of the high 
temperature equilibrium concentration of 75 
percent ortho 25 percent para (normal hydrogen). 
When normal hydrogen was desired in a cell, 
this gas was condensed directly into it. When 
parahydrogen was desired, the normal hydrogen 
was first condensed in a tube containing a con- 
version catalyst in granular form. The catalyzing 
chamber was kept at 20.4°K by surrounding it 
with liquid hydrogen. At this temperature the 
equilibrium concentration is 99.8 percent para- 
hydrogen.” 

Hydrogen from the catalyzing chamber was 
either condensed directly in the cells or stored 
until needed in a tank which could be filled at 
liquid-air temperature and emptied at room 
temperature. 

A promoted iron catalyst" of the kind used in 
the synthesis of ammonia was placed at our 
disposal by Dr. P. H. Emmett. The catalyst, 
received in the form of an oxide, was prepared 
for use by reducing in hydrogen at 500°C, out- 
gassing while hot and cooling to: 20.4°K before 
admitting hydrogen again. Liquid hydrogen 
condensed on the catalyst was allowed to remain 
in contact with it 15 minutes or longer before use. 


$08 


Fic. 2. Apparatus used in the preparation and analysis of 
the liquid hydrogen samples. 


"R. B. Scott, F. G. Brickwedde, H. C. Urey, and M. H. 
Wahl, J. Chem. Phys. 2, 454 (1934). 

2 R. W. Harkness and W. E. Deming, J. Am. Chem. 
Soc. 54, 2850 (1932). 

®8P. H. Emmett and R. W. Harkness, J. Am. Chem. 
Soc. 57, 1624 (1935). 


The ortho-para composition of the contents of 
any cell was determined immediately after trans- 
mission measurements on it had been completed. 
Brickwedde and Scott'* have measured the vapor 
pressures of normal hydrogen, parahydrogen, and 


mixtures of several intermediate concentrations, 


using a differential oil manometer to measure 
the pressure differences. With the aid of their 
results the composition of an unknown sample of 
hydrogen may be determined by comparing its 
vapor pressure with that of normal hydrogen. 
In the neighborhood of their boiling points the 
vapor pressures of para- and normal hydrogen 
differ by about 26 mm Hg; that of para is higher. 
This is approximately 1 mm fall in vapor pressure 
for each 3 percent increase in ortho content. In 
making an analysis about one-fourth of the 
contents of a cell was removed; then a sample 
from the cell and a sample of freshly prepared 
normal hydrogen were condensed in separate 
compartments in the vapor pressure comparison 
block. Their vapor pressure difference was then 
read on the differential mercury manometer, the 
total pressure on the normal hydrogen being 
given by the absolute manometer. The copper 
comparison block was connected to the glass 
system by 3 mm copper-nickel tubes. 

When parahydrogen (99.8 percent) is con- 
densed in a cell there is no change in composition 
with time since the material is in equilibrium. 
With normal hydrogen this is not true, for even 
with no catalyst present the ortho content de- 
creases at the rate of about 1 percent per hour." 
When analysis of a sample from a cell filled with 
normal hydrogen showed a composition of, for 
example, 72 percent orthohydrogen, it was 
assumed that the composition had changed 
linearly with time from 75 percent at the time 
of condensation until the time of analysis. The 
compositions given in Table I correspond to the 
mean time of transmission measurement. Strictly 
the change with time is not quite linear since it 
depends on collisions between two orthohydrogen 
molecules. Because the frequency of this type of 
collision is proportional to the square of the 
orthohydrogen concentration the corrections 
were negligible for all ‘para’ samples. 

It was necessary to reduce the pressure above 
the liquid hydrogen bath until its temperature 


4 To be published. 
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fell almost to the triple point before liquid 
hydrogen would collect in the cells. With the 
bath at this temperature and the gas pressure in 
the cell slightly above one atmosphere rapid 
condensation would begin after 15 or 20 minutes. 
Liquid was kept in the cells under a pressute of 
200-600 mm Hg, the mean temperature of the 
samples being taken as 18.4°K. To remove liquid 
a small heating coil of No. 38 constantan wire 
(50 ohms) was placed in the bottom of each cell, 
below the edge of the cadmium diaphragm. 


MEASUREMENTS 


Each transmission value was computed from 
the results of from 3 to 6 ten minute counts. 
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Roughly one-third of the counts were taken with 
a cadmium sheet 0.08 cm thick placed just jp 
front of the howitzer, so that the effect of the 
high energy cadmium-penetrating neutrons could 
be determined. After a set of counts a cell was 
either filled or emptied and another set taken. 
The normal order was: (1) both cells empty, 
(2) first (thicker) cell full, (3) both cells full, 
(4) second cell full (first empty), (5) both cells 
empty. With both cells empty the approximate 
number of neutrons counted in 10 minutes was 
2500, of which all but approximately 18 percent 
were cut off by interposing the cadmium sheet, 
After a series of cell fillings with normal hydrogen 
a similar series was taken with parahydrogen, 


TABLE I. Transmissions and cross sections (scattering+capture) of para- and normal hydrogen for neutrons of four different 
energy distributions. T is the neutron temperature; C is the composition in percent of orthohydrogen of the scattering medium; 
x is the thickness; I/Io, the transmission; o, the cross section per molecule. 


CROSS SECTIONS OF 
—In I/To a(obs.) PuRE COMPONENTS 
> % o—He x I/lo =nox a(obs.) —o(calc.)§ | (By EXTRAPOLATION) 
cm cm?} cm? cm? 
~300°K 74.2 0.234 0.780 0.249 48.6 —0.2 Forthy =55.8 
72.5 0.32 0.705 0.351 50.1 +1.7 
73.9 0.35 0.695 0.363 47.4 —1.4 
74.0 0.58 0.538 0.620 48.7 —0.1 
73.3 0.93 0.371 0.991 48.6 0.0 
1.7 0.234 0.871 0.139 27.1 —2.3 S para = 29.0 
1.4 0.32 0.809 0.212 30.1 +0.7 
3.0 0.58 0.671 0.40 31.4 +1.6 
2.5t 0.814 0.593 0.521 29.3 (—0.4) 
~120°K 74.5 0.234 0.720 0.329 64.2 +1.1 Gortho = 78.7 
73.8t 0.32 0.656 0.422 60.1 (—2.6) 
72.6 0.58 0.449 0.80 63.0 +1.1 
74.0 0.814 0.34 1.079 60.6 —2.2 
1.4 0.234 0.911 0.092 18.0 —0.4 para = 17.6 
4.1 0.32 0.879 0.128 18.3 —1.8 
6.5 0.58 0.740 0.301 23.6 +2.0 
5.0 0.814 0.690 0.371 20.8 +0.2 
5.9 0.58 
743) 25.6 0.234 0.814 | 0.50* 0.693 38.9 +5.7 
3. .58 
}s7.2 0.58 0.814 | 0.39* 0.941 | 52.7 | +02 
~120°K filtered through { 74.2 0.58 } — 
0.58 cm n—Hy, 73.4 0.234 0.769t 0.263 51.1 
16.2 0.234 0.878t 0.130 | 25.4 — 
~120°K filtered through 5.9 0.58 \ +“ 100.1 
0.58 cm p— He \ 74.3 0.234 0.671+ 0.399 77.6 Certhe 
{$4 0.58 \ =12.6 
1.4 0.234 0.932t 0.070 13.8 


* Transmission of both cells considered as a unit, one containing p— He, the other n— He. 


Transmission of second cell for filtered beam from first cell. 
Estimated. 
Deviations from linear equations of Fig. 3. 
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TABLE Il. Transmissions of liquid O2 and and cross 
sections (scattering +-capture) per atom for ~300°K neutrons. 
The thickness is designated by x. I/Io is the transmission and 
¢ is the cross section. 


—In 
MATERIAL x AVERAGES 
cm 107*4 10774 cm? 
Liquid | 0.35 | —— | —— o(0) =4.05 
Oz 0.58 0.905 0.995 3.98 
0.93 0.847 0.166 4.12 
Liquid 0.35 0.870 0.140 11.5 o(N)=12.7 
N, | 0.58 | 0.761-| 0.273 | 13.5 
0.93 0.657 0.420 13.0 


all with the howitzer at room temperature. 
Then the howitzer was cooled with liquid air to 
approximately 100°K and the process repeated. 
In this condition the transmission was also 
measured with parahydrogen in the first cell, 
normal hydrogen in the second, and vice versa. 
Data on the cadmium penetrating background 
have not been included in the tables. This back- 
ground has in all cases been subtracted to give 
the transmission values reported. Data for hydro- 
gen are given in Table I. In calculating the cross 
sections, n the number of molecules per cm* was 
taken from the results of Scott and Brickwedde."® 
It is interesting to note that because of the large 
thermal expansion of liquid hydrogen this number 
is approximately 3 percent greater at 18.4°K than 
at 20.4°K, the boiling point of normal hydrogen. 

The procedure with liquid Ne, and 
was similar to that with hydrogen except that 
no measurements on Oz and Ne were made while 
the howitzer was cooled with liquid air. Data for 
liquid O2 and N: are given in Table II; for water 
in Table III. 


DISCUSSION OF RESULTS 


The dependence of cross section on composi- 
tion and on neutron energy is shown in Fig. 3, 
where cross section per hydrogen molecule has 
been plotted against percent orthohydrogen in 
the scattering layer. Each of the four :straight 
lines refers to neutrons of a different energy 
distribution. The equations of the lines were 
determined simply to make the sum of the 
deviations zero for the para group of points and 
likewise zero for the points taken with normal 
hydrogen. The intercepts of these lines on the 


®R. B. Scott and F. G. Brickwedde, J. Chem. Phys. 5, 
736 (1937). 


TABLE III. Transmission of water and the cross sections 
(scattering+capture) of (1) water per molecule and of (2) 
the protons of water per proton, for neutrons of ~300°K, 
and ~ 120° K. T is the neutron temperature; x is the thick- 
ness; is the transmission; }(o(H2O) —a(0)). Average 
for “room temperature” neutrons ¢p=42.2 X The 
0.93 cm thickness was not used in obtaining this value. 


—In 
T x I/lo =nox op 
cm cm? | cm? 

~300°K | 0.23, | 0.495 0.703 
~270°K | 0.32 | 0.370 0.994 92.0 44.0 
~300°K | 0.35 | 0.349 1.05 89.1 42.5 

- 0.58 | 0.190 1.66 85.0 40.5 

sty 0.58 | 0.170 1.77 90.5 43.2 

se 0.93 | 0.079 2.54 81.2 38.6 
~120°K | 0.23, | 0.400 0.915 116.0 56.0 


0 and 100 percent axes give the cross sections of 
pure para- and pure orthohydrogen for neutrons 
of each energy distribution. The values of these 
intercepts are given in the last column of Table I. 
The use of linear extrapolation requires some 
justification. For very thin layers a linear relation 
between cross section and composition is to be 
expected, since this requires only that the cross 


. section of an individual molecule be independent 


of whether its neighbors are para- or ortho- 
molecules. For thick layers a nonlinear relation 
is to be expected because of multiple scattering. 
The geometry of the present apparatus was such 
as to minimize any complications due to multiple 
scattering, since a neutron was effectively re- 
moved from the beam by its first collision, its 
subsequent history being of no importance. The 
transmission curves for water (Fig. 4) which has 
approximately the same proton density as liquid 
hydrogen showed a reasonably smooth expo- 
nential decrease which is experimental evidence 
that multiple scattering was not significant. 

The energy distribution of the neutrons is the 
most uncertain factor in the interpretation of the 
results. This depends on the effectiveness of 
the howitzer in cooling down the neutrons from 
their high energy and establishing something of 
the nature of a Maxwellian distribution. A 
howitzer similar to the one used in these experi- 
ments had previously been tested with a velocity 
selector.'®!* With the howitzer at room tempera- 

%*J. R. Dunning, G. B. Pegram, G. A. Fink, D. P. 
Mitchell, E. Segré, Phys. Rev. 48, 704 (1935). 

7G. A. Fink, J. R. Dunning, G. B. Pegram, D. P. 


Mitchell, Phys. Rev. 49, 103 (1936). 
8G. A. Fink, Phys. Rev. 50, 738 (1936). 


272 BRICKWEDDE, DUNNING, HOGE AND MANLEY 


NEUTRON TEMPERATURE 
© 300°K 
@ 120° K 
& 120° K, FILTERED 
«120° K, p-H, FILTERED 


On. 
| 


, CM? x 


50 
a 
re) 
oO 
uJ 
"95 

0 | | 
(e) 20 40 


60 100 


PERCENT ORTHOHYDROGEN 


Fic. 3. Dependence of neutron cross section (scattering + capture) per molecule on ortho-para 
composition of the liquid hydrogen scattering sample. Each straight line refers to neutrons of a 
different energy distribution. The indicated temperatures are only approximate. 


ture the emergent neutrons absorbable in cad- 
mium were found to possess a distribution whose 
maximum corresponded approximately to ther- 
mal equilibrium with the paraffin. However, it 
was shown that more high energy neutrons were 
present than in a true Maxwellian distribution. 
The effect of all high energy neutrons above the 
cadmium absorption limit (about 0.4 ev!®: ”) is 
eliminated by the cadmium difference method, 
but there still remains a long high energy “‘tail’’ 
extending up to about 0.4 ev. 

The obtaining of a beam of neutrons with a 
Maxwellian distribution of velocities charac- 
teristic of liquid-air temperatures is of greater 
difficulty than for room 2 
The cooling efficiency of the howitzer used in 
these experiments was investigated by Mr. P. N. 

190. R. Frisch and G. Placzek, Nature 137, 357 (1936). 

20 J. G. Hoffman and M. S. Livingston, Phys. Rev. 52, 
1228 (1937). 


21Q. R. Frisch, H. von Halban, Jr., Jorgen Koch, Nature 


139, 922 (1937). 
22 P. B. Moon and J. R. Tillman, Proc. Roy. Soc. A153, 


476 (1936). 


Powers of Columbia University, who measured 
the absorption of the emerging neutrons in boron. 
Using the cadmium difference method, he found 
that the absorption cross section of the neutrons 
in boron was 1.6 times greater when the howitzer 
was cooled with liquid air than when it was at 
room temperature. 

Experiments have indicated that absorption 
in boron follows the 1/z law®*: *4 hence, if 7; and 
T, are the “effective’’ temperatures of the 
neutron beams obtained when the howitzer is at 
room temperature and when it is cooled with 
liquid air, respectively, we have 


(T;/T2)'=1.6. (1) 


The two imbedded thermocouples indicated that 
the average temperature of the paraffin was 
approximately 105°K while liquid air was being 
circulated. In the ideal case in which the neutron 


23F, Rasetti, E. Segré, G. Fink, J. R. Dunning, G. B. 


Pegram, Phys. Rev. 49, 104 (1936). 
*%F. Rasetti, D. P. Mitchell, G. A. Fink, G. B. Pegram, 
Phys. Rev. 49, 777 (1936). 
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beams would be in thermal equilibrium with the 
paraffin at either temperature and the BF; 
receiver counted all neutrons passing through it, 
this ratio would be (300°/105°)!=1.69. Eq. (1) 
cannot be used to calculate 7; or 72 unless some 
assumption is made concerning the relative 
effectiveness of paraffin at 300°K and at 105°K 
in reducing the effective temperatures of the 
neutrons. The effective temperatures must be 
higher than the temperature of the paraffin. 
However, it does not seem probable that they 
are very much higher since with Eq. (1) this 
would lead to a much larger difference in neutron 
temperatures than the change in the temperature 
of the paraffin. Since 7; cannot be less than 300°, 
T; cannot be less than 117° which corresponds to 
a mean energy of 0.015 ev. This difference in 
temperature will make a considerable difference 
in the fraction of the neutrons having energies 
greater than 0.023 ev, which is the minimum 
energy required to excite the O- 1 transition, 
giving rise to inelastic scattering in p— He. We are 
therefore referring to the neutrons as ~300°K 
and ~120°K neutrons to emphasize the fact 
that the effective temperature of the cold neu- 
trons is above the temperature of the howitzer. 

Since the absorption of neutrons in boron 
depends upon their velocity, the fraction of the 
neutrons entering the BF; receiver that are 
counted depends upon the velocity distribution 
of the neutrons. The receiver used in this investi- 
gation would absorb and count about 40 percent 
of an ~300°K neutron energy distribution. In 
accordance with the 1/v law a larger percentage 
of the slower neutrons than of the faster ones 
should be counted. Since in each transmission 
measurement counts are taken both with and 
without the scattering sample in the beam, the 
detector will introduce an error only if the energy 
distribution of the neutrons is altered by passage 
through the scattering sample. Since the scatter- 
ing is to some extent selective this will necessarily 
be the case. The calculation of this correction is 
difficult even when the velocity distribution of 
the neutrons is known. In our experiments the 
velocity distribution was not exactly known and 
no correction of this kind has been applied to the 
measurements. The correction, however, is be- 
lieved to be small. 

This uncertainty in energy distribution is un- 


fortunately characteristic of experiments involv- 
ing slow neutrons, and is almost always the chief 
limiting factor in the interpretation of results. 
It seems probable that the possibility of using 
parahydrogen as a neutron filter brought out by 
the present experiments may help to clear up 
this situation. The use of two separate cells for 
liquid hydrogen (suggested by Professor Teller) 
made it possible to study the effect of filtration. 
By placing parahydrogen in the first cell, a 
filtered beam was obtained which could then be 
investigated by placing scattering material in the 
second cell. As a matter of interest measurements 
were also made using a normal hydrogen filter, 
so that in all four different neutron energy 
groups were investigated. These are designated 
as ~300°K, ~120°K, ~120°K para-filtered, and 
~120°K normal-filtered. 

It is evident from Fig. 3 that the scattering 
cross section is sensitive both to hydrogen com- 
position and to neutron temperature. The de- 
pendence of cross section on composition in- 
creases with decreasing neutron energy. The 
slope of the ~300°K neutron curve is least, that 
of the ~120°K neutron curve is greater, while 
the slope of the para-filtered neutron curve is 
still greater. To see that this last result fits into 
the picture it should be remembered that accord- 
ing to Schwinger and Teller’s theory para- 
hydrogen should be practically transparent for 
neutrons of less than 0.023 ev. Hence the effect 
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Fic. 4. Transmission of neutrons by water, liquid O, 
and liquid Ne. Neutron temperature ~300°K except as 
noted. 
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of para-filtration should be to remove high 
energy neutrons from the ~120°K beam and 
effectively cool it still further. The curve for 
normal-hydrogen-filtered neutrons is less inter- 
esting and more difficult to interpret. The 25 
percent of parahydrogen in the filtering layer 
would tend to remove high energy neutrons from 
the beam, while the 75 percent of orthohydrogen 
would work somewhat in the opposite direction. 
Less weight is to be given to the dotted curves 
than to those for ~300°K and ~120°K neutrons, 
both because of the fewer experimental points 
involved, and because of the larger number of 
measurements required to determine a single 
point. 

Figure 3 shows that the para cross section is 
consistently smaller than the ortho cross section. 
This one fact is sufficient to eliminate the possi- 
bility of an interaction independent of spin 
(case 1 discussed in the introduction) for if 
there were no spin dependence the differences 
should be in the opposite direction. The problem 
becomes therefore one of distinguishing between 
the two cases involving spin dependent forces, 
(2a) a real singlet deuteron level, and (2b) a 
virtual singlet deuteron level. Schwinger and 
Tellert have shown that the dependence of the 
para cross section on neutron energy will be very 
different in the two cases. If the singlet level is 
real the para cross section should be slightly 
greater (about 30 percent) for “‘liquid air”’ 
(0.012 ev) neutrons than for ‘‘room tempera- 
ture” (0.037 ev) neutrons, while if the singlet 
level is virtual the cross sections for “‘liquid air” 
neutrons should be ~100 times smaller. More- 
over if the singlet state should be real the ortho- 
hydrogen cross section would be only about 40 
percent greater than the parahydrogen cross 
section for liquid air neutrons, whereas for a 
virtual state the difference may be very large. 
Fig. 3 shows (1) that there is a marked decrease 
in the para cross section with decreasing neutron 
energy and (2) that the ratio of the ortho to para 
cross sections for the ~120°K neutron distribu- 
tion was 4 : 1. Our results therefore prove con- 
clusively that the singlet state of the deuteron 
is virtual. This confirms the work of Halpern, 
Estermann, Simpson and Stern.*® Taking a more 


% J. Halpern, I. Estermann, O. C. Simpson and O. Stern, 
Phys. Rev. 52, 142 (1937). 
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general view of the matter it may be said that 
the virtual level hypothesis is the only one which 
leads to large differences in the para and ortho 
cross sections, and that none of the other 
possibilities can account for differences as large 
as those observed. 

As Schwinger® pointed out, these results also 
require that the neutron spin be 2 rather than 3, 

The presence of the tail of high energy ney- 
trons unfortunately makes it impossible to gain 
any information concerning the range of the 
neutron-proton forces from our experiments, 
The fact that the observed cross section of para- 
hydrogen for ~120°K neutrons is not as small 
as the theory predicts is to be attributed to these 
high energy neutrons, whose inelastic scattering 
completely masks the small elastic para scatter- 
ing (0—0) which depends sensitively on the 
range of the neutron-proton forces. 

It should be borne in mind that Schwinger and 
Teller’s calculations give cross sections not as a 
function of temperature, but as a function of the 
energy of a monochromatic neutron beam. Also 
their calculations are strictly applicable only to 
hydrogen in the gaseous state, whereas our 
measurements apply to the liquid. These facts, 
together with the excess of high energy neutrons 
in the distributions with which we worked, are 
sufficient to account for the differences between 
our results and the predictions of Schwinger and 
Teller. Specifically, they are sufficient to account 
for the fact that the ortho cross section increases 
somewhat more rapidly with decreasing neutron 
‘temperature’? than was predicted, while the 
para cross section does not decrease as rapidly. 

The possibility referred to earlier of using para- 
hydrogen as a filter to obtain slow neutrons does 
not at first sight appear extremely promising. 
The ratio 29.0/17.6 found for the cross sections 
of parahydrogen for ~300°K and ~120°K 
neutrons is far from the 100 to 1 ratio”® predicted 
by Schwinger and Teller for 0.037 and 0.012 ev 
neutrons. However, the mean energy (3KT7/2) 
corresponding to 120°K is not 0.012 but 0.015 ev, 
and it must be remembered that 300°K and 
120°K are to be interpreted not as actual tem- 
peratures but as lower limits. It seems reasonable 


26 As the authors point out this figure cannot be con- 
sidered as a precise estimate, since it depends on the 
assumed range of the force of interaction. 
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therefore that a considerable number of the 
~120°K neutrons should have energies greater 
than 0.023 ev, which is the point at which 0-1 
scattering begins to take place and the para 
cross section begins to increase rapidly. It seems 
probable therefore that the ratio of para cross 
sections for 0.037 and 0.012 ev neutrons is 
actually much larger than 29.0/17.6. A layer of 
parahydrogen of suitable thickness should there- 
fore transmit neutrons in the low energy range 
of a Maxwell distribution, and cut off rather 
sharply at an energy of 0.023 ev. Such a low 
energy beam shou’. be quite useful in experi- 
mental work. 

In the actual transmission measurements and 
calculations of cross sections, the chief errors 
involved are in cell thickness, composition of 
the hydrogen, and statistical fluctuations of the 
counts. Of these the largest is probably in the 
determination of‘ cell thicknesses. Although this 
will affect the absolute values of the cross 
sections it will have very little influence on the 
relative values for different ortho-para composi- 
tions. Boiling may reduce the effective thickness 
of material in a cell. An estimate of the rate of 
heat transfer to the cells by radiation indicated 
that the amount of gas produced by evaporation 
in one second would not fill more than 2 percent 
of the cell volumes. We have observed in other 
experiments that the bubbles in liquid hydrogen 
are small and rise rapidly. Errors due to the 


presence of bubbles should be of the order of 1 
percent. Uncertainty in the measurement of the 
composition is 1 percent. The number of counts 
was sufficient to make the probable error due to 
statistical fluctuations less than 2 percent. The 
combined error due to all the above factors 
certainly cannot be greater than 10 percent. 

The results of the measurements on liquid Ne, 
Os, and H.O are in substantial agreement with 
previous values.*: Log transmission has been 
plotted against scattering thickness in Fig. 4. 
In obtaining proton cross sections from the water 
measurements the value 4.0 10~* cm? for the 
cross section of oxygen was taken from Table II. 
The same value was used for both the ~300°K 
and ~120°K neutrons since no information on 
the temperature dependence was available. The 
ratio of the proton cross sections in water thus 
obtained is o120/¢300= 1.29. Errors in these meas- 
urements are substantially the same as in the 
hydrogen measurements. 

The authors wish to express their appreciation 
of the cooperation of Director Lyman J. Briggs 
of the National Bureau of Standards and Dean 
George B. Pegram of Columbia University which 
made these experiments possible. They are also 
indebted to Professor Edward Teller of the 
George Washington University and to Mr. 
Julian Schwinger of Columbia University for 
valuable discussion and advice. 
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A mass-spectrographic study of the isotopic constitution of Sr, Ba, Bi, and Tl has been made. 
Although the sensitivity and resolving power of the apparatus were sufficient to detect ex- 
tremely rare isotopes, no new ones were found. Relative abundances of the isotopes of Sr, Ba, 
and Tl were determined. It was possible to show that if Hg?® exists it must be with an abun- 
dance less than one-tenth of that at one time assigned to it by Aston. 


LTHOUGH a knowledge of the existence or 
nonexistence of rare stable isotopes is of 
great value in the study of nuclear physics, 


* National Research Fellow. 


relatively few elements have been investigated 
with apparatus capable of detecting isotopes 
present in only small amounts. The present 
paper is a report on the continuation of a sys- 


tematic program, the aim of which is to establish 
the existence or nonexistence of possible rare 
isotopes, as well as to determine the relative 
abundances of the accepted stable ones. Because 
of the extremely high sensitivity and resolving 
power of the apparatus it is possible to extend 
the search for rare isotopes beyond previously 
attainable limits. 
APPARATUS 


A description of the mass spectrometer em- 
ployed has already been given.! It may be 
described briefly as an apparatus in which 
positive ions are formed by the collision of a 
controlled beam of low energy electrons with the 
gas or vapor of the elements to be studied. The 
positive ions are drawn from the ionizing region 
by suitable electric fields and the m/e analysis is 
accomplished by means of a 180° magnetic 
analyzer. Inasmuch as the analyzed ion current 
is measured by means of an electrometer tube 
it is relatively easy to obtain accurate determina- 
tions of relative abundances. Because of the fact 
that the metal parts of the mass spectrometer 
tube are constructed entirely of Nichrome, 
copper, and tungsten, and are surrounded by a 
Pyrex envelope without grease or wax joints, the 
tube may be thoroughly baked, and hence, little 
difficulty is encountered with impurities. In the 
study of relatively nonvolatile substances such 
as Sr, Ba, Bi, and TI, the element is placed in a 
small Nichrome furnace inside of the spec- 
trometer tube. The furnace may be heated to a 
temperature as high as 1200°C and the stream of 
vapor emitted passes into the electron beam 
where ions then are formed. For the investigation 
of volatile substances the furnace is, of course, 
unnecessary, and the gas or vapor can be ad- 
mitted to the tube through a capillary leak or by 
some other means. 

In the present work the slits shown in Fig. 1 of 
reference 1 had the following widths: S; =0.007”, 
S2=0.006", and S3, the exit slit of the 180° 
analyzer (not shown in the figure) =0.006” for 
Hg, Tl, and Bi; 0.018” for Sr and Ba. 


RESULTS 


Strontium 
A clean sample of strontium metal obtained 
from Eimer and Amend, labeled 99.9 percent 


1 Nier, Phys. Rev. 53, 282 (1938). 
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Fic. 1. Mass spectrum showing isotopes of strontium. 
The experimental points are obtained by plotting positive 
ion current passing through the exit slit of the analyzer as 
a function of the energy of the ions while the magnetic 
field is held constant. As there is a definite relation be- 
tween the energy and the mass of the ions collected, 
namely, energy X mass=a constant, the voltage scale can 
be readily changed to a mass scale, ‘as was done in the 
figure. The width of the peaks is, of course, due to the 
finite width of the slits used and the imperfections in 
focusing. 


pure, was placed in the small furnace, shown in 
Fig. 1 of reference 1. After thorough baking of 
the tube to remove impurities, work was begun. 
It was found that a furnace temperature of the 
order of 700°C gave sufficiently large positive ion 
currents without too great a rate of consumption 
of strontium. A typical mass spectrum is shown 
in Fig. 1. The isotopes observed were the same 
as those reported by other investigators.?~* In 
Table I are given the relative abundances found 
in the present work as well as those given by 
Sampson and Bleakney.‘ The results are seen to 
be in good agreement. The writer believes his 
relative abundances as measured for Sr‘? and 
Sr** relative to Sr** to be correct within one 
percent, while that given for Sr“ may be in 
error by two percent. A possible source of error 
not taken into account in these measurements is 
the effect of the free evaporation of Sr atoms 
from the small furnace. In the case of ideal 
distillation this would mean that the Sr*/Sr* 
ratio as measured would be (88/84)!=1.023 
times greater than the true ratio. However, as 
the Sr atoms are evaporated from a solid surface 
2 Blewett and Sampson, Phys. Rev. 49, 778 (1936). 
3 Dempster, Phys. Rev. 50, 186 (1936). 


4Sampson and Bleakney, Phys. Rev. 50, 456 (1936). 
5 Mattauch, Akad. Wiss. Wien Ber. 145, 461 (1936). 
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and not from a liquid (where mixing of particles 
on the surface with those inside can take place) 
it does not seem likely that one would have 
appreciable discrimination. 

If one assumes the packing fraction of stron- 
tium to be that found by Mattauch® for Sr**, 
namely, —9.0+0.5, and if one uses a conversion 
factor 1.000275 in going from the physical to 
the chemical scale, an atomic weight 85.61 is 
computed for strontium from the writer's as well 
as from Sampson and Bleakney’s relative abun- 
dances. On the other hand, if one assumes the 
packing fraction of Sr to be the same as that 
found by Aston’ for Kr**, namely —7.40, the 
atomic weight 85.62 is computed. From the 
packing fraction — 6.3, deduced from Dempster's* 

cking fraction curve, one obtains the atomic 
weight 85.63. The chemically determined value 
is accepted as 85.63. 

The real reason for studying strontium was for 
the purpose of searching for possible rare iso- 
topes. By analogy to most of the heavier even 
elements one might expect to find strontium to 
consist of a sequence of isotopes having suc- 
cessive mass numbers flanked by a light isotope 
two units removed from the lightest isotope in 
the sequence and a heavier isotope two units 
above the heaviest isotope in the sequence. 
Thus, on this basis, one should expect to find a 
stable isotope, Sr®*°. However, a careful investiga- 
tion of this position as well as the other vacant 
positions around the known stable isotopes failed 
to show even the faintest trace of any new 
isotope. It was possible to set the following 
upper limits for abundances in these positions 
relative to Sr®*: Sr®?, Sr®, and Sr%, 1/300,000 ; 
Sr**, 1/100,000; Sr*, 1/50,000; 1/100,000 ; 
Sr®, Sr*!, and Sr*®, 1/200,000. 


Barium 
Barium metal obtained from Eimer and 
Amend, labeled 99 percent pure, was employed 


TABLE I. Relative abundances of isotopes of strontium. 


MEAN Mass 


| 


87 86 84 NUMBER 
Abundance 100 | 8.50 | 11.94 | 0.68 87.710 
S. & B. values 100 | 9.10 | 11.64] 0.61 87.713 


® Mattauch, Naturwiss. 25, 170 (1937). 
7 Aston, Nature 140, 149 (1937). 
* Dempster, Phys. Rev. 53, 869 (1938). 


as a source of barium. The procedure was the 


same as that used for strontium. Fig. 2 shows a 


typical mass spectrum obtained. An idea of the 
sensitivity of the apparatus is gained when one 
considers that the points for the rare isotopes, 
Ba and Ba™®, are plotted to a scale 100 times 
that for the other isotopes. The isotopes found 
are the same as those reported by Dempster® 
and Sampson and Bleakney.‘ Table II gives the 
relative abundances found in the present work, 
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Fic. 2. Mass spectrum showing the isotopes of barium. 
Current scale in region of Ba! and Ba!™ is multiplied by a 
factor of one hundred. 


together with the abundances as given by 
Sampson and Bleakney. The writer believes his 
abundances relative to Ba"** all to be correct 
within two percent except the ones for Ba and 
Ba, which may be in error by four percent. 
As in strontium the evaporation of the atoms 
takes places from the surface of a solid, and 
hence, one should expect little, if any, separation 
to take place in the evaporation process. 

If one assumes the packing fraction of barium 
to be the same as that given recently by Aston’ 
for Xe, —4.4, and if one uses a conversion 
factor 1.000275 in going from the physical to 
the chemical scale, an atomic weight 137.33 is 
obtained from the writer’s data. With a packing 
fraction —4.0, as deduced from Dempster’s* 
curve, the value 137.33 is also found. The 
chemically determined value is given as 137.36. 

As in the case of strontium, the main purpose 
in investigating barium was to search for possible 


® Dempster, Phys. Rev. 49, 947 (1936). 


278 ALFRED O. NIER 


TABLE II. Relative abundances of isotopes of barium. 


MEAN 

138! 137 | 136 | 135) 134| 132 | 130} Mass No. 
Abundance 100} 15.8) 10.9} 9.2 | 3.37,0.136)0.141 137.42 
S. & B. values 100} 14.9} 11.7) 7.8 | 2.35|0.022/0.22 137.47 


rare isotopes. However, none was found, and it 
was possible to set the following upper limits for 
abundances of hypothetical isotopes relative 
to Ba™§: 1/200,000; and 
1/100,000; Ba’, 1/35,000; Ba™!, Bal, 
and Ba!’, 1,100,000. 


Bismuth 


Aston’ found bismuth to consist of a single 
isotope of mass 209. This fact was verified by 
Bainbridge and Jordan" who showed that no 
other isotope existed with an abundance greater 
than two percent. From the present study one 
can set the following upper limits of abundance 
for the existence of hypothetical isotopes relative 
to Bit, and Bi*", 1/100,000; Bi?!° 
and 1/50,000; Bi?*, and Bi?®, 
1/100,000. Bismuth is thus seen to be single to 
a very high degree, indeed. 

Thallium 

Thallium, an element having an odd atomic 
number, was found by Schiiler and Keyston™ to 
consist of two isotopes having the masses 203 
and 205. The abundance ratio 203/205 was given 
as 0.435. The existence of these two isotopes was 
verified by Aston.” As it appears to be a general 
rule that odd elements do not have more than 
two stable isotopes it did not seem likely that any 
additional ones would be found. However, very 
few elements have been investigated with appa- 
ratus capable of detecting truly rare isotopes, 
and hence a study of thallium is of interest if for 
no other reason than to better establish the 


validity of the rule. 


1 Aston, Phil. Mag. 49, 1197 (1925). ; 

" Bainbridge and Jordan, Phys. Rev. 50, 282 (1936). 

® Schiiler and Keyston, Zeits. f. Physik 70, 1-2, 1 (1931). 
13 Aston, Proc. Roy. Soc. Al34, 571 (1932). 


Chemically pure thallium metal, obtained 
from Eimer and Amend, was used as a source of 
thallium. The 203/205 ratio was found to be 
0.410+2 percent, a number in excellent agree. 
ment with the ratio 0.417 found by Aston, 
As the evaporation of the metal takes place from 
the liquid phase the abundance ratio as measured 
may be high by a factor, (205/203)!= 1.005, due 
to the effect of ideal distillation. Account was 
taken of this factor in arriving at the number 
0.410. 

The abundance ratio found corresponds to a 
mean mass number 204.42+0.01. With a packing 
fraction of +2.0 derived from Dempster's" data 
and a conversion factor 1.000275 in going from 
the physical to the chemical scale, an atomic 
weight 204.40 is computed from the writer's 
data. This is to be compared with a chemically 
determined value 204.39. 

No new isotopes were found, and it was 
possible to set the following upper limits of 
abundance for the existence of possible rare ones 
relative to TI, TP°7, and 
1/50,000; TI? and TI, 1/30,000; TH, 
and TI!*, 1/50,000. 


Mercury 

Hg?" was at one time reported by Aston" to 
exist as a stable isotope with an abundance 
approximately 1/5000 that of Hg?. Recently 
the writer!® undertook a search for this isotope, 
but because of limited resolving power was able 
to conclude only that if it existed at all it must 
be with an abundance less than that assigned to 
it by Aston. In the present work on Bi and TI 
the resolving power of the apparatus was some- 
what higher than was the case in the earlier work, 
mainly because of the narrower slits used. With 
a pressure of around 10~° mm of Hg in the tube, 
the sensitivity was sufficient to show definifely 
that if Hg? existed it must be with an abundance 
less than 1/50,000 that of Hg?. 

4 Dempster, Phys. Rev. 53, 64 (1938). 


% Aston, Nature 130, 847 (1932). 
® Nier, Phys. Rev. 52, 933 (1937). 
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Domain Theory of Ferromagnetics Under Stress* 


Part III. The Reversible Susceptibility 
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It is shown that the statistical theory of Parts I and II can be extended so as to give formulas 
for the magnetization curve and other properties of an ideal reversible specimen, over the 
range of field values in which magnetization proceeds by a displacement of the boundaries 
between domains. The application to actual specimens is made by assuming that the reversible 
properties of an actual specimen are identical with the properties, at the same magnetization, 
of an ideal specimen having the same initial susceptibility. This leads to a formula for the 
reversible susceptibility of nickel and iron that was originally suggested by Gans and verified by 
him experimentally, but for which it is believed no satisfactory derivation has hitherto been 
offered. In the case of one nickel specimen and of cobalt, where the data do not follow the Gans 
curve, formulas agreeing better with experiment are obtained by taking account of the aniso- 
tropy of the domains. Reasons are given for believing that the replacement of variable by fixed 
domain volumes, an arbitrary step taken in order to simplify the problem, is probably not 


essential to the result. 


HE behavior of ferromagnetic crystals at 
fields of order of magnitude 10? oersteds has 
been satisfactorily explained! by supposing that 
an increase of field causes a “‘rotafion’”’ of the 
magnetization vector toward the field direction. 
This vector has a constant magnitude, the 
saturation value J,; and except in special cases, 
it has also a uniform direction throughout each 
crystal—the direction for which the free energy, 
consisting of a term due to the field and a term 
due to the crystalline anisotropy forces, is a 
minimum. The special cases just referred to are 
those in which there are several such directions, 
as is the case when the field direction is [110] or 
[111] in an iron crystal. In these cases the 
crystal is saturated not as a whole, but only 
throughout ‘‘domains,’’ equal volumes being 
magnetized in each of the equivalent directions. 
The properties of large grained polycrystalline 
specimens in this same range of field values can 
be calculated by averaging those of crystals.* 

At fields of order of magnitude 1 oersted, the 
process is different and a successful quantitative 
theory is still lacking. Certain fundamental 
ideas, however, have proved useful and are 

* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1 An outline of the theory and references to the original 
papers are given in E. C. Stoner, Magnetism and Matter 
(1934), pp. 392 ff., and F. Bitter, Introduction to Ferro- 


magnetism (1937), pp. 194 ff. 
*R. Gans, Ann. d. Physik 15, 28 (1932). 


probably essentially correct. It is supposed that 
even in the demagnetized state the specimen is 
saturated throughout domains; equal volumes 
are magnetized in each of the possible directions, 
which are now the directions of minimum free 
energy determined by the anisotropy forces 
alone—the ‘‘directions of easy magnetization.” 
When a field is applied and gradually increased, 
the electron spins responsible for the mag- 
netization undergo transitions in the boundary 
region between domains; the effect is a motion of 
the boundary, so that the more favorably 
oriented domains grow at the expense of their 
less favorably oriented neighbors. This displace- 
ment of the boundaries is in part reversible, 
proceeding an infinitesimal distance dx when the 
field is changed by an infinitesimal amount d//; 
and in part irreversible, for when a small change 
of field enables the boundary to pass a local 
energy barrier, the motion may continue a 
considerable distance, with dissipation of energy 
through the setting up of local eddy currents and 
in other ways. 

The most direct evidence of the correctness of 
this picture is obtained from measurements of 
the Barkhausen effect*® and from the study of 
magnetic powder and colloid patterns.‘ The 
regions whose reversal of magnetization is 


* For references see Stoner, reference 1, p. 436. 
‘F. Bitter, reference 1, p bp: 59 ff.; W. ral ‘Elmore, Phys. 
Rev. 51, 982 (1937); 53, 75 ( 


279 


1938). 
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observed in the Barkhausen effect, although not 
necessarily identical with the domains, at least 
set a lower limit to their order of magnitude. It is 
consistent with the order of magnitude of the 
dimensions observed in surface powder patterns, 
but the latter method of study shows that the 
actual shapes and arrangements of the differently 
magnetized regions may be very complicated. 
Various types of force are operative in this 


process. The exchange forces are important only . 


in the boundaries between the domains, where 
their partial annihilation has the effect of a 
positive surface energy proportional to the total 
area of the interdomain surfaces, and therefore 
plays a role in determining the size of the 
domains. This size in fact seems to be determined 
primarily by a balance between these forces and 
magnetic forces dependent upon the size and 
shape of the specimen, so that the domain size 
varies approximately as the square root of the 
shortest dimension of the specimen.’ The 
anisotropy forces are important chiefly at higher 
fields, but according to one theory® they are also 
instrumental in determining the thickness of the 
interdomain space and hence indirectly the size 
of the domains. The theory of the behavior of an 
ideal crystal in an applied field has been attacked 
by several authors,5~? but even in the simplest 
cases the problem is difficult because nonlinear 
partial differential equations occur at the 
outset. 

In an actual crystal, though the order of 
magnitude of the domain size is probably de- 
termined by the factors already mentioned, the 
exact behavior in a field is dependent very 
largely on forces that vary irregularly from one 
point to another. Chief among these are internal 
stresses due to impurities or to imperfect an- 
nealing. When these have been reduced to a 
minimum, there remain unavoidable stresses 
resulting from the unequal magnetostriction of 
differently oriented domains. Becker* has shown 
how the internal stresses limit the initial sus- 
ceptibility in a somewhat idealized model, and 
Kersten® has shown that the magnitude of the 


i Frenkel and J. Dorfman, Nature 126, 274 (1930). 
°L. Landau and E. Lifshitz, Physik. Zeits. Sowjetunion 
8, 153 (1935). 

7F, Bitter, reference 1, pp. 185-192. 

§R. Becker, Physik. Zeits. 33, 905 (1932). 

»M. Kersten, Zeits. f. Physik 12, 665 (1931). 
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highest attainable initial susceptibility in iron. 
nickel alloys of various compositions is ade. 
quately accounted for by setting the internal 
stress equal to the lower ‘limit determined by 
magnetostriction. In addition to this effect of 
mild stresses in controlling the reversible bound- 
ary displacement, there is the irreversible effect 
of local high stresses, which set up energy barriers 
that lead to Barkhausen reversals and hence to 
hysteresis.” 

The importance of irregular forces suggests the 
possibility of handling the problem by a sta- 
tistical method, in which the exact nature of the 
process need not be considered in detail. This 
method of attack was originated by Heisenberg. 
The quantities that have so far been calculated 
by it are the magnetostriction as a function of 
magnetization ; the normal component of mag- 
netization as a function of the component parallel 
to the field ; and the elastic constants as functions 
of magnetization, given the magnetization curve. 
In parts I and II under the present title," 
previous work in this field was summarized, more 
general formulas were developed, and some 
errors in earlier calculations were pointed out and 
corrected. No new principles were introduced, 
except the device of using thermodynamic rela- 
tions to eliminate certain unknown quantities; 
by this means it was possible to avoid an 
incorrect approximation made in some of the 
earlier work. 

It is possible to extend the theory so as 
actually to calculate the magnetization curve of 
an ideal reversible specimen, and from this all the 
reversible properties of an actual specimen, given 
a single reversible property—the initial sus- 
ceptibility. This extension is the subject of the 
present paper; but in order to make the presen- 
tation convincing and to insure against misin- 
terpretation of the results, it will be necessary to 
review the entire theory briefly, from a somewhat 
different point of view. 

A few general remarks will make the aim 
clearer. The Becker-Kersten rotation theory, 
applicable only when the internal stresses o; are 
large and fairly uniform over a domain, gives for 


Kondorsky, Physik. Zeits. Sowjetunion 11, 597 
(1937). 

1 W. F. Brown, Phys. Rev. 52, 325 (1937); 53, 482 
(1938): hereafter cited as I and II, respectively. 
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the initial susceptibility xo=(2/9)J.2, Aw is 
the saturation magnetostriction. Becker's sim- 
plified theory of boundary displacement against 
small, inhomogeneous stresses of amplitude o; 
gives xo=(1/7)J.2/d.0;. Although the mecha- 
nisms assumed are entirely unlike, the results are 
practically identical. Again, it is possible to 
develop a theory similar to that of Becker and 
Kersten, except that internal stresses are replaced 
by internal fields ; the result of either theory may 
be written in the form xo =J,*?/e, where ¢ is an 
energy density associated with the “random’’ 
forces; and in each case the area to the left of the 
magnetization curve, /HdJ, is approximately 
equal to the ¢ of the xo formula. In other words, 
the most important characteristics of the theo- 
retical magnetization curve—its initial slope, 
and the area between it and the J axis—are 
determined by energy magnitudes rather than by 
the exact nature of the internal forces or the 
exact mechanism through which they make their 
effect felt. It should therefore be possible to 
develop at least an approximate theory in which 
no specific assumptions are made in regard to the 
nature of the forces or the mechanism. The 
advantage of such a theory is that it avoids the 
danger of misinterpreting an agreement between 


‘theory and experiment as evidence for the 


particular forces or mechanism assumed in the 
theory. 
GENERAL THEORY 


The theory will take no account of irreversible 
phenomena, and therefore the properties pre- 
dicted by it will be those of an ideal reversible 
specimen. It is necessary to relate these in some 
way to the properties of an actual specimen. Now 
there are certain properties of a ferromagnetic 
specimen which, to a first approximation, are 
determined by the magnetization alone and are 
independent of the previous history and hence of 
the field strength associated with that mag- 
netization in the particular state under observa- 
tion. The quantities of which this is true are those 
measured in a small, approximately reversible 
change from any initial state: for instance, the 
elastic constants” and the reversible suscepti- 
bility’ (ratio AJ/AH when is small and 

® E. Giebe and E. Blechschmidt, Ann. d. Physik 11, 905 


(1931); O. von Auwers, Ann. d. Physik 17, 83 (1933); 
W. T. Cooke, Phys. Rev. 50, 1158 (1936). 
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negative). This suggests that such quantities 
should have the same value for the actual 
specimen and for a reversible specimen at the 
same magnetization. We assume, therefore, that 
so far as a specimen’s properties are functions of 
the magnetization alone, they are the same 
functions of magnetization as are the properties 
of an ideal reversible specimen having the same 
initial susceptibility. 

The correct model for a single crystal is now a 
group of many domains with fixed orientations of 
their magnetization vectors, but with variable 
volumes; the positions of the boundaries are 
determined by equilibrium between internal and 
external forces, that is by the minimum of the 
total free energy due to external fields and 
stresses and to internal ‘‘random”’ forces. On the 
supposition that the mechanism by which these 
internal forces act is unimportant, we proceed 
to substitute for this model a more manageable 
one, consisting of a group of N domains of fixed 
and equal volumes, but with variable orientations 
of their magnetization vectors : the magnetization 
of each domain is to be along one of the “easy 
magnetization” directions, but may change from 
one of these directions to another when the field 
or stress is changed. The exact microscopic state 
of the model under given external conditions is to 
be that for which the total free energy is a 
minimum. A macroscopic state of the model, of 
course, will be described by merely giving the 
fractional number ,=N,/N of domains mag- 
netized in each of the possible directions o, 
without regard to which domains these are. 

The actual microscopic state, for given field 
components #7; and stress components X ;;, is 
that corresponding to the minimum value of" 


V=V'+V", (1) 


where V’ is the ordered free energy per cm’ due to 
the applied field and stresses, and V” is the 
energy due to the random internal forces. V’ is 


given by 
— LX eis’ = DX jne jn’ 
— LAS ie— LX istic’ 
DX jul ine’) = Lene Ve. (2) 
1%3R. Gans, Physik. Zeits. 12, 1053 (1911); Ann. d. 


Physik 61, 379 (1920). 
4], Eq. (7). The constant term Vo has been dropped. 
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In a large number of samples of the same 
material, microstates of equal V’ will occur 
equally often if the internal forces are random," 
but microstates of lower V’ will occur oftener. 
On the other hand there are many more micro- 
states of zero ordered energy than of highly 
negative. Somewhere between 0 and — ~, there 
is a “‘normal”’ value of V’ that will occur most of 
the time; and among samples with about this 
value of V’, certain values of the partition 
numbers N, will be normal. Just what value of 
V’ is normal for given 77;,’s and X ;;’s depends on 
the average magnitude of the random forces: the 
larger these are, the larger will be the external 
field or stress required to produce a specified 
magnetization or strain. 

Assuming that V’ has some value, we may 
proceed to determine the normal values of the 
N,’s, and hence the macroscopic properties of the 
specimen, by methods similar to those used in 
statistical mechanics. The result is 


ng=exp (—LoV.)/dY.exp (—LoV.), (3) 


where Ly is determined by 
XV. exp exp (—LoV.)=V"’. (4) 


This is in essence the method that was introduced 
by Heisenberg and used by Bozorth and by Gans 
and v. Harlem. These writers, however, did not 
actually assume a given value of V’, but rather a 
given value of the component of magnetization 
in the field direction, J,. In the problems treated 
by them the two quantities differ only by a factor 
H, and moreover it is really J, and not V’ that is 
known and must be used to eliminate the 
Lagrangian multiplier Lo; the two procedures 
therefore lead to the same result. In the work of 
Akulov and Kondorsky stresses as well as fields 
are supposed to be present, and a given value of 
V’ must be assumed. There is now a difficulty, 
however; for V’ is not really known, and so Lo 
remains an undetermined quantity, whose de- 
pendence on the field and stress are unknown. 
Akulov and Kondorsky endeavored to overcome 
this difficulty by supposing the magnetization 
curve at zero stress to be given, and seeking only 


% More accurately, this is our definition of ‘‘random” 
as applied to these forces. 


a formula for the N,’s good to the first order in 
the stresses. But in doing this they failed to take 
account of the effect of stresses upon the mag- 
netization curve, and thus obtained an incorrect 
result. It was shown in I that this error could be 
corrected, and the result they sought obtained, 
by using the thermodynamic relations dJ;/aX ,, 
=e ,/0H;. The equations thus obtained are 
rather complicated, it is still necessary to know 
the magnetization curve, and the physical 
significance of the troublesome Lagrangian 
multiplier never becomes apparent. The theory 
at this stage is clearly incomplete. 

A similar problem is encountered in statistical 
mechanics, in calculating the normal partition 
numbers for a given energy of an isolated system, 
The energy is really unknown, and s0 jts 
Lagrangian multiplier ZL remains indeterminate 
until it has been given some thermodynamic 
interpretation. This may be done in a number of 
ways, for instance by considering ‘wo systems 
capable of exchanging energy. The same equa- 
tions as before are obtained for each system, but 
with the same L for both; L therefore represents 
a quantity which is the same for two systems 
when they are in thermodynamic equilibrium— 
namely, the temperature. (A more thorough 
investigation of course shows that L is pro- 
portional to the reciprocal of temperature on the 
usual absolute scale.) 

This and other methods, suitably modified, 
may be applied here. The phrase “suitably 
modified”’ is important, for it must be emphasized 
that this is not a problem in statistical mechanics. 
Our model is not a group of domains continually 
jumping from one state to another as a result of 
exchanges of energy. It is a group of domains ina 
definite, permanent state of lowest total free 
energy, under the joint action of ordered and 
random forces—all constant. The direct effect of 
temperature agitation on the domains as a whole 
is completely negligible because of their size; its 
indirect effect in bringing about the equilibrium 
of the boundaries is fully taken into account in 
the temperature-dependent ‘“‘constants”’ (such as 
J, and \..) appearing in the free energy function. 
Another way of saying this is that the domain 
order-disorder term in the free energy is negligible 
compared with the other terms. This term is 
—kT log W, where W is the number of ways of 
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arranging the domains in the given macroscopic 
state. For unit volume it is of order of magnitude 
ETN, where N is the number of domains per cm‘. 
For a wire 1 mm in diameter, V=10° according to 
Barkhausen effect measurements, so that this 
term in the free energy density is of order of 
magnitude 10~* erg/cm*. The energy associated 
with the internal forces opposing magnetization 
is measured by the area {F/dJ to the left of the 
steep part of the magnetization curve, or roughly 
by J.2/x0; even in well-annealed specimens, its 
order of magnitude is between 10? and 10 
ergs/cm*. Thus thermal agitation of the domains 
as a whole is negligible by a factor of about 10-°. 

Consider, then, two adjacent volumes 7, v2 of 
the same material, subject to different values of 
the fields or stresses. If we now seek the normal 
N,'s for a given total ordered free energy 
= £;'+ £2’, we obtain equations 
of the form (3) for each part, but the same 
Lagrangian multiplier Lo occurs in both sets of 
equations. Lo is determined by the total Z’, but 
also it is equally well determined by the £’ for 
either part. If action at a distance by the random 
forces is negligible, a variation of the fields or 
stresses in one part does not affect the state of 
affairs in the other, and so does not affect Lo, 
which must therefore be a constant of the 
material (at the given temperature). Its value 
depends on the random forces. For cubic crystals 
or isotropic material, it follows from Eq. (41) of I 
and the corresponding equation for [100] 
domains, and from Eqs. (15), (23), and (36) of II, 
that 


xo= 3LoJ.*; (5) 


hence 1/L is an energy density of the same order 
of magnitude as /J/dJ (if the number of direc- 
tions of easy magnetization is finite), and 
measures the energy associated with the internal 
forces. If we assume that these are internal 
stresses of mean value o;, and set 1/Lo=A..0;, we 
get xo=4J.2/d.0;, a formula differing’ from 
Kersten’s only in having a numerical factor 3} 
instead of 2/9. 

Although the discussion has been limited to 
single crystals for simplicity, the extension to 
polycrystalline material is easily made by the 
method of I. Thus we get for the magnetization 
curve of the ideal reversible specimen of iron or 


nickel Eq. (12), (20), or (34) of II, according to 
the type of domain (isotropic, [100], or [111 ]); 
in the definition of » preceding each of these 
equations, Lo is now to be interpreted as a 
constant, related to the initial susceptibility by 
Eq. (5) above. For materials with other types of 
domain, Eq. (3) or (4) of IT must be used. 


THE REVERSIBLE SUSCEPTIBILITY 


The magnetization curve thus obtained in 
terms of xo or Lo, for an ideal reversible specimen, 
is different from any J—// curve of the actual 
specimen. But the differential susceptibility 
dJ/dH calculated from it corresponds to a 
measurable susceptibility for actual specimens, 
namely the ‘‘reversible susceptibility’’ mentioned 
before. The value of the reversible susceptibility 
x, at any magnetization should therefore be 
equal to the differential susceptibility of the 
ideal specimen at the same magnetization, if Lo is 
chosen so as to give the observed initial suscepti- 
bility. (The corresponding value of /7, however, 
bears no simple relation to the // of the actual 
specimen.) 

Thus for material with isotropic domains, x, is 
determined as a function of J by the parametric 
equations 


J/J,=L(n), xr/x0=3L'(n), (6) 
L() =coth »—1/n. (7) 


These equations were first proposed by Gans in 
1911, and were supported by experimental data 
agreeing very well with the calculated curve." 
Gans gave no derivation of his formula; he was 
probably led to it by an attempt to apply the 
Langevin formula for the susceptibility of a 
paramagnetic gas. A derivation along that line 
was very apologetically offered by Debye'® in 
1925. According to this interpretation, 7 should 
be equal to wH/kT, where uw is the magnetic 
moment of an elementary magnet, in Debye’s 
picture a crystal. Today we must interpret the 
elementary magnet as a domain. Its magnetic 
moment is J,vo, where vop10-*. The initial 
susceptibility thus calculated is 10° times too 
large, as was to be expected from our earlier 
discussion. 


where 


4 P. Debye, Handbuch der Radiologie, Vol. 6 (1925), 
p. 721. 
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Fic. 1. Reversible susceptibility of polycrystalline nickel. 
Experimental points (D. Kirkham): open circles, 21.6°C; 
circles black at left, 99.0°C; black at right, 171.5°C. 
Theoretical curves: 1, J, in any direction; 2, J, along any 
[111] direction; 3, J, along +H; 4, J, along [111] direc- 
tions nearest +H. 


If the domains are crystalline, in iron or nickel, 
the functions are more complicated in form, but 
the curves differ very little from this one; in 
fact the series expansions for J are identical to 
the sixth order in H (see Eqs. (15), (23), and (36) 
of II). This explains Gans’s success in fitting his 
curve to data on various specimens of iron, steel, 
and nickel. 

Figure 1 shows several attempts to fit 
Kirkham’s" results on annealed nickel. These 
measurements were not taken on the same 
ellipsoid as the magnetostriction data analyzed in 
II, but on a thinner ellipsoid, 73’’ long and }” in 
diameter. For this specimen the domains are 
probably crystalline at room temperature; but 
according to the value of xo the anisotropy forces 
are not enormously larger than the random 
forces, so that the neglected rotation effect may 
not be entirely negligible, and the specimen may 
be expected to be intermediate in its properties 
between one with crystalline domains of high 
anisotropy and one with isotropic domains. 

Curve 1 is calculated for isotropic domains, 
curve 2 for crystalline domains with [111] as the 
direction of easy magnetization. There is a 
difference only at high J, and here 2 is unreliable 


17D. Kirkham, Phys. Rev. 52, 1162 (1937). 


because of rotation. The data do not fit at all 
well. 

In the discussion of the magnetostriction of an 
iron ellipsoid (II, p. 488), it was suggested that 
when an elongated ellipsoid is demagnetized by a 
longitudinal field, more domains are left oriented 
along directions nearly parallel to the ellipsoid 
axis than along directions nearly perpendicular 
to it. This applies to the present case. Curve 3 js 
calculated for the extreme case in which all the 
domain magnetizations are either parallel or anti- 
parallel to the field (as might be the case if the 
specimen were long and thin and the domains 
magnetically isotropic). Curve 4 is calculated for 
domain moments oriented along only the [111] 
directions nearest the field.'* Most of the data fall 
between these curves. As the temperature rises, 
the anisotropy constant decreases, the effect of 
rotation (or the tendency toward _ isotropy) 
increases, and in the range where rotation js 
important the points should fall nearer to the 
isotropic curve—as in fact they do. 

In Fig. 2, the points are experimental data of 
Samuel’ on soft cobalt. Curve 1 is the Gans 
formula, which does not fit at all. Cobalt crystals 
have two opposite directions of easy magneti- 
zation ; the anisotropy constant is large, so that 
rotation should be negligible except at high 
fields, and the calculation for crystalline domains 
should be a good approximation. This calculation 
gives 


J/J,=E(n), xr/xo=3E'(n), (8) 
where E(n) =(1/n) f u tanh udu. (9) 
0 


Series expansions for the function E(m) may be 
obtained from the expansions in II, Eqs. (22)- 
(27). The function (8) is plotted as curve 2 in 
Fig. 2. 

The deviation at high J is expected. The 
deviation at low J may be partly due to a poor 
value of xo. The experimental results deviate 
considerably among themselves; the agreement 


18 The formulas are: curve 3, J/J,=tanh 7, x-/xo=sech*9 
=1—(J/J,)?; curve 4, J/J,= (12/2) fo" "dv tanh anda, 
(12/m) fo" dy Jala? sech? anda, a’ =cos 6, 
cot 6’=v2 cos y; the result of the integration with respect 
to a can be expressed in terms of the functions £ and E’ 
defined by Eq. (9). 

19 M. Samuel, Ann. d. Physik 86, 798 (1928): specimen H. 
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with theory is probably as good as can be 
expected. Gans* found a considerably worse 
deviation between the experimental results at 
higher fields and the theoretical curve calculated 
by the rotation theory. 

Cobalt is interesting because, with only two 
directions of easy magnetization, it is easy to 
work out a more specialized theory, with more 
definite assumptions about the internal forces. 
Internal stresses, being impartial with respect to 
opposite directions, should have no effect on the 
reversible magnetization, and only internal fields 
need be considered. Curve 3 is calculated on the 
assumption that the magnetic moment of each 
crystalline domain is along that one of the two 
possible directions favored by the resultant //, 
external plus internal; it is assumed that the 
values of the relevant component of internal 
field are constant and are distributed according 
to a Gaussian error curve. This gives 


0 


0 


The difference between curves 2 and 3 is unim- 
portant at the present stage of development of 
the theory, and the wider applicability of the 
general method is strongly in its favor. 


MODIFICATION OF EARLIER FORMULAS 


The theorem that Lo(= —ZL of I) is a constant 
enables us to replace Eq. (16) of I by the simpler 
relation, which now follows immediately from 


(14), 
aY,/dye=LP,». (12) 


Hence only the first terms remain in the right 
members of Eqs. (34) when the differentiation is 
at constant 7/7. From the second of Eqs. (44), 
x1=0. This explains why good agreement with 
experiment was obtained by neglecting x; in Eq. 
(54). It is now clear, also, that the coefficients 
x1 and xz in Eq. (40) have no relation to the 
actual magnetization curve of a crystal, but are 


It is possible to compute the variation of the elastic 
constants with magnetization without introducing the 
approximations of I, but this work is too far from com- 
pletion to discuss at present. 
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determined by the curve of the equivalent 
reversible specimen, which may be obtained 
from the x, vs. J curve by the formula 


J 
(1/x,)dJ. (13) 


Here //, is the field strength corresponding to the 
magnetization J for an ideal reversible specimen 
whose reversible properties coincide with those of 
the actual specimen. 


CRITICISM OF THE THEORY 


The fact that the theory gives only reversible 
properties is a serious limitation, but less serious 
than might have been feared, provided the 
concept of the ‘equivalent reversible specimen” 
is used as a guide in interpreting the theoretical 
results. The assumption of complete randomness 
of the internal forces is an oversimplification ; 
whether it is too drastic an oversimplification can 
only be decided by comparison of the predictions 
of the theory with experiment, and it would not 
be at all surprising if the theory proved useless 
for pure, well-annealed crystals. The replacement 


Fic. 2. Reversible susceptibility of polycrystalline cobalt. 
Experimental points (M. Samuel): open circles, virgin 
curve; solid circles, ideal curve; circles black at bottom, 
descending magnetization curve; black at top, ascending. 
Theoretical curves: 1, general theory, isotropic domains; 
2, general theory, crystalline domains; 3, special theory of 
internal fields. 
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of variable volumes with fixed magnetization 
directions by fixed volumes with variable 
magnetization directions is the step that is most 
difficult to justify. If this step is defensible, then 
the final results should be capable of derivation 
without it. Some efforts in that direction are 
presented in the appendix. 

Despite these limitations, the theory seems 


WILLIAM FULLER BROWN, JR. 


successful enough to justify its further study; jt 
may help to give a partial understanding of the 
complex processes that occur at low fields, until a 
more detailed theory becomes possible. 

The writer wishes to express his appreciation to 
Professor Francis Bitter, of the Massachusetts 
Institute of Technology, for helpful comments 
and criticisms. 


APPENDIX 


(a) Alternative derivations of the distribution 
formulas. 


These parallel familiar arguments of statistical 
mechanics. 

Method 1.—Assume, as before, that in a large 
number of copies of our simplified model, all 
subject to the same fields and stresses, microstates 
of equal ordered energy occur equally often ; and 
that this remains true even if the fields and 
stresses have different values in two adjacent 
parts v1, 2 into which the model may be divided. 
Then the relative frequency of occurrence of any 
microstate, for either part separately or for the 
model as a whole, is a function only of the corre- 
sponding ordered energy Eo’ =v2V2’, 
or V;'+v2V 2’: let this function be f(y’), 
g(E2’), or h(E’). The probability that a copy 
selected at random is in a specified microstate, 
with ordered energies E,’ and E,’, is 


h(E\' + E2’) 


For two microstates with the same total energy 
E,'+£,’, the values of f(£,')g(E2’) are therefore 
equal ; or, f(x)g(y) =f(z)g(x+y—z), for arbitrary 
x, y, and z. Differentiation with respect to x and y 
gives f’(x)g(y) =f(2)g’(x+y—2) =f(x)g’(y), or 
f'(x)/f(x) = g'(y)/g(y) = constant = A. Hence 
f( Ey’) = Be-4""", g(E2’) = Ce~4**’. By varying the 
fields or stresses in v2 without varying them in 7, 
and assuming action at a distance negligible, we 
infer as before that A is a constant of the 
material. Thus for unit volume, the probability 
of a microstate of ordered energy V’=}>on,V, is 


Be-4’’; B is determined by }>Be~4”’ =1, where 


> denotes a summation over all microstates. The 


mean n,’s for the whole collection of copies are 
now given by (”.)~"=B>n,e-4"’. When the sum- 


mation is carried out,*! Eq. (3) is obtained, with 
Li=A/N. 

Method 2.—For given fields and _ stresses, 
the actual microstate of any copy of a volume 
v, containing N domains, is that for which 
E'+E”=0(V'+V”) is a minimum. If we knew 
E” for each microstate of each copy, we could 
solve the problem by first finding the normal 
minimum E”’, say E,’’, attainable by varying the 
microstate, for given N,’s; and then minimizing 
E’+E,” by variation of the N,’s. The £," 
obtained in the first step is a function of the 
N,’s, but if the internal forces are random it 
depends on the N,’s only because the normal 
minimum V” to be found in W microstates jis 
lower, the greater the number W of microstates 
there are to choose from. Thus £y”’ is a function 
of W=N!/N,!N2!---, and hence of log W, or of 
—N>n, log n,. With negligible action at a 


distance, Ey’ must be proportional to N for large 
N and for given n,’s, hence E,’"’=AN>n, log n,, 


or Vo’ =e>-n, log n,. Upon minimizing V' + V;" 


with respect to the m,’s we obtain the same 
formulas as before, with Lo replaced by 1/e. In 
the demagnetized state each n,=1/k, where k is 
the number of directions of easy magnetization, 
and V,’’=—elog k; when the field is large one 
n,=1 and the others vanish, and V’’=0. Hence 
the work done against internal forces in mag- 
netizing unit volume of the specimen, or /J/dJ, 
is € log k. 
(b) Boundary displacement against random 
forces. 

In the following discussion we abandon the 

somewhat arbitrary replacement of variable 


"Cf. J. Frenkel, Wave Mechanics: Elementary Theory 
(1932), p. 199. 
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yolumes by fixed volumes. We assume that a 
fraction v, of the total volume has J, in direction 
¢, and that each of the quantities 7, v2, --- isa 
function of the quantities 1), V2, --- 

Consider the ratio v;/v;. When V;= Vj, this is 
1, whatever the values of the other V,’s. When lV’; 
js much larger algebraically than V7, v;/v; is 
practically zero; and when V’; is much smaller 
algebraically than V’;, v;/v; is very large, what- 
ever the values of the other V,’s. Thus when 
V;-V;=— ©, 0, or +, v;/v; has a definite 
value, independent of the other V,’s and also of 
the absolute value of V’; or V; (with respect to 
any zero). Let us assume that in general v;/v; isa 
function only of V;—Vj, say v;/v;=f(Vi-—V)). 
Then for any three directions 7, 7, k we have 


1=(v;/v;) 
Vi), 
or if ¢=log f, 
o(x—y) +o(y—2) +o(s—x) =0. 
Differentiating with respect to z and then setting 


s=0 gives ¢’(y) = —¢’(—x) =constant = —A, or 
f(y) =Be~4". When y=0, f=1, therefore B=1. 


or finally 
=Ce4"1 etc., where C is determined by 


This derivation fails when there are only two 
directions of easy magnetization, but this case 
may be regarded as a limiting case attained by 
letting the value of V, for the other directions 
become very large. 

It is not obvious that the assumption made in 
this argument is equivalent to the assumption 


that the internal forces are random. The deri- 
vation is therefore incomplete; but it does 
suggest that the introduction of fixed instead of 
variable volumes, and the resulting superficial 
analogy to statistical mechanics, are probably 
not essential to the theory. 


Note added in proof: The following additional 
results of the theory are interesting because of 
their relation to some recently published work 
of Kondorsky.” In an iron crystal, when the 
field is along a specimen axis for which the 
demagnetizing factor is negligible, it may be 
assumed that only the [100] directions nearest 
this axis are occupied. Let the cosine of the acute 
angle between such a direction and the specimen 
axis be /, and let j=J,/J;,, n»=LoJ,//,: then the 
magnetization curve of the equivalent reversible 
specimen is given by j=/tanh yl, so that 
Xr = sech? yl = 2P(1—j?/P). 
It follows that for three crystals of the same 
material with their specimen axes along [100], 
[110], and [111] (?=1, 4, and }, respectively), 
and with the field in each case along the specimen 
axis, the initial susceptibilities should be in the 
ratios 6: 3:2, and x,/xo should be given by 
1—2j?, and 1—3j?, respectively. The 
ratios 6:3:2 were observed in Williams’ * 
measurements and are predicted by Kondorsky’'s 
theory when «;>xe, and also by earlier theories 
of Akulov and of Bozorth.” The expressions for 
xr/xo are identical with those of Kondorsky’s 
theory. 


2 E. Kondorsky, Phys. Rev. 53, 319 and 1022 (1938). 
2H. J. Williams, Phys. Rev. 52, 1004 (1937). 
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Measurements of Magnetic Viscosity in Iron 
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| Magnetic viscosit y—the change of induction occurring after the magnetizing field has ceased 
| to vary—has been measured in a bar of commercial iron, with a measuring circuit designed to 
eliminate spurious effects arising from sparking at switches. Magnetic viscosity is found to 
| depend on previous magnetic states, in agreement with observations of Mitkevitch and con- 
; trary to Preisach’s results. It appears, therefore, as if reversible domains could retain their 
lagging propensity while the magnetic force is varying over a considerable range. The viscous 
effect is found to depend on the magnitude of the previous change of induction if the latter is 
small, but is independent of this change if it is greater than 60 gauss. In small, subsidiary 
hysteresis loops, such as are used in determining reversible permeability, it is found that the 
Barkhausen effect is not present although there is hysteresis and magnetic viscosity. It is 
therefore concluded that the movement of the boundaries of saturated domains in the material 


is subject to time-lag. 


HEN a ferromagnetic material is subjected 
to a magnetizing field and this field is 
suddenly changed, the magnetic induction in the 
material will, in general, continue to change for 
an appreciable time after the magnetizing field 
has ceased to vary. Eddy currents in the specimen 
are responsible for part of this time-lag in 
magnetization. By a method given by Rayleigh 
it is possible to calculate the magnitude of this 
part of the phenomenon. In addition to the effect 
arising from eddy currents there is another time- 
lag which may persist long after the eddy currents 
have died away. The cause of this additional 
time-lag of induction, usually called magnetic 
viscosity, is not clearly understood at present. In 
the experiments described below measurements 
of magnetic viscosity in commercial iron wire 
have been made, and certain factors which 
modify this viscosity have been investigated. 


EXPERIMENTAL TECHNIQUE 


Figure 1 gives a diagram of the apparatus used. 
A pendulum closes switch K,, and_ shortly 
afterwards, Ke. When K;, is closed the current 
through the solenoid C is short circuited and the 
field in C falls practically to zero. The induction 
in the iron specimen S continues to change after 
K, closes, and with the closing of K» the galva- 
nometer G begins to deflect as a result of the 
induced e.m.f. in the pick-up coil P. This 
deflection is quickly brought back to zero by 


1Lord Rayleigh, Rep. Brit. Assoc. 52, 446 (1882). 


closing the switch K;, thus adding through the 
mutual inductance M enough flux in the galva- 
nometer circuit to balance the flux removed from 
P. As the magnetic induction of S continues to 
decrease the rheostat R; is varied by hand so as to 
keep the galvanometer deflection zero. From the 
final reading of the ammeter A? (after no more 
induction decrease can be observed), the number 
of turns in P, the mutual inductance M, and the 
cross-sectional area of S, it is possible to calculate 
the change of induction which occurred in § 
after the closing of the switch Ke. 

This method is superior to the ballistic method 
because the inductive changes may be followed 
over a long period of time. Another advantage of 
the circuit shown in Fig. 1 is that it eliminates 
uncertainties connected with the spark produced 
when a switch is opened in an inductive circuit. 
An extended series of measurements was made 
with a circuit in which the solenoid current was 
reduced to zero by opening a switch in series 
with the battery. It was found that when proper 
values of resistance and capacity were used as 
shunt across the solenoid so as to produce 


ps 


Fic. 1, Circuit diagram. 
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theoretically a nonoscillatory decay of current in 
the solenoid, nevertheless a high frequency 
oscillation of short duration occurred at the 
instant of break. This oscillation was detected 
with a cathode-ray oscillograph connected across 
a noninductive resistance in series with the 
solenoid. Under certain conditions no magnetic 
viscosity whatever could be observed, probably 
because of the spark oscillations superposed on 
the unidirectional decay current. Under other 
conditions of resistance and capacity these 
superposed spark oscillations did not appear to 
affect the viscosity. Pending a more complete 
investigation of this phenomenon it was decided 
to eliminate uncertainties by the use of the short- 
circuiting device of Fig. 1. There is a very definite 
probability that some of the conflicting results of 
previous experimenters may be due to the 
unsuspected presence of spark oscillations in their 
circuits. 

Further details of the circuit are as follows: 
The galvanometer was of the high sensitivity 
type, of period 19 seconds, used with a scale 
distance of 470 cm. It was necessary to insulate 
the circuits with considerable care in order to 
avoid spurious effects. The battery E was a 6 
volt storage battery. The solenoid C, wound on a 
micarta tube, was 30.6 cm long, contained 7000 
turns, and had a resistance of 42.5 ohms. Its 
mean diameter was 9.26 cm. The coil P had 
10,000 turns and was supported in the center of 
the solenoid. The specimen S was centered with 
respect to C and P. It was a wire of commercial 
iron, 0.30 cm in diameter and 15.1 cm long. It was 
annealed by slow cooling from a bright red heat. 
Figure 2 gives the descending branch of the 
hysteresis loop of this specimen, without cor- 
rection for the factor of demagnetization. 

In order for the eddy currents to have time to 
die away in a specimen of this diameter it is 
necessary to have a fairly large time interval 
between the closing of switches K; and Ke. In 
most of the work the switches were set so that 
this time interval, by calculation, was 0.61 
seconds. During this time the current in the 
solenoid must fall to a negligible value and the 
eddy currents in the specimen must die away. If 
R; in Fig. 1 is taken to be 100 ohms, the value 
usually used, a simple calculation shows that 
0.13 second after the closing of K, the current in 
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C is reduced by the factor 1.1X10-*. The 
maximum current used in the experiments was 
0.125 ampere. Thus after 0.13 second the current 
has fallen to a value 1.4X10-? ampere. Tests 
showed that with the specimen in the solenoid the 
induction changes caused by a current of 6X 107-7 
ampere produced no detectable effect on the 
galvanometer. We may therefore assume the 
solenoid current to be zero 0.13 second after the 
closing of K,. The eddy currents—and therefore 


Fic. 2. Descending branch of the hysteresis loop of a 
commercial iron wire which had been annealed by slow 
cooling from red heat. No correction has been made for 
demagnetization. 


the magnetic induction—in the specimen will 
require some time to die away after the solenoid 
current has become zero. Rayleigh’s formula for 
the decay of eddy currents may be put in the 
form log b= —1.44t/(Acu), where d is the fraction 
of the total induction change which occurs after a 
time ¢, and ¢ is the time after the magnetic field is 
removed; o and A are the conductivity and 
cross-sectional area, respectively, of the speci- 
men. The value of the permeability u is set equal 
to dB/dH and can be obtained from the B-J/ 
curve of Fig. 2 after the appropriate factor of 
demagnetization has been applied. By using the 
least favorable value of dB/dH and _ taking 
t=0.48 second, which is the time remaining after 
the solenoid current becomes effectively zero 
before Ke is closed, b is calculated to be 2.6 X 10~°. 
The maximum value of the total induction 
change in the experiments where viscosity was 
measured was 6400. Thus we obtain for the 
induction change (6B = 64000) which occurs after 
Kz is closed, and which is caused by eddy 
currents, the value 0.017. The smallest induction 


| 
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t sec. 


Fic. 3. The magnetic viscosity of iron as a function of the 
time it was kept at the point D of the hysteresis loop. 


change which could be detected with certainty by 
the apparatus was about 0.2 gauss. Thus any 
effect arising from eddy currents is too small to be 
detected in any of the present experiments. 


INFLUENCE OF PREVIOUS MAGNETIC STATES 


The small hysteresis loop shown in Fig. 3 was 
traversed several times so as to stabilize the loop. 
The viscosity at E following the drop from D was 
then measured and found to depend on the 
length of time the iron was kept at point D. The 
experimental points on the large curve of Fig. 3 
were obtained as follows. The iron was left at 
state A of the small loop for about one minute—a 
sufficient time for all appreciable viscous effects 
to disappear—then by a reversing switch carried 
to state D. After a time ¢ the pendulum was 
allowed to close key K,, and 0.61 second later, 
K;. Thus the iron was carried to state E and the 
viscous effect, B,, measured. In order to measure 
t a noninductive resistance of 40 ohms was 
connected in series with the solenoid and po- 
tential leads from the ends of this resistance were 
led to the filament and grid, respectively, of a 
triode amplifying tube. When the reversing 
switch was thrown the potential change across 
the 40 ohms caused a change of plate current and 
thus actuated, through a sensitive relay, the 
stylus of a revolving drum chronograph. When 
the key K;, closed the stylus was again actuated. 
In this way a sufficiently accurate determination 
of ¢ could be made. 

Figure 3 shows that the magnetic viscosity, B,, 
at E is diminished when ¢ is reduced below about 
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3.5 seconds. For larger values of ¢ it appears that 
B, is constant as ¢ changes. In interpreting the 
curve it is to be noted that when ¢ is less than 
about 0.3 second the iron at D may still be the 
seat of eddy currents of appreciable magnitude at 
the instant the key K, is closed. However, it js 
quite evident that after all eddy currents have 
died away there is still a magnetic readjustment 
going on at D which must be completed if the fy! 
value of the viscosity at E is to be obtained, 
When the iron is carried to D we may say that 
the magnetic viscosity at that point consists jn 
the reversal of lagging magnetic domains, so that 
the induction at D continues to increase for some 
time after point D is reached by the magnetic 
field. If all of these lagging domains are not 
allowed to reverse before the field is reduced to 
the value at E we do not find so big a viscous 
effect at E. The decrease of viscosity at E is due 
to the fact that certain members of a group of | 
viscous domains have not had time to reverse at 
D, hence these members are oriented in the right 
direction at E and do not need to reverse under 
the influence of the demagnetizing field. 

Let Np and Nz be the numbers of domains 
which show magnetic viscosity in reversing at the 
respective points D and E. Assume (1) that 
members of the group Np which do not have 
time to reverse their magnetism at D are 
detracted from the viscous group Ng at E; 
(2) that the number of unreversed viscous 
domains at D varies exponentially with the time { 
after the field becomes constant. That is, let 
N=Npe-™, where N is the number of domains 
still unreversed ¢ seconds after the field has 
ceased to vary, and m isa constant. 

The measured magnetic viscosity at E is then 
given by B,=k(Ne—N)=k(Ne— Nope) =h, 
—kse-™', where k, k;, and ke are constants. The 
curve of Fig. 3 is drawn with k,=5.75, ke=4.25, 
and m=0.75 in the above equation. The experi- 
mental points fit the curve as well as could be 
expected. 

A consequence of the above equation for B, is 
that when point D is chosen far out on the 
saturated part of the magnetization curve the 
value of B, should be independent of ¢; for in this 
case the viscosity at D is found experimentally to 
be very small, so that Np may be set equal to 
zero and thus B,=k,. The constant hk, in this 
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case, however, may be different from the &, for 
other positions of D, because, in general, Ne 
depends on the previous magnetic history of the 
specimen. Experimental tests showed that when 
point D is taken to a field of 34 gauss the value of 
B, at Cis about the same, whether ¢ is 0.3 second 
or one minute. 

In case Np is larger than Nx the value of B, 
would be negative for small values of ¢. The 
specimen of iron used in this experiment did not 
give negative values of B, for the part of the 
hysteresis loop examined. Previous experiments 
disagree to some extent in this matter. Mitke- 
vitch? and Richter* were able to observe the 
negative effect predicted by the above equation. 
Preisach* not only failed to report the negative 
effect but also he states that the time during 
which a specimen is subjected to a field has no 
effect on the viscosity appearing when the field is 
cut off. Preisach states that the ‘‘superposition 
principle” is not applicable to magnetic viscosity. 
The conflict between the results of Preisach and 
those of Mitkevitch, Richter, and the writer 
could very well be due to the presence of spark 
oscillations in the solenoid of Preisach. He does 
not give a diagram of his circuit nor does 
he mention any precautions taken to stop 
oscillations. 


RELATION BETWEEN THE MAGNITUDE OF AN 
INDUCTION CHANGE AND THE ENSUING 
MAGNETIC VISCOSITY 


The procedure in this experiment was as 
follows. From a saturating field of 35 gauss the 
field was decreased to a definite small value, 77’, 


B, 
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Fic. 4. Magnetic viscosity of iron as a function of the 
induction decrease AB, 


2A. Mitkevitch, J. de phys. 7, 133 (1936). 
°G. Richter, Ann. d. Physik 29, 605 (1937). 
‘F. Preisach, Zeits. f. Physik 94, 277 (1935). 


VISCOSITY 291 


and kept constant at that value for one minute. 
Then the pendulum was allowed to operate the 
switches K, and Ke, thus removing the field 17’, 
producing the induction decrease AB, and 
allowing the resulting magnetic viscosity B, to be 
measured. Curve I of Fig. 4, obtained in this 
way, shows that B, is constant for large values of 
AB but diminishes rapidly when AB is decreased 
below about 60 gauss. 

If the time of holding /7’ constant is diminished 
from one minute to 5 seconds curve II is ob- 
tained. Curves I and II coincide for large values 
of AB but below 50 gauss curve II is the higher. 

These curves show that the number of viscous 
domains which reverse after the induction jump 
AB has occurred is practically independent of the 
magnitude of AB when the latter is large, but as 
AB becomes small this number decreases. Curve 
II is higher than I because for the former curve 
some of the viscous domains which had no time 
to reverse during the pause at H/’ remain 
unreversed during the change AB and thus 
contribute to the viscosity after AB has occurred. 
The curve II does not pass through the origin 
because when the saturating field of 35 gauss is 
brought to zero for 5 seconds and no further 
change made, i.e., AB=0, then even after 5 
seconds there are still numerous domains which 
have not yet had time to reverse. These domains 
give the viscosity measured when Kz is closed at 
the end of the 5 seconds. 


RELATION BETWEEN MAGNETIC VISCOSITY AND 
BARKHAUSEN DISCONTINUITIES OF 
MAGNETIZATION 


In the previous experiments if the galvanometer 
is replaced by an amplifier and telephone re- 
ceivers it is possible to detect Barkhausen 
discontinuities which persist and can be heard in 
decreasing numbers for some time after the 
magnetic field has ceased to vary. In line with 
this observation, therefore, it is usually con- 
sidered that the entire phenomenon of magnetic 
viscosity is due to the existence of delayed 
Barkhausen reversals of magnetization. The 
cause of this delay in the reversal of Barkhausen 
domains has been attributed® to persisting eddy 
currents and to the resulting microscopic strains 


5H. Kuhlewein, Physik. Zeits. 32, 860 (1931); C. W, 
Heaps, Phys. Rev. 49, 409 (1936). 
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set up by the local inequalities of temperature. 
The experiments described below indicate that 
magnetic viscosity is not to be explained entirely 
in terms of these delayed Barkhausen reversals. 

A small subsidiary hysteresis loop similar to 
that shown in the inset of Fig. 3 was examined, 
with an amplifier and telephone receivers, for a 
lagging Barkhausen effect at the points A, C, D, 
and E. At these same points the magnetic 
viscosity was also measured with the galvanome- 
ter as described above. It was found for a loop of 
large area that Barkhausen jumps and magnetic 
viscosity were both present at all the specified 
points. When the loop was made smaller and 
stabilized by repeated traversing of the cycle the 
Barkhausen jumps disappeared at certain of the 
points. 

For example, when the limits of variation of B 
and // were, respectively, 390 and 1.6 gauss no 
Barkhausen effect could be detected at any 
of the points. The magnetic viscosity, however, 
amounted to 1.0 gauss at each point. This figure 
represents, as before, the total induction change 
occurring 0.61 second after the solenoid current is 
switched on or off. For points A and D this time 
is only approximate, as here the key K, could not 
be used and it was necessary to close the reversing 
switch, and subsequently Ke, by hand. In 
observing the Barkhausen effect the amplifier 
system which had been substituted for G was 
kept connected at all times, i.e., Ke was kept 
closed. Thus when the reversing switch was 
operated there was always a loud, sharp noise 
caused by the big inductive jump. This noise was 
followed, in case of lagging Barkhausen jumps, 
by a crackling, or popping noise which gradually 
died away as the individual impulses became less 
frequent. 

For the loop of limits specified above this 
lagging Barkhausen effect was at first large at 
point A. After about five circuits of the loop it 
vanished at this point. The other points of the 
loop showed no lag at any time. The magnetic 
viscosity as measured with the galvanometer 
showed a similar initial large effect at A— 
approximately 9 gauss—which diminished to the 
constant value 1.0 gauss after about eight 
circuits of the loop. 

It seems improbable that this viscous effect of 
1.0 gauss is due to Barkhausen jumps too small to 


HEAPS 


be heard in the telephone receivers. As the 
Barkhausen effect diminishes with successive 
trips around the loop it is not the intensity of the 
individual impulses which seems to weaken, it js 
their number. It appears, therefore, as if the 
members of a definite group of lagging, reversible 
domains were gradually removed from. this 
lagging group by repeated trips around the 
loop. If the magnetic viscosity which remains 
after several cycles is due to discontinuities of 
magnetization these discontinuities must be of a 
different order of magnitude from those ordinarily 
observed because none whatever could be 
detected even when four stages of amplification 
were used. 

Small loops of the kind here described are used 
in determining reversible susceptibility. The area 
of these loops may be made very small by 
decreasing the amplitude of the field variation 
but the statement which is sometimes made‘ that 
these small cyclic processes involve no dissipation 
of energy appears to be incorrect according to the 
best experimental evidence.’ It is only in the 
limit, as the amplitude of the field variation 
approaches zero, that the area becomes zero. For 
the loop as used above the width along the B axis 
was about 45 gauss. Accordingly the disappear- 
ance of lagging Barkhausen jumps is not to be 
attributed to the disappearance of hysteresis. 

The question arises as to whether the Bark- 
hausen discontinuities merely fail to lag behind 
the field or disappear for all points around the 
loop. An experimental test of this question was 
made as follows. The field was caused to vary 
slowly and continuously around the loop and 
simultaneously a check was made, with the use of 
the telephone receivers, to see whether Bark- 
hausen discontinuities occurred. The slow con- 
tinuous field variation was obtained by passing 
the plate current of a radio tube, type 89, 
through the solenoid. When the heater current 
of the tube was switched on or off the plate 
current rose or fell smoothly and slowly. After 
several cycles had stabilized the loop no Bark- 
hausen jumps could be detected at any point. 
Accordingly we conclude that in these small 


‘°F. Bitter, Introduction to Ferromagnetism (McGraw- 


Hill Book Co., 1937), p. 182. 
7E. Spuhrmann, Zeits. f. Physik 39, 332 (1926); K. 


Uller, Zeits. f. Physik 38, 72 (1926); G. J. Sizoo, Ann. d. 
Physik 3, 270 (1929). 
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subsidiary loops there is hysteresis and magnetic 
viscosity without any Barkhausen effect. 

According to the domain theory of ferro- 
magnetism a material may change its intensity of 
magnetization by three different processes: 
(1) reversal of the saturated magnetization of 
small domains, (2) reorientation of the mag- 
netization of the domains, (3) increase of size of 
one domain at the expense of another. The first 
process is supposed to produce the sudden 
inductive jumps associated with the Barkhausen 
effect. The last two processes may take place 
slowly and continuously, so would not be 
detected by the use of telephone receivers. It 
appears, therefore, in view of the experiments 
just described, that the phenomenon of magnetic 
viscosity is involved in one or both of the two 
last named processes, as well as in the first. 

It is usually considered that process (2) occurs 


in strong fields, process (3) in weak fields. For 
the loop used in the present experiment there is 
always a demagnetizing field present ; however, it 
seems reasonable to suppose that process (3) is 
largely responsible for the induction changes of 


these small loops. 


In this case the conclusion to be drawn is that 
when the boundaries of domains move under the 
influence of applied fields a certain time is 
required for equilibrium to be attained after the 
magnetic field ceases to change. These moving 
boundaries would be the seat of eddy currents 
and magnetostrictive strains, so that magnetic 
viscosity would make its appearance in process 
(3) as a result of the time required for the 
disappearance of microscopic temperature gradi- 
ents and the readjustment of the local strains. 


’F. Bloch, Zeits. f. Physik 74, 333 (1932); R. Becker, 
Physik. Zeits. 33, 905 (1932). 
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Methods are described for growing single crystals of sodium in the form of rods 4.7 mm in 
diameter and 10 cm long, for handling the crystals in an atmosphere of helium, and for meas- 
uring the principal elastic moduli at low temperatures. Values of the adiabatic and isothermal 
moduli and constants are tabulated at ten degree intervals between 80°K and 210°K. The 
value of the Debye characteristic temperature, calculated from the values of the elastic con- 


stants at 80°K, is 164°K. 


EXPERIMENTAL METHOD FOR MEASURING THE 
ELastic CONSTANTS 


COMPLETE description of the dynamical 

method employed in this research has ap- 
peared in previous issues of this journal.' Ac- 
cordingly, it will suffice here briefly to review its 
essential features, and their adaptation to the 
present experimental problem. The specimen is 
in the form of a right circular cylinder 4.7 mm in 
diameter and a few centimeters long. The re- 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1 Balamuth, Phys. Rev. 45, 715 (1934); Rose, Phys. Rev. 
49, 50 (1936); Durand, Phys. Rev. 50, 449 (1936). 


quired data are deduced from the observed 


‘behavior of a separately excited composite 


piezoelectric oscillator constructed by cementing 
to one end of the specimen a suitably cut cylinder 
of crystalline quartz of identical cross-section. 
Silver electrodes are chemically deposited in 
proper position on the quartz, and the oscillator 
is suspended vertically by delicate supports 
attached at the middle of the quartz cylinder. 
One or more harmonic frequencies of free longi- 
tudinal or torsional vibration of this system are 
measured by observing the variation of the 
electrical impedance of the composite oscillator 
with the frequency of the applied voltage, and 
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from these the fundamental frequencies, f, of free 
longitudinal and torsional vibration of the speci- 
men cylinder alone are calculated. The latter 
quantities are related, respectively, to the 
Young’s modulus, £, and the torsion modulus, G, 
of the specimen material for the direction of the 
cylinder axis by the formula 


2fL=(M/p)}, (1) 


where L is the length of the specimen cylinder, p 
is the density, and M the elastic modulus. It 
remains only to relate EZ and G to the principal 
elastic moduli and elastic constants of the 
crystal. 

Young’s modulus for an arbitrary direction in 
a cubic crystal is related to the principal elastic 
moduli by the formula? 


1/E=s,,—2sT, (2) 

S=S11—Si2— 3544, 

P= 

and a, 8, y are the direction cosines with respect 
to the principal crystal axes. 

When torsional surface tractions are applied 
across the end faces of a crystalline cylinder, the 
resultant twist is, in general, accompanied by a 
bend. This additional strain can be, and in static 
measurements usually is, prevented by additional 
contraints, but the relation between the torsion 
modulus and the principal elastic moduli is quite 
different according as this is or is not done. If 
bending is not prevented 


1/G=sut+4sT, (3) 
while if bending is prevented* 
1/G=sy+4sT (4) 
where x 

At first it was thought that the torsional 
vibration of the specimen cylinder would be 
governed by the formula (3), and the experi- 
mental results on the first three crystals grown 
were not inconsistent with this view. The results 


of subsequent experiments, however, showed 
conclusively that the high frequency‘ torsional 


where 


2 Voigt, Lehrbuch der Kristalphysik, p. 739. 
3 Goens and Schmidt, Naturwiss. 19, 520 (1931). 
4 The frequencies employed lie between 28 and 55 kilo- 


cycles. 


vibration of the specimen constituent of the 
composite oscillator is governed instead by the 
formula (4). The earlier error arose from the 
fortuitous circumstance that the cylinder axes of 
the three crystals lay very nearly in the (100), 
(110) and (111) directions, for which directions 
the two formulae are identical. 

The elastic moduli and elastic constants are 
related by the formulae, 


C11 = (Str +512)/ (S11 — S12) (S11 +2512), 
C12= —S12/(S1r— S12) (Sir +2512), (5) 
Cy = 1/S45. 


Lastly, the compressibility, x, is given by the 
formula, 
k=3(Si: +2512), 
=3/(€1+2¢12). (6) 


The method yields immediately the adiabatic 
moduli and constants. These are related to the 
isothermal quantities by the formulae® 


(Si1)aa. — (S11) is. = (Si2)aa. — (Si2)is.= — Ta*/ pcp, 
(C11) aa. — (C11) is. = (C12) aa. — (C12) is. (7) 


= Ta?/(pc,) (S11 +2512)’, 


where T is the absolute temperature, a is the 
coefficient of linear expansion, c, is the specific 
heat at constant pressure, and p is the density. 
The adiabatic and isothermal s4, and c4, are the 
same. 

To recapitulate: The essential measurements 
of the method are the fundamental frequencies of 
the specimen, its length and density, and the 
orientation of the cylinder axis in the crystal 
lattice. In addition, the specimen cylinder must 
be homogeneous and of uniform cross section. 


PREPARATION OF THE CRYSTALS 


The crystals are grown of Mallinckrodt 
analytic reagent grade sodium triply distilled in 
vacuum, with the rejection of the first and last 
quarters of the distillate in the first distillation. 
The still is constructed of Pyrex glass. Spectro- 
scopic analysis shows the final product to contain 
traces of -copper, vanadium, and aluminum, 
together with an amount of potassium which 
defies quantitative chemical analysis. The second 
distillation leaves the sodium in an ampoule, 
which is sealed off and placed in a Pyrex glass 


5 Voigt, reference 2, p. 789. 
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GS 
— STEEL BALL 


vA 


He CRUCIBLE 


HEATING COIL 


Fic. 1. Cross section of apparatus for growing 
the sodium crystals. 


tube connected with the upper chamber of the 
apparatus diagrammed in Fig. 1. This apparatus 
is evacuated and kept at about 125°C with an 
enveloping furnace which has a glass front. 

The ampule is broken in vacuum and the 
sodium distilled into the upper chamber. The 
bottom of this chamber is drawn out into a long 
tube of about 1 mm bore, which projects down- 
ward into a crucible and is closed at the top by a 
nickel plated steel ball. The crucible is made 
from a very thin walled glass tube, of the sort 
used for the stems of hydrometer bulbs, carefully 
selected for uniformity of bore.* The interior of 
the crucible is thinly coated with Socony Super 
Hecla mineral oil, previously boiled in vacuum in 
the presence of sodium. The crucible hangs in a 
nickel frame, suspended by a nickel chain which 
passes over pulleys to an iron counterweight. The 
position of the counterweight is controlled by a 
solenoid, which can be raised or lowered either 
manually or by a geared down electric motor. 

When the distillation into the upper chamber 
is complete, helium at about 2 mm pressure is 
admitted above the molten sodium and the steel 
ball is raised by an electromagnet. The crucible is 
lowered steadily, by hand, as it fills, so that the 
bottom of the filling tube remains only a few 
millimeters below the surface of the sodium. 
After the crucible is filled it is lowered electrically 


* The writers are indebted to the firm of Kessler and 
Mumberg, of New York, for ‘ene to make this 
selection directly from their stock. 


out of the furnace, past the auxiliary heating coil 
shown at the bottom of the figure, at the rate of 
5 cm per hour. Helium is then admitted to the 
lower chamber, the apparatus is opened at a 
ground joint level with the top of the lowered 


- crucible, and the crucible is quickly capped with 


soft wax in the emerging stream of helium. The 
yield of good crystals, single over the entire 10 cm 
of length, is about one in three. The crystals 
appear to have no tendency to grow in any 
preferred orientation.’ 


LocATION OF THE CYLINDER AXIS IN THE 
CrYstTAL LATTICE 


The orientation of the cylinder axis with 
respect to the principal crystal axes is obtained 
by the measurement of four Laue x-ray trans- 
mission photographs. The diffraction camera has 
a G. E. Co. model CA 2 molybdenum x-ray tube 
operated at 35 kv and 20 ma. The specimen-to- 
film distance is 4.86 cm. Dupont film backed by a 
Patterson fluorazure screen is used, and the 
exposure time is 10 minutes. 

The crystal, in the crucible, is set with its 
cylinder axis normal to the x-ray beam, and the 
pictures are taken after successive 90-degree 
rotations of the cylinder about its own axis. The 
spots are identified by comparison with the 
admirable set of Laue transmission photographs 
of a body-centered cubic lattice published by 
Majima.* The three direction cosines of each 
beam with respect to the crystal axes are 
calculated by the method of least squares from 
measurements of the radial distances of 8 to 12 
spots from the center spot. The specimen-to-film 
distance used in these calculations is obtained by 
analysis of a picture taken through a cylinder of 
crystalline NaCl, in which the cylinder axis 
coincides accurately with a principal axis. 

The four x-ray beams lie on the surface of a 
cone whose axis is the cylinder axis and the 
supplement of whose vertex angle is twice the 
error in setting the cylinder axis normal to the 
x-ray beam. If the latter were zero, the direction 
cosines of the cylinder axis would be those of the 

7This method may be compared with that devised 
independently by Andrade and Tsien, ‘‘The Glide of 
Single Crystals of Sodium and Potassium,”’ Proc. Roy. 
Soc. 163, 1 (1937). 


8 Majima, Scientific Papers of the Institute of Physical 
and Chemical Research, Vol. 7, p. 249. 
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Fic. 2. Cross section of the chamber in which the crystals 
are handled. 


vector cross product of any two mutually 
perpendicular x-ray beams. The values which are 
adopted are the arithmetic means of the four 
corresponding direction cosines obtained from the 
four possible vector cross products among the 
four x-ray beams. 

The writers were rather surprised to discover 
that the orientation of the cylinder axis can be 
obtained in this manner with an accuracy of 
about five minutes of arc. 


HANDLING THE CRYSTALS 


After the x-ray analysis the crystal is placed in 
the handling chamber diagrammed in Fig. 2.’The 
chamber, which is cylindrical in form and built of 
brass, is evacuated and filled with dry helium at 
atmospheric pressure. The crystals are handled 
through opposite holes in the chamber (shown by 
the dotted lines of Fig. 2), with the hands and 
arms encased in long obstetrical rubber gloves 
whose openings are securely clamped about the 
peripheries of the holes. The auxiliary chamber, 
L, is a lock which permits the introduction of 
crystals and tools to the chamber with a mini- 
mum expenditure of helium. 

The crystal is removed from the crucible and 
cleaned with anhydrous ether kept standing in 
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the presence of sodium. It is then laid in ag 
V block and cut to the desired length with a thin 
strip of nickel under tension. The strip moves in a 
transverse slot in the V_ block. The quartz 
cylinder is prepared by coating one end with a 
very thin film of Empire stopcock lubricant, 
which has been boiled in vacuum immediately 
before use. The quartz and specimen cylinders 
are laid end to end in the V block, and the join js 
made with the specimen held gently in place, by 
imparting a forward rotatory motion to the 
quartz. 

The quartz cylinder has minute diametrically 
opposed transverse nicks at its center, and the 
composite oscillator is mounted by engaging 
these nicks with the jaws of a clip, shown at C, 
Fig. 2. The jaws are constructed by stretching 
fine wires horizontally across the forked lower 
termini of adjacent light vertical strips of spring 
metal. Two such strips are necessary for a 
longitudinal and four for a torsional oscillator. 

The oscillator is enclosed by joining the glass 
tubes AA (Fig. 2), and the mount lock is closed 
by joining the glass tube B, to the brass tube Bs, 
Lastly the envelope containing the mounted 
oscillator is disengaged from the handling 
chamber at the point W, where it rests on a 
rubber washer, is removed and inserted in a 
dewar of boiling liquid nitrogen or liquid oxygen, 
or in the cryostat described in a_ previous 


paper.”® 


CALCULATION OF THE ELAstTic MopuLI AND 
CONSTANTS 


The quantities which follow immediately from 
the measurements made in the cryostat are the 


TABLE I. Typical data on a single sodium crystal. 


QUARTZ VELocITY 
No. n 10-5 cm/sec. 

I 5 11004 1.443 
L II 3 13107 1.468 
I 3 18187 1.472 
II 2 18263 1.478 
T III 2 20237 2.267 
IV 1 27837 2.253 


* Empire Laboratory Supply Co., New York. 

Siegel and Quimby, Phys. Rev. 54, 76 (1938). The 
temperature of this cryostat is measured with a copper- 
constantan thermocouple, calibrated by comparison with 
a platinum resistance thermometer certificated by the 
National Bureau of Standards, 
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frequencies of longitudinal and torsional vibra- 
tion of the specimen crystal cylinder. The length 
of the cylinder at room temperature is measured 
on a comparator. The value adopted for the 
density of sodium at 20°C is 0.9712." The lengths 


and densities at other temperatures are calculated . 


from these numbers and the data on the thermal 
expansion of crystalline sodium previously pub- 
lished in this journal.’ Corresponding values of 
E and G are computed with Eq. (1). 

Twelve crystals in all were examined in the 
course of the present research. Of these, six were 
investigated very thoroughly at the temperature 
79°K, and contributed twenty-four values of 


and G at this temperature. Where possible - 


measurements were made with two different 
longitudinal and two different torsional quartz 
cylinders, and at two or more different harmonic 
frequencies. A typical set of data is given in 
Table I. The first column designates the quartz 
cylinder by number, the second contains the 
harmonic numbers, i.e., the approximate number 
of half-waves of vibration in the specimen 

TABLE II. Observed and computed values of the longi- 


tudinal and torsional velocities in cylinders of crystal sodium 
at 79°K. 


LONGITUDINAL VELOCITY TORSIONAL VELOCITY 
cm/sec. cm, sec. 

r OBSERVED ComPUTED OBSERVED COMPUTED 
0.0004 1.399 1.440 2.364 2.423 
0.0130 1.465* 1.457 2.260* 2.342 
0.1175 1.793 1.833 
0.1476 1.795* 1.807 1.656* 1.680 
0.1488 1.766 1.812 1.678 1.708 
0.2128 2.104* 2.102 1.406* 1.389 
0.2286 2.146 2.198 1.259 1.257 
0.2491 1.148 1.149 
0.2661 2.472 2.492 1.132 1.131 
0.2907 2.780* 2.765 1.187* 1.181 
0.3044 2.962 2.962 1.073 1.089 
0.3145 3.206 3.137 


* Indicates the mean of several observations. 


cylinder, the third the measured fundamental 
vibration frequency, and the last the calculated 
velocity of propagation, V(=2fL). 

The remaining six crystals contributed ten 
values of E and G at 79°K,” and all of the data 


" Richards and Brink, J. Am. Chem. Soc. 29, 117 (1907). 

2 Instead of twelve. The occurrence of a parasitic mode 
of vibration, feebly excited, vitiated the torsional measure- 
ments on one of these crystals, and another was, after the 
removal of a flaw, too short to permit reliable measurement 
of E. One of the first six crystals was dropped before E 
was measured. 


+2.00 


ri90 


FREQUENCY SQUARED 


+185 TEMPERATURE °K 
80 100 120 140 160 


Fic. 3. The variation of f? with temperature for the 
torsional vibration of one of the crystals. 


which show the variation of the elastic moduli 
with temperature. The number is adequate, for 
it is characteristic of the present experimental 
method that the temperature variation of the 
several elastic moduli can be measured with 
considerably greater accuracy than their abso- 
lute values. Imperfections in the specimen and 
other sources of error affect the former measure- 
ments less than the latter. 

The adiabatic elastic moduli at 79°K were 
evaluated by solving thirty-four equations of the 
forms (2) and (4) simultaneously, by the method 
of least squares, for the three unknown quantities 
S11, Sig and S44. The concordance of the data on 
all twelve crystals is exhibited in Table II, in 
which the observed propagation velocities may 
be compared with those computed from the 
values of the elastic moduli obtained in this 
manner. The maximum deviation is 3.6 percent 
and the average deviation is 1.3 percent. It will 
be noted that, since the value of the orientation 
factor, ', must lie between 0 and 3, the range of 
possible values is well covered by the twelve 
crystals. 

The variation of the elastic constants with 
temperature is obtained as follows: The squares 
of the observed frequencies, f, of longitudinal and 
torsional vibration are plotted as functions of the 
temperature. The observations on the six crystals 
are thus represented by ten curves, of which that 
shown in Fig. 3 is typical. Values of f? at ten 
different temperatures between 80°K and 210°K 
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i 


Fic. 4. The variation of the adiabatic elastic moduli of 
sodium with temperature. 


are read from each curve, and corresponding 
values of E and G are calculated with Eq. (1). 
These data yield ten sets of ten equations of the 
forms (2) and (4), each set corresponding to a 
single temperature. Each set of equations is 
solved simultaneously, by the method of least 
squares, for S11, Siz and S44. The result of this 
computation is shown by the curves of Fig. 4. 
The absolute values of the elastic moduli at 
79°K read from the curves of Fig. 4 differ slightly 
from those previously computed, since the 
former are based on ten observations and the 
latter on thirty-four. Thus the values of si1, — S12 
and s44 at 80°K deduced from all observations 
are, in units of 10-" cm?/dyne, 4.83, 2.09 and 
1.69, respectively, while those deduced from the 
ten observations alone are 4.91, 1.98 and 1.69. 
However, these curves are used simply to 
determine the ratios of the values of the moduli at 
T°K to the values at 79°K. The values hereinafter 
recorded for the elastic moduli at T°K are 
obtained by applying the ratios obtained from 
the curves of Fig. 4 to the values at 79°K 
deduced from the thirty-four observations at 
79°K in the manner previously described. The 
concordance of the data from which the ratios 
are evaluated is exhibited in Table III, in which 
the observed ratios of the velocities at 80°K to 
those at 210°K for the six crystals may be 
compared with the ratios computed with the 
values of the elastic moduli at 80°K and 210°K, 


thus obtained. The results confirm the remark 
made above, that the present method yields a 
more accurate determination of the temperature 
variation of the elastic moduli than of their 
absolute values. The maximum deviation jp 
Table III is 1.5 percent and the average deviation 
is 0.4 percent. 

The isothermal moduli and constants are 
computed from the adiabatic quantities with 
Eqs. (7), in which, since the difference is small, 


TABLE III. The observed and computed ratios of the 
longitudinal and torsional velocities at 80°K to those at 210°K 
for six crystals. 


LONGITUDINAL TORSIONAL 
OBSERVED COMPUTED OBSERVED ComPUTED 
1.073 1.065 1.060 1.060 
1.052 1.053 1.086 1.101 

1.072 1.075 

1.061 1.061 1.060 1.064 
1.061 1.059 

1.070 1.071 1.086 1.078 


the values p= 1.0, cp>=1.1X 107 and a=6.0X 10-5 
are substituted. 


RESULTS 


The elastic moduli and constants 


_ The adiabatic and isothermal elastic moduli 
and constants, and the isothermal compressibility, 
are given as functions of the temperature in 
Tables IV and V. The estimated uncertainty in 
the absolute values of 51; and sq, is about 2 
percent and of s;2 about 3 percent. The precision 
with which the tabulated values are given is 
significant only for the calculation of the temper- 
ature variation. 

Fuchs" has given a quantum-theoretical calcu- 
lation of the quantities (¢i:—¢i2) and cs, for 
sodium at 0°K, which yields the values 0.141 and 
0.58, respectively, in units of 10" dynes/cm?. The 
corresponding experimental values at 80°K are 
0.145 and 0.59. 


The compressibility 


The accuracy of the present measurement of 
the compressibility is very poor. A significant 
figure is lost in the evaluation of the difference 


13 Fuchs, Proc. Roy. Soc. 157, 444 (1937). 
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quantity (si+2s12), and the uncertainty in $12 is 
doubled. Thus a change of one percent in the 
assumed values of si; and $)2 will, if the errors are 
additive, change the calculated value of « from 
1.98K10-" to 1.70X10-". The writers have 
little doubt that the values given in Table V are 
high. 

Bridgman" has measured the variation of the 
compressibility of polycrystalline sodium with 
pressure between 2000 and 20,000 kg ‘cm’. As the 
pressure is lowered the compressibility increases 
toward the (extrapolated) value 1.710-" 
cm?/dyne at zero pressure and room tempera- 
ture. If the linear variation of compressibility 
with temperature here found to hold between 
g0°K and 210°K be assumed to hold over the 
entire range from 0°K to 293°K, then Bridgman’s 
value reduces at absolute zero to 1.210-", and 
the value here found to 1.8X10-". A quantum- 
theoretical calculation of the compressibility at 
0°K" gives k= 1.6 10-". 


TaBLE IV. The adiabatic elastic moduli and elastic con- 
stants of sodium. The number of significant figures does not 
indicate the accuracy of the absolute value. 


ADIABATIC CONSTANTS 
dyne, cm? 


ADIABATIC Moputti 
X10" cm? dyne 


T°K su —Si2 Sua C12 cus 


80 | 4.826 | 2.087 | 1.685 | 0.608 | 0.463 | 0.593 

90 | 4.861 | 2.102 | 1.705 | 0.603 | 0.459 | 0.586 
100 | 4.898 | 2.118 | 1.726 | 0.598 | 0.456 | 0.579 
110 | 4.937 | 2.134 | 1.747 | 0.593 | 0.451 | 0.572 
120 | 4.975 | 2.151 | 1.770 | 0.589 | 0.448 | 0.565 
130 | 5.014 | 2.168 | 1.793 | 0.585 | 0.445 | 0.558 
140 | 5.054 | 2.186 | 1.818 | 0.581 | 0.443 | 0.550 
150 | 5.094 | 2.205 | 1.846 | 0.579 | 0.442 | 0.542 
160 | 5.135 | 2.224 | 1.874 | 0.576 | 0.440 | 0.534 
170 | 5.177 | 2.243 | 1.904 | 0.572 | 0.437 | 0.525 
180 | 5.220 | 2.262 | 1.936 | 0.568 | 0.434 | 0.516 
190 | 5.263 | 2.281 | 1.968 | 0.564 | 0.431 | 0.508 
200 | 5.308 | 2.301 | 2.002 | 0.560 | 0.428 | 0.499 
210 | 5.354 | 2.321 | 2.037 | 0.555 | 0.425 | 0.491 


4 Bridgman, Proc. Am. Acad. Arts and Sci. 70, 94 (1935). 


% Mott and Jones, The Theory of the Properties of Metals 
and Alloys, p. 143. 
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The Debye characteristic temperature 


A calculation of the Debye characteristic 
temperature from the elastic constants requires 
the evaluation of the reciprocal cubes of each of 
three propagation velocities, averaged over a unit 
sphere in the crystal.'* The latter problem can 
only be solved by approximate methods, whose 
difficulty increases enormously when the ani- 
sotropy of the crystal is large.“ In the present 
instance the averages are calculated by simply 
evaluating the velocities at 386 points on the 
sphere and taking the arithmetic means. The 
result at 0.80°K is, =164°K. The value ob- 
tained from specific heat data is 9 = 159°K."" 


TABLE V. The tsothermal elastic moduli, constants, and 
compressibility of sodium. 


ISOTHER- 
ONSTANTS 
BILITY 
T°K su x10" 
80 4.829 2.084 0.601 0.456 1.98 
90 4.864 2.099 0.597 0.453 2.00 
100 4.901 2.115 0.592 0.449 2.01 
110 4.939 2.131 0.586 0.445 2.03 
120 4.977 2.147 0.581 0.441 2.05 
130 5.018 2.164 0.575 0.436 2.07 
140 5.058 2.181 0.571 0.432 2.09 
150 5.099 2.199 0.567 0.430 2.10 
160 5.140 2.218 0.564 0.428 2.11 
170 5.183 2.237 0.560 0.425 2.13 
180 5.226 2.256 0.556 0.422 2.14 
190 5.269 2.275 0.552 0.419 2.16 
200 5.315 2.294 0.546 0.415 2.18 
210 5.361 2.313 0.541 0.410 2.20 


In conclusion, the writers gratefully acknow]l- 
edge their indebtedness to Dr. Don Kirkham and 
Dr. Seymour Rosin for their help in the design 
and construction of the diffraction camera and 
handling chamber, and to the Ernest Kempton 
Adams Fellowship Committee of Columbia 
University for the financial assistance which 


made this research possible. 


16 Reference 15, p. 3. 
aa and Zeidler, Zeits. f. physik. Chemie 123, 403 
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Direct Determination of Crystal Structure from X-Ray Data 


MELVIN AVRAMI 
University of Chicago, Chicago, Illinois 
(Received May 23, 1938) 


The empirical similarity in the behavior of scattering- 
factor curves for different atoms is used to replace the 
structure-factor formulae by an approximating system. 
This system is shown to lead to m’th degree algebraic 
equations, designated fundamental equations, whose coeffi- 
cients are calculable in terms of crystal structure and 
atomic scattering factors. The roots of these equations 
give the coordinates of the unknown atomic positions. 
Thus if the values of a sufficient number (2m in general; 
m+1 where there is a center of symmetry) of sequential 
structure-factors are known, the problem of direct structure 
determination is solved. If only the absolute magnitudes 
are known, one may derive reciprocal equations of the 


m(m—1)th degree whose coefficients are calculable jp 
terms of the squared moduli of the structure factors. We 
can then write fundamental equations of degree m(m—1) /2 
whose roots give the coordinates of the interatomic dis. 
tances from which the unknown parameters themselves 
may readily be determined. If the number of unknown 
parameters in any projection is m, the minimum number of 
structure factor moduli required in this case is m(m—1)+1, 
The method developed in the paper is applied to the 
KH,PO, crystal, and the results are found to compare very 
well with the values found by West on the basis of a 
trial-error and Fourier synthesis investigation. 


N the problem of determining crystal structure 
from x-ray data, only the absolute values of 
the ‘“‘crystal structure factors’’ associated 
with Laue indices Akl, are known." ? The 
are related by theory to the distribution of scat- 
tering centers (atoms or ions) in the unit cell as 
follows: 


Here f,;(hk/) is the scattering factor in the Akl 
direction of the 7’th center, xyz are its coordinates 
in fractional parts of the unit axial ratios, and m 
is the number of scattering centers in the unit 
cell. The f;(Ak/) may be taken as known from 
independent experimentation. From these and 
the absolute values of the Fi,.: we seek to deter- 
mine the x;y; and 2;. In the past crystallographers 
have resorted to methods of trial and error for 
solving (1). We shall here show how these re- 
peated trials may be avoided and the structure 
obtained directiy from the experimentally given 
quantities. This method can be seen to be con- 
nected with a quasi-direct method of Patterson,’ 
but has, however, many advantages over the 
latter. A more detailed treatment will be given 
in a forthcoming paper in the Zeitschrift fiir 
Kristallographie. 


1W. L. Bragg and J. West, Zeits. f. Krist. 69, 118 (1928). 

2A. H. Compton and S. Allison, X-Rays in Theory and 
Experiment. 

3W. L. Patterson, Zeits. f. Krist. 90, 517 (1935). 


I. THE APPROXIMATING SYSTEM 


To carry through our treatment we make use 
of a simplification permitted by the experimental 
data. Examination of any table of atomic scat- 
tering factors‘ will show that the general variation 
at low temperatures of the f;(k/) with hl is simi- 
lar for different atoms.’ Thus to an approxima- 
tion (which is better the lower the temperature 
and the larger the values of hk/) we may put 
fi(hkl) =f f(hkl) where f; is a positive constant 
characteristic of atom j, and f(hk/), the same for 
all atoms, may be taken as 


/ 


Thus (1) may be written 


f(hkl) 


and it is this approximating system which we 
shall solve for the x;y,2;; these and the phases 
of the F’s are the only unknowns. It can be 
shown that the quantities in (2) may be regarded 
as the formal Fourier coefficients of a step func- 
tion whose highest order singularities (poles) 
coincide in position with the centers of scattering 
in the crystal. 

It is possible to set up equations depending 
only on the (F/f)’s whose roots are the x;y; and 


4 James and Brindley, Phil. Mag. 12, 81 (1931). 


300 


UME 54 


ulable in 
‘tors, We 
(m )/2 
omic dis. 
lemselves 
unknown 
umber of 
1) +41, 
1 to the 
are Very 
isis of a 


ike use 
mental 
iC scat- 
riation 
is simi- 
OXima- 
-rature 
ay put 
ynstant 
me for 


DETERMINATION OF CRYSTAL STRUCTURE 301 


z; respectively. These cannot be used in their 
simple forms to determine the x,y; and 2; in the 
general case because the phases of the (F f)'s 
are unknown. Nevertheless, as we shall show later, 
it is possible to eliminate the unknown phases 
and get equations entirely similar in form, so 
that it is advantageous to work out these equa- 
tions and their solutions. 

Suppose that it is the x; which we wish to 
determine. Then keeping & and / fixed and setting 


we may write from (2) 


An m’th degree equation of which the a’s are 
the roots may be determined as follows. Into 


we successively substitute the various a's and 
multiply by the associated d's. Adding and using 
we find 

Q0A p +amA pim=O, (4) 


which may be regarded as a recursion formula 
for the A’s applicable for all p. Choosing a set 


™ of m+1 successive equations of the form (4) we 
nd 
A, A A A pim 
a,=(—1)"|A A ptr A (5) 
A p+m-l A pim+r—2 A ptm+r A p+2m—1 
or (3) may also be written I]. APPLICATION TO THE CASE OF 
UNKNOWN PHASE FACTORS 
1 a 
A, GB ) cannot be used in its simple form to deter 
« mine the x;y; and zs; because the phase of the 
A pim-1 Apim A p+2m—1 (F/f)’s, and hence of the A’s, are unknown. 


While (3’) may be evaluated for different com- 
mencing indices p and different fixed values of 
kand J, it is actually independent of these except 
perhaps for a common proportionality factor. 
(3) or (3’), and its analogs for the determination 
of y; and z;, may be designated fundamental 
equations. Once these are solved, the identifica- 
tion of the roots, i.e. the correct association of 
each a; with its d;, may be carried out by solving 
the linear system obtained upon selecting any 
m equations from the system (2’). 

In cases where the fundamental equations are 
too high in degree to be solved by elementary 
methods, general procedures for the approxima- 
tion of roots must be resorted to. With the aid 
of the recursion formula (4) a valuable method 
of this kind may be established. This simple 
method, which is applicable to any system of the 
form (2'), with a's general, and therefore to the 
solution of any algebraic equation (since such 
may always be associated with a system (2’)), 
will be described elsewhere. 


In certain cases, i.e., where the distribution has 
a center of symmetry, the (F/f)’s are all real 
and the ambiguity is merely one of sign' which 
can often be determined by other methods. In 
such cases the distribution of scattering centers 
may be obtained independently of our method 
by a Fourier synthesis. There are also important 
cases where (as in KH2PO,; see below) the phases 
are known of a sequence of F's sufficient in 
number to determine the unknown parameters. 
In the most general case, however, we are faced 
with an ignorance of the phases and the necessity 
of eliminating them. To do this we multiply (2) 
by its conjugate complex. We then get 


file 


Introducing the abbreviations 


(2) 
ch we 
phases 
‘an be 
rarded 
» func- m m 
poles) 
tering 
nding 
y; and 
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| 


j=1 


Eis = Ew, 


6 
= ( ) 


m(m—1)=n, 


where w=1, 2, ---” runs over the entire range 
for j, s j#s, we get 


w=L 


which is exactly the same as (2) in form, and may 
be treated in the same way to get fundamental 
equations whose roots are the é’s, y’s and ¢’s. 
These roots occur in conjugate complexes, i.e. 
the fundamental equations are reciprocal and may 
be replaced by equations of half their degree. 
The roots give the coordinates of the interatomic 
distances rather than those of the atomic posi- 
tions. The latter may easily be deduced from the 
former when the associated pair of atoms is 
identified by the procedure outlined above for 
(2). 

Upon examination it may be seen that these 
roots correspond very closely to the positions of 
the peaks in Patterson’s ‘‘sharpened up”’ series.* 
These positions, which Patterson’s method can 
rarely locate because of the slow convergence of 
his series, are here determined in a simple alge- 
braic way. 

In case some of the positions of the scattering 
centers are known, their contributions may also 
be subtracted from | Fyx:|?/f?(hkl) and an ex- 
pression of the form (7) with m smaller than 
m(m—1) will result. 

When the experimental data are inaccurate, it 
will be advisable to get the coefficients of the 
fundamental equations by averaging the results 
from a number of the available sequences of F's 
(G's). Since both our assumptions and the ex- 
perimental data are generally more valid for F's 
with higher indices, it will be wise to restrict 
ourselves to the use of these wherever possible. 
In the case of symmetry, F,=F_,, and the 
minimum number of sequential F’s (including Fy) 
which must be known to determine an n’th 
degree fundamental equation is n+1. 


TABLE I. 
— 

hkl SIN 6 F (corr.) | fSx+fp to G 
060 0.248 40 10.8 2.3 —14 
260 | 0.262 49 10.3 | 24 37 
460 0.298 39 9.1 By 1.5 
660 | 0.351 18 7.5 14 | ~8% 
860 0.413 20 6 1.05 —3.8 
1060 0.482 24 4.4 0.75 85 
1260 | 0.554 14 2.35 | 0.6 77 


III. APPLICATION TO KH2PO, 


We shall now apply our method to the simple 
case of the KH2PQO, crystal which has been in- 
vestigated very carefully by J. West.® The space 
group is V4". The unit cell contains four 
molecules, 


4 P atoms at 000; 301; 
4 K atoms at 003; 303; 


and 16 O atoms in the general positions xyz; 
VEZ; 9x2; (2—x), ¥, (4-3); (2-9), (2+}); 
(x+2), 9, (4-2); (y+), x, and these 

323- It is the parameters xyz which we wish to 
determine. 

With West we may neglect the effect of the 
eight H atoms. Then, considering first the hk, 
reflections, it is found that “Owing to the 
special positions occupied by the potassium and 
phosphorus atoms, the (hk) spectra reduce to 
two types. . . viz. 

Type I. (A even) Contribution of atoms is potassium 
+ phosphorus + oxygens, 
[Type II.] (h odd) Contribution from oxygen atoms only, 


The spectra of type II are all weak compared 
with those of type I which receive a full con- 
tribution from potassium and phosphorus. We 
conclude that the signs of the phase factors (and 
therefore of the structure amplitudes F(hk»)) 
for type I spectra are governed by the potas- 
sium and phosphorus atoms. It follows that 

. . we shall know the signs of . . . F(hko), for 
all planes with # even though not with h odd.” 

The knowledge of the signs of a sufficient 
sequence of F’s permits us, in this case, to de- 
termine a fundamental equation for the atomic 
coordinates themselves. The contribution to 
F(hko), with h even, from K and P is 4[fi(hko) 


+fpr(hko)]. Subtracting this and dividing by 


5 J. West, Zeits. f. Krist. 74, 306 (1930), 
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fo(hko) we get a set of numbers, Gix,, which, 
from (8), can be seen to be equal to 


= 16[ cos + 
+cos (9) 


Thus, keeping & constant, and varying h over 
even values, we see that we have a system, with 
a center of symmetry, in which the roots are 
etiz, and their conjugate complexes. The 
fundamental equation for the determination of 
these roots is then a reciprocal equation of the 
fourth degree. To determine its coefficients we 
here use only one of the available sequences of 
F’s since this yields adequate accuracy. The 
values, as calculated from Fyso and f curves 
given by West, are given in Table I. These values 
yield 


1 a ae a‘ 

3.7 —1.4 3.7 15 —8.6 
-1.4 3.7 15 -86 —3.8/= 

3.7 15 -86 —3.8 8.5 

15 -86 —3.8 8.5 7.7 


2130a4 — 608a* + 3178a?+211a+2020=0. (10) 


Since this fundamental equation should be re- 
ciprocal and the symmetric coefficients equal, 
i.€., We may average these. Doing 
this, and dividing through by the first (last) 
coefficient, we have 


(10’) 


This. reciprocal equation may be solved by 
elementary methods and the roots are found to 


be 


or in degrees x=31°, y=53°, 


which agrees well with West's values of 29° and 
52°. Once x and y are known, the third parameter 
is readily found in the manner given by West. 

The remarks of West concerning the trial-error 
and double Fourier series methods apply with 
equal force to our procedure. ‘“The extent to 
which these methods are applicable without 
modification to complex structures must depend 
on the structures e.g. whether some of the param- 
eters of the heavier atoms are fixed by sym- 
metry considerations. In any case, however, it 
should not often be necessary to use a laborious 
trial and error method throughout to determine 
completely the structure. By the aid of certain 
reasonable assumptions and a qualitative exam- 
ination of the experimental data a complex struc- 
ture may often be partially solved. This partial 
solution may well be sufficient to determine the 
signs of a considerable portion of the F values 
for the reflections from planes in some important 
zone.” 

In conclusion we may make one general re- 
mark. Our method bears the same relation to 
the ordinary trial and error method as the 
systematic solution of a polynomial equation 
(e.g. the fundamental equation) bears to a cut 
and try attempt to determine its roots from linear 
combinations of homogeneous powers of these 
roots (e.g. 2’). 

The author wishes to acknowledge his in- 
debtedness for helpful advice and suggestion to 
Dr. W. H. Zachariasen and Dr. W. Bartky of 
the University of Chicago. 


3.7 

1.5 

—8.6 

—3.8 

8&5 

7.7 

simple 

een in- 

space 

s four 


AUGUST 15, 1938 


PHYSICAL REVIEW 


VOLUME 5, 


LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Evidence for the Production of Penetrating Secondary 
Cosmic-Ray Particles in the Atmosphere 


At the suggestion of Professor A. H. Compton, the fol- 
lowing experiment of the Hsiung! type was performed to 
obtain some evidence concerning the production of pene- 
trating cosmic-ray secondaries in the atmosphere. This 
subject was formerly discussed by A. H. Compton,’ and 
now the problem has renewed interest because of the dis- 
covery of the heavy electron and a recent paper by Bowen, 
Millikan and Neher? in which they suggested that all radia- 
tions observed at sea level are secondaries produced in the 
upper atmosphere. 

For the present experiment, four Geiger-Miiller tubes 
(2.5 cm diam. and 30 cm long) were arranged for fourfold 
vertical coincidence according to Fig. 1. The writers are 


NY 
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indebted to Mr. Francis Shonka for these tubes. 8.1 cm of 
lead were placed between tubes 2 and 3 to absorb the soft 
radiation, thus enabling one to record only the penetrating 
rays. The bottom tube was shielded from soft shower 
particles by 1.25 cm of lead, and a 0.625 cm lead shield was 
placed about tube number 3 to absorb Compton recoil elec- 


trons produced in the lead plate below. Any penetrating 
ionizing rays produced in the lead between the counters by 
a nonionizing incident ray will not be recorded since the 
upper two tubes will not be activated. However, if part of 
the lead (in this case 2.2 cm) is moved to a position above 
all four tubes, penetrating ionizing rays produced in the 
upper piece of lead will be recorded if they traverse the 
four tubes. Therefore, by alternately recording fourfold 
coincidences with the same piece of lead first bet ween tubes 
two and three (position B) and then above all four tubes 
(position A), one may determine the relative number of 
penetrating rays produced in that piece of lead. 

The experiment was performed in an aeroplane up to an 
altitude of 25,000 feet. The results are shown in Table |, 


TABLE I. Results of measurements. 


13,000 ft. to 20,000 ft. to 
ALTITUDE 20,000 it. 25,000 ft. 25,000 ft. 
Position of 
a A B A B A B 
See (Fig. 1) 

Number of 7 5 s& 4 11 2 
Fourfold 4 5 9 7 10 4 
Coincidences 6 4 9 § 10 11 
in 3 min. 6 7 11 5 10 2 
intervals. 1 ) 10 6 

15 6 
Average 4.8+0.7/5.2+0.7) 9.2+1.0 |6.0+0.7|11.0+0.9 |5.2406 
Ratio A/B=0.9| Ratio A/B=1.5 | Ratio A/B=2.1 


Three-minute readings were made with the lead alternately 
in positions A and B. Up to 20,000 feet no appreciable 
differences between A and B were observed. From 20,000 
feet to 25,000 feet, the excess with the lead above is quite 
apparent, and at 25,000 feet the ratio A/B=2.1. This is 
well outside of the rather large statistical errors. 

This experiment indicates that above 20,000 feet pene- 
trating ionizing rays (heavy electrons) are produced by 
nonionizing agents (photons). It is of interest to note that 
the ‘‘hump” in the cosmic-ray absorption curve at about 
25,000 feet is made more prominent when the apparatus is 
shielded with 6 cm of lead.‘ This indicates that this hump 
is possibly related to the production of penetrating secon- 
daries. 

The results of the present experiment are in good agree- 
ment with Heitler’s® assumptions that the heavy electrons 
are formed by secondary photons which are abundant at 
high altitudes and are strongly absorbed by 2 cm of lead. 
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With the lead in position B the fourfold vertical coinci- 
dences at 25,000 feet were three times as numerous as at 
sea level. This represents the ratio of the intensities of the 
vertical penetrating component, which agrees well with 
Heitler’s® calculations. Braddick and Gilbert’ found a 
ratio of 9 to 1 at 35,000 feet. This high ratio is explained by 
Heitler as due to insufficient shielding of the soft rays 
(3 cm of lead). 

Twofold coincidences for the upper two tubes at 25,000 
feet and at sea level give a ratio of 12 to 1 for the total 
vertical radiation. This is in good agreement with Pfotzer’s® 
data obtained by balloon flights with threefold coincidences. 
Special care was taken to have the efficiency of the ap- 
paratus a maximum. By activating the tubes with a 
radium source, conditions similar to those at 25,000 feet 
were reproduced in the laboratory, without any apparent 
change in the efficiency. The above agreements seem to 
indicate this also. 

The writers wish to express their appreciation to 
Professor A. H. Compton not only for suggesting the 
exp2riment, but also for his continued support. 

MARCEL SCHEIN 
Votney C. WILSON 
Ryerson Physical Laboratory, 
University of Chicago, 


Chicago, Illinois, 
August 2, 1938, 


1D.S. Hsiung, Phys. Rev. 46, 653 (1934). B. Rossi, Proc. Lond. Conf. 
on Nuclear Phys. (1934). 
2A.H. Compton, Proc. Phys. Soc. London 47, 747 (1935). 
3T. S. Bowen, R. A. Millikan, and H. V. Neher, Phys. Rev. 53, 217 
1938). 
' 4A, H. Compton and R. J. Stephenson, Phys. Rev. 45, 442 (1934). 
§W. Heitler, Proc. Roy. Soc. A166, 529 (1938). 
6W. Heitler, Proc. Roy. Soc. A161, 261 (1937). 
7T. Braddick and C. W. Gilbert, Proc. Roy. Soc. A156, 570 (1936). 
8G. Pfotzer, Zeits. f. Physik 102, 23 (1936). 


A New Rydberg Series in N, 


The absorption spectrum of nitrogen below 1000A shows 
some interesting features that appear to have been over- 
looked by previous investigators because of insufficient 
dispersion. The lower spectrum of Fig. 1 (negative; ab- 
sorption appears white) shows the general aspect between 
785 and 850A at a pressure of 0.1—-0.2 mm and path length 
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130 cm. The complex of overlapping band systems at longer 
wave-lengths definitely terminates at about 800A with a 
converging series of band heads. These are shown with 
greater enlargement in the upper spectrogram. They ob- 
viously constitute the higher and intermediate members of 
an electronic series, the lower members of which (1 <6) are 
presumably confused with the complex of bands at longer 
wave-lengths. Diffuse absorption (O2?) overlaps the bands 
n=6, 8, 10, 11, 12, 13, blotting out the first three and also 
the fourth except on two plates taken at lower pressures. 

The frequencies of the heads follow a Rydberg law closely. 
With n=7 to 25 and a=0.34 (n+a=Rydberg denomina- 
tor) only the lowest measured member deviates appreciably, 
being shifted about —10 cm™'. The limit of the series lies 
at 125,670+2 cm“, or 15.503 volts. This is just the first 
ionization potential of No computed indirectly by Mulliken! 
from the data on Hopfield’s Rydberg series. The latter 
converge to the _ ah excited state of N»*, while the limit 
of our series is the 2~,* normal state. The latter has an 
internuclear distance of 1.11A as compared to 1.09 for the 
normal state of neutral No», a circumstance which is favor- 
able to the observation of such series. Further evidence 
that this is actually an electronic series comes from the 
fact that no discrete bands are observed beyond the limit, 
but a region of continuous absorption extends from here 
with gradually diminishing intensity to at least 600A. The 
Hopfield series just referred to does not appear on these 
plates, while at higher pressures the continuum below 
800A is almost completely absorbed. 

That the high dispersion and resolution of a 3-meter 
grazing-incidence spectrograph were required to show the 
existence of this series will be apparent from the spectro- 
grams. We are using the Lyman continuum, which covers 
the longer wave-length region not accessible with the 
helium continuum. A thorough analysis of the whole ab- 
sorption spectrum below 1000A is in progress. 

R. E. Worry 
I. A, JENKINS 
Department of Physics, 


University of California, 
July 19, 1938. 


1R.S. Mulliken, Phys. Rev. 46, 144 (1934), 
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Remarks on the Article by H. R. Crane and J. Halpern— 
“New Experimental Evidence for the Neutrino” 


The beautiful 8-ray tracks of Cl’, obtained by H. R. 
Crane and J. Halpern, show at their initial points clusters 
of droplets which the authors attribute to the ionization 
produced by the recoiling atom of A**, From the statistical 
relation between the number of ions and the energy of the 
B-ray they calculate the value of the kinetic energy of the 
atom which they find to be in general agreement with the 
hypothesis of the simultaneous emission of an electron and 
a neutrino. 

I should like to point out that the relation between the 
energy of an atom and the ionization produced by this atom 
is much too complicated to allow any direct test of this 
kind. The average energy expenditure per ion pair increases 
rapidly when the energy of the ionizing charged or un- 
charged atom tends towards zero. Even in the case of the 
a-recoil, when the energy of the atom is more than 1000 
times larger than the energy of the recoiling A** atom, the 
average energy spent for a pair of ions is of the order of 
300 ev. Another example is the ionization produced in argon 
by moving argon ions, which starts at 240 ev and even at 
400 ev is many times smaller than the maximum specific 
ionization, corresponding to about 30 ev per ion pair.! 

It seems very probable that the A atoms, having an 
energy of about 400 ev, are either uncharged or undergo 
frequent losses of charge. The 8-recoiling atoms of naturally 
radioactive elements are, as is well known, uncharged. 

On the other hand, even if the energy of the A atoms were 
converted into ionization with the maximum efficiency, 
one must not consider the absolute value of the kinetic 
energy, but the energy of the movement relative to the 
center of gravity of the system formed by the argon and, 
say, a nitrogen atom. This last energy is, however, equal 
only to 14/(14+38) =0.27 of the total kinetic energy. 

All these considerations show that the number of ion 
pairs produced by a 400 ev A atom must be very small, 
probably of the order of unity. I think that the clusters ob- 
served by Crane and Halpern cannot be considered as giving 
support to the hypothesis of the neutrino, unless the energy 
of the neutrino, and, consequently, the total energy evolved 
in the 8-transformation is much larger than the energy 
deduced from the upper limit of the 8-ray spectrum. 

L. WERTENSTEIN 

Laboratory of Atomic Physics, 

Free University of Poland, Warsaw, 
July 10, 1938. 

a and Scheele, Handbuch der Physik, second edition, Vol. 22, 

p. 


A Reply to the Above Note by L. Wertenstein 


Professor Wertenstein has very kindly given us the 
opportunity of commenting upon his remarks published 
above. Professor Wertenstein is entirely justified in calling 
attention to the danger in our assumption! that one ion 
pair is formed for every 30 ev lost by the slowly moving 
argon atom in air. This value was adopted by us only as a 
working hypothesis, because of the lack of experimental 
data on the total ionization produced by argon atoms (or 
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ions) of 400 ev and less, in air. We are now undertaking a 
measurement of this quantity, and prefer to postpone an 
opinion until the results have been obtained. In the mean- 
time, we wish to point out (1) that our conclusion that a 
neutrino participates in the beta-ray disintegration can be 
arrived at without detailed consideration of the compli- 
cated mechanism by which droplets are produced in the 
cloud chamber by the recoiling atom, and (2) that a small 
efficiency of ionization by the recoiling atom does not 
necessarily imply a small efficiency in production of 
droplets. 

(1) Let us consider the comparison between the results 
obtained with radiophosphorus (upper limit ~1.7 Mev) 
and radiochlorine (upper limit ~5 Mev). More specifically, 
let us compare the recoils resulting from the emission of 
beta-rays of the same energy in the two cases, namely be- 
tween | and 1.5 Mev. The energy of recoil given to the atom 
by the beta-ray alone, is in the two cases essentially the 
same, about 25 ev, which is too little to produce a cluster of 
droplets on any reasonable assumption about ionization 
efficiency. The quantity which we may expect to differ 
radically in the two cases is the energy of the neutrino, 
which is 3.5 to 4 Mev in the chlorine case and 0.2 to 0.7 
Mev in the phosphorus case. The experiment shows that 
the chlorine disintegrations of this group produce clusters 
of from 15 to 34 droplets, while the corresponding phos- 
phorus disintegrations produce no visible clusters at all. 
Entirely similar results have been obtained by comparing 
radiochlorine with radionitrogen (upper limit 1.25 Mev). 
The most reasonable explanation of this effect seems at 
present to be that the much greater energy of recoil in the 
chlorine case is due to the escape of a neutrino of high 
energy (3.5 to 4 Mev). 

(2) We have recently become aware of a mechanism 
other than ionization, by which a moving atom may bring 
about the formation of droplets in the cloud chamber, 
namely dissociation. Whether or not this is an important 
effect in our experiment we are not vet prepared to say. It 
is well known that a moving atom of several hundred ev 
energy strongly dissociates the molecules of the gas through 
which it passes, for example, in air, O.—20(~6 ev) and 
N.—2N(~9 ev). The most probable fate of the free atoms 
is combination (in triple collisions) with molecules to form 
NO:, N2O, and other combinations, including those due to 
reactions with the alcohol vapor, which is present in the 
cloud chamber. We have attempted a preliminary test of 
the ability of such molecules to act as centers of condensa- 
tion under normal operating conditions in a cloud chamber. 
The injection of minute amounts of NO» gas was found to 
produce an extremely dense cloud. Injection of air by the 
same method produced no cloud. It is therefore possible 
(but must be proved) that every molecule of NO» and of some 
other compounds resulting from dissociation by the recoil 
atom is effective in bringing about the condensation of a 
droplet.? If this effect should prove to be at all large, it 
would mean that the interpretation of droplet counts in 
terms of energy of recoil would have to be made on the 
basis of a direct experiment, and not on the basis of the effi- 
ciency of ionization alone. It will be fortunate if such an 
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effect is large, because it will increase the sensitivity for 
detection of slow heavy particles in a cloud chamber. 

To sum up, (1) the conclusions of the experiment are ar- 
rived at by a comparison of Cl*8, P® and N"™, rather than 
by an absolute estimate of the energy of recoil, and (2) at 
least one possible mechanism in addition to ionization exists 
for the production of droplets by slowly moving atoms, but: 
its importance relative to ionization remains to be de- 
termined. 

H. R. CRAneE 


J. HALPERN 
University of Michigan, 
Ann Arbor, Michigan, 
July 30, 1938. 


1 Crane and Halpern, Phys. Rev. 53, 789 (1938). 

2 The gases in the chamber were air, alcohol and ethylene dichloride 
in the experiment with Cl; nitrogen, alcohol and PHs; in the P® ex- 
periment; air and alcohol in the N experiment. 


The Hard and Soft Component of Cosmic Radiation in 
Terms of a Common Primary Background ' 


The fundamental assumptions which control the view 
here discussed are the following: 

1. There is a single primary background of charged, 
but nonelectronic type, with energy distribution function 
F(E) of the form A /E?. 

2. Each primary ray experiences, per unit of path, a 
number of ‘‘events’’ or collisions, the average number of 
these per unit of path being 8. This number we assume 
independent of the energy, and we assume the events to be 
distributed individually and collectively at random in any 
interval concerned. It results from the latter part of the 
assumption that the probability of occurrence of » events 
within the interval 8x is (8x)"e**/n! 

3. In each of these events the ray loses an amount of 
energy ¢, which is the same for all energies. No ray of 
energy E can, of course, experience more than m losses e, 
where me= E., 

4. If a ray has energy less than e¢, it has, nevertheless, 
the same chance of experiencing an event; but in such an 
event it loses all of its energy. 

5. When a primary ray loses energy, the energy loss 
goes into secondaries (possibly in part heavy electrons) and 
some of it gets converted directly or indirectly into ordinary 
electrons, which then proceed according to the Bethe- 
Heitler theory, the energy being lost, finally, in ionization 
and in radiative losses in all cases. 

6. The story of the relationship between the number of 
secondaries per primary and the loss of energy of the 
primaries is controlled by 


ndn/dx+an=dW/dx, (1) 


where 7 is the average energy of all the particles other than 
primaries, a is the loss per particle per unit of path in 
ionization and in radiative losses, ” is the number of secon- 
daries and higher order particles per primary, and W is the 
average energy lost by a primary between the top of the 
homogeneous atmosphere and the depth x. 

On the basis of the foregoing assumptions, expressions 
have been developed for the intensities of the primaries and 
secondaries—using the latter term to include all the higher 
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order particles. The particular case considered is one where 
8 is less than the energy for entry through the earth's 
magnetic field at the magnetic equator, but is equal to the 


energy for entry at some higher magnetic latitude ¢, taken 
as 42° in the example. 


2 2 d 
Metecs of walter 


Intensity 


Fic. 1. 


For both the primaries and the secondaries, the intensity 
comes out in the form of a series whose terms correspond to 
the probabilities of successively larger and larger numbers 
of events. The earlier terms dominate the story at small 
values of x, and these terms show a latitude variation. At 
greater depths, the later terms of the series dominate the 
situation. These terms show no latitude effect. Moreover, 
they correspond to an ever increasing ‘‘hardening"’ of the 
radiation until finally an absorption coefficient is attained 
which depends entirely upon the distribution function 
F(E), and is independent of all other elements of the theory. 
This absorption coefficient varies with x and is given by 
u=(p—1)/x. 

We thus realize a condition in which the /atitude effect 
increases with altitude, and in which the measured radiation 
at any latitude ‘‘hardens”’ with increase of depth. 

The growth of secondaries as determined by (1) results in 
a maximum for the intensity of the secondaries. Moreover, 
this maximum occurs at depths which are greater the smaller 
the latitude. 

Figure 1 shows, for magnetic latitude 42° and for the 
equator, the intensity of the secondaries plotted against 
depth for 8=4 per meter of water, «=6X 10° ev, and p=3. 
These curves are only drawn for illustration and not 
necessarily with the best choice of constants to fit all the 
experimental data. 

The large value of Be determines a large number of 
secondaries per primary, which results in the phenomena 
measured being mainly secondary. 

W. F. G. Swann 

Bartol Research Foundation of the Franklin Institute, 


Swarthmore, Pennsylvania, 
July 25, 1938. ° 


1 Part of a paper presented at the Cosmic-Ray Symposium held at 
the University of Chicago, June 27-30, 1938, which paper is in process 
of publication in full in the Journal of the Franklin Institute. 
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Search for Gamma-Rays from the Deuteron-Deuteron 
Reaction 


The discovery by Bonner! of a second group of neutrons 
from the reaction 


H?+ H*+He'+n! (1) 


has suggested the possibility that He*® may be formed in an 
excited state about 2 Mev above the ground state. The 
difficulty of reconciling such a state with present nuclear 
theory has been emphasized by Share? and Schiff.* If such 
a state exists we should expect to observe gamma-rays due 
to the transition of He? to its normal state. Evidence for a 
gamma-ray has been given by Kallmann and Kuhn,‘ who 
used coincidence counters to investigate the reaction. We 
have made a search for this gamma-ray by means of a 
cloud chamber. 

The cloud chamber was placed about 25 cm from a heavy 
phosphoric acid target bombarded by 0.5 Mev deuterons. A 
lead collimator was placed between the target and the 
cloud chamber and both were provided with thin aluminum 
windows. A sheet of carbon 0.15 cm thick extended across 
the center of the chamber as a source of Compton electrons. 
A magnetic field of 715 gauss was used to determine the 
electron energies. Many recoil electrons were observed, but 
most of these appeared to be due to the gamma-rays from 
slow neutron capture and from the inelastic scattering of 
fast neutrons in the material around the chamber. In order 
to determine whether any gamma-rays came directly from 
the target, the following experiment was performed. On 
alternate expansions the hole in the lead collimator was 
filled with a lead block 4 cm thick. This would have de- 
creased the intensity of a gamma-ray from the target by a 
factor ten. Actually no difference in the numbers of elec- 
trons between 1 and 2 Mev energy was observed in the 
two cases. 

In order to determine the ratio of the number of gamma- 
rays to the number of neutrons from the reaction, the 
target holder was placed directly within the cloud chamber 
and surrounded by a foil of Cellophane 0.1 mm thick. From 
the number of proton tracks coming from this foil the 
neutron intensity could be estimated using the known cross 
section for neutron-proton scattering. The gamma-ray in- 
tensity was determined simultaneously by counting the 
Compton electrons produced in a sheet of carbon in the 
center of the chamber. Since the time of sensitivity of the 
cloud chamber for electron tracks is smaller than the time of 
sensitivity for proton tracks, a correction factor was 
necessary. An estimate of this ratio was found by placing in 
the chamber a natural radioactive source which emitted a 
known ratio of beta-rays and alpha-particles. From the ratio 
of Compton electrons to protons obtained in this way it was 
calculated that there is not more than one gamma-ray for 
every 200 neutrons. This figure agrees with the estimate of 
Kallmann and Kuhn;‘ the intensity of Bonner’s second 
neutron group! is, however, one-tenth of the 2.6 Mev group. 

During the course of this investigation many protons 
were observed which penetrated a carbon sheet 0.15 cm 
thick, and consequently had an energy greater than 15 
Mev. Since these were present when the chamber was 
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separated from the target by 4 cm of lead, they must be 
due to neutrons, possibly from the secondary reaction 


1'+17.6 Mev (2) 
due to recoiling H* nuclei from the reaction 
H?+ H*+H3+H!. (3) 


We understand that energetic protons from the analogous 
reaction of He* with H? have been reported by Oliphant.s 
The ratio of the number of these very energetic recoil pro- 
tons to those of the 2.6 Mev group was of the order of one 
to one thousand; consequently reaction (2) must be an 
exceedingly probable one. 

I wish to express my thanks to Dr. H. R. Crane for his 
helpful advice and encouragement. 


ARTHUR J. RUHLIG 


Department of Physics, 
University of Michigan, 
Ann Arbor, Michigan, 
August 1, 1938. 


! Bonner, Phys. Rev. 53, 711 (1938); also Baldinger, Huber and Staub, 
Helv. Phys. Acta 11, 245 (1938). 

2? Share, Phys. Rev. 53, 875 (1938). 

3 Schiff, Phys. Rev. 54, 92 (1938). 

*Kallmann and Kuhn, Naturwiss. 26, 106 (1938). 

* Unpublished; private communication from Professor H. A. Bethe, 


Observation of H' and H® Ranges from the Disintegration 
of Deuterium by Deuterons 


Recent cloud chamber experiments! on the energy 
spectrum of recoil protons produced by neutrons from the 
reaction 

have shown the presence of a low energy group of neutrons, 
and indicate the existence of an excitation level in He* at 
1.89 Mev. Ionization chamber studies? have led to this 
same conclusion. In view of this excitation level in He® 
it seems reasonable to expect a similar level in H*. Conse- 
quently we have looked for a short range group of protons 
from the reaction 

We used apparatus previously described. The protons 
were observed in a cloud chamber filled with air and 
alcohol vapor. They were allowed to pass into the chamber 
through an aluminum foil which had a stopping power of 
approximately 0.5 cm. A slit system defined the direction 
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Fic. 1. Range distribution of protons from the disintegration of 
deuterium by deuterons. 
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of the disintegration particles to an angle of 99.5°+2° to 
the direction of the incident 100 kev deuterons, which 
bombarded a D;PO, target. The stopping power of the 
gas in the chamber (1.18) was determined from the mean 
range (14.7 cm) of the proton group of maximum energy. 

Figure 1 gives the range distribution of the protons 
observed, and, contrary to what one might expect, there is 
no indication of more than a single group. Preliminary 
experiments with a D,O target frozen by liquid air gave 
similar results, though a higher background was present. 
From this experiment and that described below we con- 
clude that H® is excited to a level between 0.4 and 2.9 Mev 
in less than 2 percent of the disintegrations, if at all. 

While looking for low energy protons from the deuterium 
reaction we observed the range of the H? particles produced 
in the same reaction. An accurate determination of the 
range of these particles is useful in checking the range- 
energy relation for H® particles, and hence the range-energy 
curve for short range protons. 

The H®* particles were observed in a cloud chamber 
filled with helium and water vapor. They were allowed to 
pass into the cloud chamber through the slit system and an 
aluminum foil with a stopping power of 0.50 cm.‘ Fig. 2 
gives the range distribution of the particles. The stopping 
power of the gas was found to be 0.30+0.01 cm for particles 
with a range of 3.80 cm (polonium alpha-particles). 

The most probable range is 1.28+0.10 cm, which corre- 
sponds to a mean range of 1.31+0.10 cm at 99.5°. Since the 
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Fic. 2. Range distribution of the charged particles as observed in the 
helium filled chamber. 


stopping power of the helium would probably be slightly 
higher for 1.3 cm particles than for polonium alpha- 
particles, this is actually a lower limit to the H# particle 
range. 

The energy of the H?® particles emitted at 99.5° was cal- 
culated from the Q value of 3.98 Mev.® The energy of the 
H? particles is 0.89 Mev at this angle. The range of a H® 
particle of this energy theoretically is three times the range 
of a proton with an energy of 0.89/3 Mev. If we use the 
new range-energy curve of Livingston and Bethe for 
protons, the mean range of the H® particles should be 
0.95 cm. This is beyond the limits of experimental error 
and we conclude that their range-energy curve for protons 
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in this region is in error. The experimental range of 1.31 cm 
is in much better agreement with the old curve of Livings- 
ton and Bethe® which predicts an H* range of 1.38 cm. 
This suggested change in the range-energy relation for 
protons is important since it would necessitate raising the 
mass of the neutron, and lowering the mass of He’ and C™, 
Emmett Hupsreru 


T. W. Bonner 
The Rice Institute, 
Houston, Texas, 
June 21, 1938. 


1T. W. Bonner, Phys. Rev. 52, 685 (1937); 53, 711 (1938). 

? Baldinger, Huber, and Staub, Helv. Phys. Acta 11, 245 (1938). 

3T. W. Bonner, Phys. Rev. 52, 685 (1937). 

4The stopping power of the foit was calculated by assuming 1.41 
mg/cm? = 1 cm air. 
on em Kempton, and Rutherford, Proc. Roy. Soc. A149, 406 
* Livingston and Bethe, Rev. Mod. Phys. 9, 245 (1937). 


Ferromagnetic Colloid for Studying Magnetic Structures 


The magnetic structure of a ferromagnetic material can 
be revealed by depositing a fine ferromagnetic powder 
from suspension onto the smooth surface of a specimen 
which is suitably magnetized. F. Bitter,’ who first studied 
magnetic structures by this method, used finely ground 
gamma-ferric oxide suspended in ethyl acetate. Stray 
fields near the surface of the specimen guide such a ferro- 
magnetic powder as it settles; the powder thereby forms a 
pattern recording regions where H? is greater than in 
near-by regions. A similar technique for obtaining patterns 
has been used by most of the other investigators? who 
have extended the work of Bitter. McKeehan and Elmore,’ 
however, used a true colloid (rather than a temporary 
suspension) of gamma-ferric oxide. They discovered that 
the colloid formed patterns consisting of swarms of sol 
particles which under favorable conditions did not settle 
out or coagulate but which retained their mobility within 
the solution. Although this latter technique possesses 
many advantages over the former, no great use seems to 
have been made of it, a fact undoubtedly due to the diffi- 
culty of securing or of preparing a suitable colloidal 
suspension. This situation has prompted the writer to look 
for a simple means of making a ferromagnetic colloid for 
magnetic powder investigations. 

A series of experiments have led to the preparation of a 
very satisfactory colloidal magnetite. The new colloid 
forms patterns rapidly, does not coagulate once the pattern 
is formed, and is superior for subjective observation and 
for photography since its color is a very dark reddish-black 
(gamma-ferric oxide is a light brick-red). Several colloidal 
magnetites which were prepared proved unsatisfactory. 
One of these, protected by dextrin, had such a fine particle 
size that it formed patterns only in regions possessing very 
intense stray fields. Another one, unprotected, produced 
instead of the expected pattern tree-like structures which 
did not disperse when the applied magnetic field was 
removed from the specimen. To be useful for forming 
patterns evidently the colloid particles must not be too 
small, i.e., they must lie in the upper range of colloidal 
size. They must also be protected by another colloid to 
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prevent their coagulation when they congregate near the 
surface of the specimen. The surface itself should be 
protected by a thin layer of lacquer made, for instance, 
by dissolving a little celluloid in amylacetate. 

To make the new colloid, first prepare a coarse, flaky 
precipitate of magnetite by Lefort’s method:* dissolve 2 
grams of FeCl,-4H,O and 5.4 grams of FeCl;-6H,O (or 
equivalent amounts of the sulphates) in 300 cc of hot 
water and add with constant stirring 5 grams of NaOH 
dissolved in 50 cc of water. Filter to remove salt and 
excess sodium hydroxide. Rinse precipitate in filter several 
times with water and finally once with 0.01 N HCl. Then 
transfer precipitate to one liter of one-half percent soap 
solution and boil for a short time. The former precipitate 
will now have become entirely colloidal with the exception 
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of a very small quantity of undispersed oxide which should 
be removed by filtering hot. It is interesting that the 
success of the method depends entirely on the peptization 
of the precipitate with HCl before it is added to the 
water containing soap which serves as a protective colloid, 
A drop or two of the colloidal magnetite placed on the 
magnetized specimen will give an ample supply of sol 
particles for forming a pattern. 
W. C. ELMore 


Massachusetts Institute of Technology, 
Cambridge, Massachusetts, 
July 19, 1938. 


1F. Bitter, Phys. Rev. 38, 1903 (1931); 41, 507 (1932). 
2F, Bitter, Jntroduction to Ferromagnetism (McGraw-Hill Book Co.,, 


1937), p. 59. 
3L. W. McKeehan and W. C. Elmore, Phys. Rev. 46, 226 (1934). 


4See, for instance, L. A. Welo and O. Baudisch, Phil. Mag. 3, 396 
(1927). 
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Proceedings of the American Physical Society 


MINUTES OF THE TORONTO, CANADA, MEETING, JUNE 24-25, 1938 


HE 222nd regular meeting of the American 
Physical Society was held at Toronto, 
Canada on Friday and Saturday, June 24 and 25, 
1938 at the University of Toronto as a joint 
meeting with the American Association of 
Physics Teachers. All scientific sessions were 
held in the McLennan Laboratory. The pre- 
siding officers were Dr. Lyman J. Briggs, Presi- 
dent of the Society and Professor John T. Tate, 
Vice President. The visiting members of the 
Society were housed in the Student Residences. 
The attendance at the meeting was about two 
hundred and fifty. 

A joint meeting with the American Associa- 
tion of Physics Teachers was held on Friday 
morning at 9:30 o’clock in the McLennan 
Laboratory with Professor F. K. Richtmyer, 
President of the American Association of Physics 
Teachers, presiding. A short address of welcome 
from the President of the University of Toronto, 
the Hon. Dr. H. J. Cody, opened the meeting. 
The invited speakers were as follows: Professor 
John Satterly, University of Toronto on A Com- 
parison of the Teaching of Physics in England 
and America; Professor E. F. Burton and Mem- 
bers of the Staff in Physics, University of 
Toronto, and Teachers of Physics in Secondary 
Schools of the Toronto Area, prepared a Demon- 
stration of Experiments; and liquid helium was 
produced and exhibited at the Cryogenic Labora- 
tory of the University of Toronto. 

A symposium of invited papers was arranged 
for Saturday morning at 9:30 o'clock in the 
McLennan Laboratory with President Lyman 
J. Briggs presiding. The speakers were: Professor 
H. A. McTaggart, University of Toronto on 
Monomolecular Films; Professor J. W. Beams, 
University of Virginia on High Speed Centrifuges 
and Some of Their Applications; and Professor 
H. Grayson Smith, University of Toronto, on 
Superconductivity. Professor H. A. Kramers of 
Leiden, Holland, who was scheduled to speak on 
Magnetic Phenomena at Low Temperatures did not 


arrive from Europe in time to make his address. 
- The local committee arranged for various 
excursions as well as for golf and swimming. 
The President of the University of Toronto 
invited all members of the two societies and their 
wives and friends to a garden party on the lawn 
of the Royal Ontario Museum at 4:30 o'clock on 
Friday afternoon, June 24. 

The dinner was held on Friday evening at 7:00 
o'clock for the members of the two societies and 
their friends at Hart House on the university 
campus. President Briggs presided and called 
upon President Cody, Vice President Tate and 
Dean F. K. Richtmyer. There were about one 
hundred and fifty guests at the dinner. 

Meeting of the Council: At its meeting held on 
Friday, June 24, 1938 the death of one member 
(Elon H. Hooker) was reported. One candidate 
was elected to membership, thirteen candidates 
were transferred from membership to fellowship 
and eighteen candidates were elected to member- 
ship. Elected to fellowship: Eugene Guth. Trans- 
ferred from membership to fellowship: Daniel S. 
Elliott, Walter B. Ellwood, Walter Gordy, 
Clarence W. Hickman, Victor E. Legg, Frederick 
B. Llewellyn, Carol G. Montgomery, Rose C. L. 
Mooney, Foster C. Nix, Rogers D. Rusk, John 
C. Schelleng, Ralph D. Wyckoff, and Jerrold R. 
Zacharias. Elected to membership: William J. 
Bailey, Jr., B. Vivian Bowden, E. Bretscher, 
George P. Brousseau, Eugene Brunner, Herman 
Feshbach, Robert C. Herman, Otto G. Koppius, 
J. K. Kostko, Karl D. Larsen, H. R. Letner, 
William W. McCormick, Edson R. Peck, Roy W. 
Prince, Jr., Solomon Schwebel, Arthur K. 
Solomon, Paul F. Spremulli and Turrell Uleman. 

The regular scientific program of the Physical 
Society consisted of twenty-two contributed 
papers of which numbers 1, 7, 9, 16 and 17 were 
read by title. The abstracts of these papers are 
given in the following pages. An Author Index 
will be found at the end. 

W. L. SeveriInGnAus, Secretary 
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ABSTRACTS TO CONTRIBUTED PAPERS 


1. Spectrochemical Analysis of Fertilizers, Particularly 
for the Micrometabolic Elements. STANLEY S. BALLARD, 
University of Hawaii.*—Samples of twelve basic commer- 
cial fertilizers have been analyzed by the spectrochemical 
methods developed in the spectroscopic laboratory of the 
Experiment Station of the Hawaiian Sugar Planters’ Asso- 
ciation.’ Analyses were performed semi-quantitatively, the 
amount of each element present in a sample being indicated 
by one of the following designations: Lots, Some, Less, 
Trace, Not Detected. Such an analysis reveals the presence 
or absence in a sample of the fifty-odd arc sensitive ele- 
ments, of which twenty-one were encountered in the present 
investigation. The various fertilizers were found to differ 
greatly in the number of mineral elements present. Par- 
ticular attention was given to the minor plant nutrients, 
the micrometabolic or so-called ‘‘less essential’ elements 
iron, manganese, boron, copper, zinc and aluminum, as well 
as to the major elements phosphorus, potassium, calcium 
and magnesium which often occurred as impurities when 
they were not present as major constituents. The fertilizers 
tested included two sulfates of ammonia, two nitrates of 
soda, two nitrates of potash, muriate of potash, two raw 
rock phosphates, two superphosphates, and phosphate of 
ammonia. The last five named proved to be the best sources 
of the minor elements. This paper will appear in full in 
The Hawaiian Planters’ Record. 


* To be read by title. 
1S. S. Ballard, Phys. Rev. 52, 253A (1937). 


2. The Effects of Pressure and Temperature on the 
Viscosity of Monosubstituted Toluenes. H. E. MorGANn 
AND R. B. Dow, Pennsylvania State College.—The viscosity 
of ortho-, meta-, and para- formsof chloro-, bromo-, nitro-and 
iodotoluene was determined at 30°, 75°, and 99°C at pres- 
sures ranging from atmospheric to 4000 kg/cm? by means 
of the rolling-ball type of high pressure viscometer. Calibra- 
tion of the viscometer enabled the absolute viscosity to be 
computed at the three temperatures of the experimental 
range of pressure. The applicability and validity of the 
log coeff. of viscosit y-pressure relation to represent the data 
are discussed. The tendency of the pressure coefficient 
both to increase and decrease as the substituted group is 
moved in the series, leads to the conclusion that no general 
rule can be stated for the effect of change of position of the 
substituted group. A plot of the pressure coefficient against 
atomic weight (M.W. for NO2) shows that while the former 
tends to increase with the latter for the ortho- and meta- 
bindings, it clearly decreases for the para- binding. Increase 
of temperature reduces these differences. It is suggested 
that this anomaly may be due to a fundamental difference 
in the vibration of the toluene ring when the substituted 
group is attached in the para position, a result which 
would be expected to influence the type of “interlocking” at 
high pressures as pictured by Bridgman. 


3. Viscosity Characteristics of Lubricating Oils at High 
Pressure at Various Temperatures. R. B. Dow, The Penn- 
sylvania State College.—The study of the physical properties 


of oils at high pressures begun at Harvard University and 
being continued at the School of Chemistry and Physics, 
The Pennsylvania State College, is reviewed briefly. The 
pressure coefficient of viscosity at 100°, 130°, and 210°F 
and its relation to composition, molecular weight, and vis- 
cosity index is discussed for data obtained in this laboratory 
on Rennsylvania, Oklahoma, and California blends, and 
solvent extracted oils. Similar analysis is given for inter- 
polated data from other laboratories for lard and castor oil, 
Comparison of a rolling ball viscometer with an absolute, 
concentric cylinder type shows that both instruments are 
capable of good accuracy and satisfactory precision at 
various pressures and temperatures. No evidence of a semi- 
colloidal or non-Newtonian behavior, as observed by the 
independence of the coefficient of viscosity and the rate of 
shear, has been observed for these oils at the pressures and 
temperatures investigated. 


4. New Survey of the Spectrum of the International 
Iron Arc in the Range 10,000 to 2000A. Dorotny W. 
WEEKs, Wilson College, AND GEORGE R. HARRISON, Massa- 
chusetts Institute of Technology —Wave-length measure- 
ments have been made on some 10,000 lines obtained in the 
international iron arc. A 10-meter concave grating having 
30,000 lines per inch was used in Paschen mounting at dis- 
persion 0.4A/mm for the range 2000 to 6900A. For 6900 
to 10,000A, a 10-meter grating having 15,000 lines per inch 
was used ina Wadsworth mounting at dispersion 3.3A/mm. 
Two hundred twenty inch plates were measured on a com- 
puting-recording comparator, the resulting wave-lengths 
being self-consistent to 0.003A or better as determined by 
repeat runs and by the combination principle. After elimi- 
nation of all impurity lines and possible ghosts the number 
of known iron lines is approximately trebled. The able as- 
sistance of W.P.A. clerical help is gratefully acknowledged. 


5. Electric Arc Reignition at High Pressures. J. D. 
Cosine, R. B. Power* anp L. P. Winsor, Harvard 
Graduate School of Engineering.—A study has been made of 
the reignition of short a.c. arcs in air at pressures from one 
to eight atmospheres, for currents under 3.5 amperes, ina 
resistance circuit. Pure graphite electrodes were used and 
were separated by distances from 1/8 mm to 3 mm. As 
previously indicated for nitrogen, the reignition voltage 
under these conditions increases with pressure when the 
r.m.s. arc current is kept constant. The reignition potential 
for constant gas pressure decreases as the current is in- 
creased. The characteristic of reignition potential vs. cur- 
rent has been found to exhibit a discontinuity at a critical 
current. For this current the reignition potential decreases 
by about 1000 volts to a new characteristic typical of the 
higher currents. This critical current depends upon the 
gas pressure and to a lesser degree upon the electrode sepa- 
ration. For pressures of atmospheric, or below, the transi- 
tion is not abrupt but is completely defined. As the r.m.s. 
current is increased from low values to values near the 
critical value, the number of cycles for which the reignition 
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potential i is on the higher characteristic decreases. For cur- 
rents equal to, or greater than the critical value, only the 
low reignition potentials are obtained. This transition is 
accompanied by a change in the burning voltage of the arc. 
The factors influencing reignition are discussed and pos- 
sible mechanisms are considered. 

* Kennecott Wire and Cable Company. 


6. Townsend Ionization Coefficients in Argon Photo- 
Tubes. W.S. Huxrorp, Northwestern University.—Values 
of the ionization coefficient a have been obtained in Cs-Ag-O 
tubes with parallel plates in which the cathodes were acti- 
vated by the usual commercial process. Six tubes having 
electrode separations of 2, 4, 6, 8, 12 and 16 mm were em- 
ployed containing argon at a pressure of 0.175 mm of 
mercury. The P.D. values lie in the neighhorhood of the 
minimum sparking potential and the measurements cover 
a range which includes the maximum value of a/E plotted 
as a function of E/p. The results are in substantial agree- 
ment with other recent measurements of the electron ioniza- 
tion coefficient in highly purified argon. Calculations of the 
coefficient y, which represents the number of electrons 
leaving the cathode for each positive ion reaching it, yield 
values which are considerably greater than those obtained 
with cathodes having a higher electron work function. In 
this calculation account must be taken of an increased 
photo-sensitivity of the cathode which occurs at high gas 
amplification and is apparently induced by the action of 
positive ions. Sparking potential measurements for pure 
argon in tubes having nickel cathodes, nickel cathodes acti- 
vated with barium, and Cs-Ag-O cathodes show a progres- 
sive lowering of the minimum sparking potential and a 
shift of this minimum to lower P.D. values as the work 
function of the cathode decreases. 


7. Ultrasonic Reflection Losses in Air, Helium and CO;.* 
R. S. ALLEMAN AND J. C. HuBBARD, Johns Hopkins Uni- 
versity —The new heat conduction theory of acoustic 
reflection developed by Herzfeld! requires a modification 
of the theory of the acoustic resonator interferometer.’ 
As a result the values of +’, the intensity reflection coeffi- 
cient, first obtained by Hubbard,? and those of Curtis,’ 
showing a dependence upon frequency, have been recalcu- 
lated, and together with new values by the authors are 
given below, showing excellent general agreement with 
Herzfeld’s theory. 


Author Gas kc/sec Exp Theo 
1.H Air 597 0.034 0.034 
2.A 1237 051 
189 015 0122 
4.H $97 0148 0216 
a © He 252 .048 6 

‘ 6. A 792 .089 O81 
> 1001 085 092 
8. A 1192 101 097 
9A 1976 124 
* To be read by title. 
F. Herzfeld, Phys. Rev. 53, 899 (19 


K. 
ae Hubbard, Phys. Rev. 38, 1011 (1930); 41, 523 (1932); 46, 525 


8. An X-Ray Test of Superstructure in FeNi;. F. E. 
HawortH, Bell Telephone Laboratories, Inc.—A search has 
been made for superstructure in permalloy of composition 


near FeNi;. Curves of x-ray scattering power near the ab- 
sorption edges of Ni and Fe show that FeX@ radiation 
should give the strongest superstructure lines, and calcula- 
tions indicate that reflection from (321) planes of a per- 
fectly ordered structure should have about one-seventh the 
intensity of the ordinary (222) reflection. The specimen was 
a powdered alloy containing 70 percent nickel, annealed at 
1000°C and baked at 425°C. It was placed in a focusing 
camera in a beam of FeX§ rays which were selected and 
concentrated by a curved rocksalt crystal, and exposed for 
one hundred hours. No superstructure lines appeared, al- 
though superstructure in 8-brass was easily detected with 
similar technique. The negative result indicates that no 
long range order exists in the region near FeNis. 


9. The Spectrum of Doubly Ionized Lead, Pb III.* M. F. 
CraAwrorp, A. B. McLay A. M. Crooxker, University 
of Toronto.—The spectrum of lead has been measured in 
the region 800A to 10,000A. The lines were excited by an 
electrodeless discharge in lead vapor. A series of spectro- 
grams of the entire region were taken for different excita- 
tion conditions in the source, and the variation of intensity 
of a line with excitation showed, in general, to which spec- 
trum it belonged. About four hundred lines have been classi- 
fied as transitions between ninety levels of the following 
configurations: 6sns, n=6 to 10; 6snp, n=6 to 8; 6snd, 
n=6 to 9; 6snf, n=5 and 6; 6sng, n=5 to 8; 6snh, n=6 to 
9; 6p"; Op 7s; 6p 6d; Sd® 6s? 6p. The hyperfine structures of 
sixty-four of the classified lines have been measured with a 
twenty-one foot grating by Crooker.'! The observed struc- 
tures of these are consistent with this multiplet analysis. 
The 6sng configurations show interesting anomalous multi- 
plet and hyperfine structures. A theoretical interpretation 
of these will be published soon by ¢ one of the authors 
(M. F. C.). 


* To be read by title. 
1Can. J. Research 14, 115 (1936). 


10. The Number of Lines in a Series as a Function of 
Electron Concentration. Frep L. Monver, National 
Bureau of Standards.—In the spectrum of an intense dis- 
charge the successive lines of a series are increasingly 
broadened until they merge into a continuous spectrum at 
an effective quantum number which can be quite definitely 
evaluated. Spectra of capillary discharges under conditions 
of nearly complete ionization have been obtained with elec- 
tron concentrations, N,, ranging from 5 X10" to 1.3 X10"*. 
In this range the maximum effective quantum number 
for the violet component of the diffuse series decreases 
from 17.0 to 8.0. Values are about 0.3 greater for the red 
component and 0.5 greater for the F series. A plot of log 
versus log N, givesa straight line with a slope indicating that 
n~7is proportional to N,. Spectra of low current discharges 
give n about 25 for N,=10" but values seem to depend on 
vapor pressure and resolving power for N, less than 10", 
Sugita has shown that if the number of states is limited by 
electron collisions then n~? is proportional to N,. On the 
view that determines the number of discrete states it fol- 
lows that the ionization potential is reduced by about 0.2 
volt when N,=10"*, 
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11. The Effect of Magnetic Fields on Persistent Cur- 
rents in a Closed Superconducting Circuit. K. C. MANN, 
H. Grayson SMITH AND J. O. WILHELM, University of 
Toronto.—A study has been made of the effect of external 
magnetic fields in preventing persistent currents from flow- 
ing in a closed superconducting circuit, in order to deter- 
mine as closely as possible the field strength at which the 
first traces of electrical resistance appear. In agreement 
with earlier work by a less accurate method,! and with the 
resistance measurements of de Haas,? but in contradiction 
to the idealized theory of the transition state,’ it is clearly 
shown that no resistance appears until the mean magnetic 
induction within the specimen becomes considerable. It has 
been found further that the amount of resistanceless cur- 
rent which the specimen can carry under these conditions is 
limited. The results are greatly affected by the addition of 
small amounts of impurity. When the current strength was 
slightly greater than its equilibrium value, it required a 
measurable time to decay, the behavior being similar to 
the relaxation effects observed in the magnetic and thermal 
properties of superconductors. Several authors have sug- 
gested that the relaxation in magnetic properties is due to 
the decay of eddy currents. A thermodynamic argument 
now shows that relaxation in the thermal properties will be 
caused in the same way. 


1K. C. Mann and J. O. Wilhelm, Trans. Roy. Soc. Canada 31, 19 


(1937). 
2W. J. de Haas, J. Voogt and J. M. Jonker, Physica 1, 281 (1934). 
3 F. London, Physica 3, 450 (1936); Nature 137, 991 (1936). 


12. A Theoretical Lower Bound for the Energy of the 
Normal State of Helium. A. F. STEVENSON AND M. F. 
CrawrorD, University of Toronto—A method of using 
Weinstein’s expression! for a lower bound for the energy, 
as given recently by one of us,? has been applied to the 
normal state of heljum, using the same sequence of func- 
tions as was used by Hylleraas in his calculation of an 
upper bound. The successive approximations obtained by 
the lower bound method shaw a satisfactory convergence 
towards the upper bound finally obtained by Hylleraas. 


1D. H. Weinstein, Proc. Nat. Acad. Sci. 20, 529 (1934). 
2A. F. Stevenson, Phys. Rev. 53, 199 (1938). 


13. Experiments and Theory on the Performance of the 
Geiger Point Counter. J. Lioyp BOHN AND JOSEPH 
MorGan, Temple University.—Preparation of sensitive and 
insensitive points, why rough points or very sharp ones will 
not work properly, and why a point and plane combination 
may cease to function as the point approaches the plane 
will be discussed. Curves showing operating characteristics 
of various combinations of points, planes and tubes at 
atmospheric pressure and for point and tube combination 
at various pressures will be shown—also curves showing 
the variation in time for the reappearance of counts, after 
a Bakelite cap on the front of the tube had been touched, 
with voltage. That a modified form of Zeleny’s! theory 
suffices to explain the action of the Geiger point counter 


is concluded. 
1 Phys. Rev. 15, 566 (1922). 


14. A Vacuum Tube High Voltage Stabilizing Circuit 
for d.c. Sources. H. L. Scuuttz,* Yale University—When 


the stabilizing circuit to be described is applied to a con. 
ventional rectified-filtered source, initial voltage variations 
are reduced to less than 0.03 percent and the apparent 
internal resistance viewed from the output is of the order 
of a few thousandths of an ohm. Continuous control of 
output voltage from several hundred to several thousand 
volts at currents as high as several hundred milliamperes js 
possible using tubes of the receiving type. The circuit js 
similar to one developed by C. D. Bock, formerly of this 
laboratory, and has the following principle of operation, 
A triode or group of triodes connected in parallel is em- 
ployed as a variable impedance in the positive line from 
the source to be stabilized. Residual voltage variations 
across part of a bleeder resistance in the output are am- 
plified and actuate the grid of the control tube or tubes in 
proper phase with respect to the incoming variations, 
Simplicity and an actual improvement in operation js 
attained when plate and screen voltages for this amplifier 
are derived from taps on the common bleeder resistance, 
Instability occurs only if the load has appreciable reactance, 


* Sterling Fellow. 


15. Threshold Energies for Some (f,n) Reactions. 
M. G. Wuirte, L. A. DELSAsso, RUBBY SHERR AND L., N, 
RIpENOUR, Princeton University —Measurements of yield 
as a function of energy have been made by the method of 
stacked foils for several nuclear reactions of the (p,) type. 
The maximum incident proton energy is 4.2 Mev; the ion 
currents employed are about 1 microampere. The results 
to date indicate that the energy thresholds of three reac- 
tions studied are as follows: 


Half-life of Threshold energy, 
Reaction Radioelement Mev 
Ni®°(p,y)Cu®! or 
Ni®(p,n)Cu®! 3.4 hours <1.9 
Ni®(p,n) Cus 12.8 hours 21 
Ag?(p,n)Cd? 6.7 hours 2.4 


16. The Dirac Equation with Spin Greater Than 1/2.* 
ARTHUR BRAMLEY, Washington, D. C.—The interesting 
feature of the equations for electron with spin greater than 
1/2 and rest mass different from zero is that besides the 
Hamiltonian there exist linear equations of constraint 
acting on the wave functions. Their invariance under 
reflection demands that the two numbers / and & character- 
izing the spin be identical. From these two sets of four 
linear equations second order_differential equations can be 
constructed. In order to make the spin components definite 
we shall take the electron as propagated in the direction 
of a constant magnetic field. The second order equation 
corresponding to the relativistic Pauli equation*® shows that. 
the magnetic interaction is proportional to the velocity 
vector, not, as we might expect, the spin vector. Another 
consequence of the equations is that the matrices a must 
be expressed in terms of hypercomplex numbers such that 
a,*=a,?=a2=k(k+1)/3. In that case it follows from the 
equations of constraint that 


(3+R)k+R h u'v'* + |2 
pt) ~ — 
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in Dirac’s notation. This extra condition among the com- 
ponents of the momentum requires that for every velocity 
there shall be a characteristic rest mass.* 


* To be read by title. 

1 Dirac, Proc. Roy. Soc. 155, 447 (1936). 

? Sakata and Yukawa, Japan Proc. Phys, Math. Soc. 19, 91 (1937). 
3 Bramley, Science 87, 281 (1938). 


17. Virtual Levels of Be* and the Alpha-Particle Model. 
Jenny E, RosENTHAL, Brooklyn College.—It is known ex- 
perimentally that the normal state of Be® is barely stable 
and that there is a broad virtual level at about 3 Mev. A 
theoretical calculation has been made assuming a “‘molec- 
ular” model with a Morse potential bet ween the a-particles 
adjusted to a 10 Mev depth of the potential well. The wave 
functions ¥ were evaluated for the continuum region; only 
the order of magnitude of the moment of inertia and not 
its exact value is important in this calculation. Theoreti- 
cally we can speak of a virtual state if, for a given value of 
the angular momentum, the ratio of the amplitude of y 
inside the potential well to that outside reaches an appre- 
ciable maximum for some energy range. The calculations 
show that no maximum exists for any value of the energy; 
a consideration of the Coulomb barrier should not modify 
the results appreciably. We conclude that the molecular 
model does not account for the existence of virtual levels. 


18. The $-Ray Spectrum of I'*. R. H. Bacon, E. N. 
GrISEWOOD AND C. W. VAN DER MERWE, New York Uni- 
versity, Washington Square College—The 8-ray spectrum 
of 1 has been measured in a Wilson cloud chamber. A 
K.U. plot shows an abrupt change of slope indicating two 
intersecting straight lines. If the assumption is made that 
a single straight line means the existence of a single group 
our results would indicate two groups with end points at 
2.1 Mev and 1.05 Mev, respectively. The end point of the 
high energy group agrees with that found by Alichanian, 
Alichanow and Dzelepow,* but the literature, to our 
knowledge, makes no mention of the other group. Details 
of the method used will be given at the meeting. 


* Physik. Zeits. Sowjetunion 10, 78-102 (1936). 


19. The Radioactive Isotope of Samarium. T. R. WIL- 
kins, University of Rochester AND A. J. DEMPSTER, Univer- 
sity of Chicago.—G. Hevesy (1934) reported that samarium 
was a natural a-emitter with an apparent half-life of 1.2 
X10" years. The real half-life can be calculated only if the 
active isotope can be identified and its percentage abun- 
dance known. Various suggestions have been made as to 
this isotope, some based on hyperfine structure abnormal- 
ities and some on relations between decay constant, nuclear 
radius and abundance ratios. We have devised a new tool 
for such studies by combining the mass spectrograph and 
the photographic emulsion track technique. A direct 
deposit of the isotopes of samarium was made on an Ilford 
R type emulsion in the mass spectrograph. The plate was 
set aside for several months and then developed and 
examined by the microscope. These tracks lie in the emul- 
sion, and are seen by focusing on the emulsion through the 
half-transparent samarium deposited on the gelatin super- 


coating of this type of plate. The typical a-ray track 
(1.13 em air range and a-ray spacing) was found under 
the isotope of Sm'*. Schumann plates have been used for 
calibration but these are actually unnecessary since the 
numbers of ‘‘spots’’ comprising the visible deposits corre- 
spond to the abundance ratios of the isotopes. Taking the 
abundance of Sm"™* as 14 percent, the real half-life will be 
14 percent of the apparent half-life. 


20. Altitude Variation of the Relative Frequency of 
Occurrence of Pairs and Larger Cosmic-Ray Showers. 
Darot K. Froman, Macdonald College, McGill University 
AND J. C. STEARNS, University of Denver.—The ratio of 
triple to quadruple frequency of coincidence of Geiger- 
Miiller counters, arranged to detect showers containing 
two or more and three or more ionizing shower rays, 
respectively, has been measured for various thicknesses of 
scatterer at two altitudes, viz. 36 meters and 4300 meters 
above sea level. In accordance with the predictions of the 
multiplicative theory of shower production, this ratio 
decreases rapidly with increase in the thickness of the 
scatterer. For very thin scatterers, the ratio is much 
greater at the higher altitude. For very thin scatterers it is 
expected, on the theory, that this ratio would be greater 
for showers produced by photons than it would for electron- 
induced showers. Experimental data on this point are 
being taken, and the preliminary results tend to confirm 
this view. The marked variation of the ratio with altitude 
suggests the possibility that the relative number of 
photons and electrons in the cosmic-ray complex is a 
function of altitude. Experiments, designed to test this 
point, are being initiated this summer. 


21. The Michelson-Morley-Miller Experiment. W. 8. 
CARTMEL, Université de Montréal.—According to Lorentz 
t=/y*/c8-(1+cos?—8@). In an oblique interferometer we 
have cos? @ in one arm and sin? (@—@) in the other and for 
their difference, 


20—sin & sin 2(0+ ]. (1) 


But to find the velocity along the sides of the triangle we 
must consider the three rays which meet at the plate. 
From Huygens’ principle V,/sin = Va-/sin ig: = V,/sinr 
from which we derive sin ig. As shown by 
Lodge,! ig: =ia+v*/c i sin 2¢, from which we derive 
for the reflected ray Va = V.[1+ }e/ce(1—sin Q) cot 
Xcos ], and finally 


Q) cot cos 2(0+ 32) 
—cos 20+sin Q sin 2(04+}32)). (2) 


Since the obliquity is slight in the Michelson-Morley- Miller 
interferometer, being only 6}°, 


Disregarding the limits we have a perfect sine curve, but 
taking them into account the fringe shift curve has maxima 
and minima alternately 63° and 83}° apart, which would 
lead to hit or miss results taking measurements every 224°, 
and this was found experimentally by Professor Miller.? 
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With a 45° interferometer the maxima and minima would 
be 45° apart. It is evident that with the very oblique inter- 
ferometer of Kennedy, only a very small part of the actual 
fringe shift would be found taking measurements every 
90°, which accounts for the smallness of the results found 
by Illingworth.* Eq. (2), however, is believed to be only 
approximate, the exact solution being 


6F=21/d sin Qv?/c? sin 2(0+ 3), 


an equation which may be derived from the integrals of 
Lorentz, also by the method given in New York two years 

1 Sir Oliver Lodge, Phil. Trans. 184A, 727 (1894). 

2 Dayton C. Miller, Rev. Mod. Phys. 5, 203 (1933). 


3K. K. Illingworth, Phys. Rev. 30, 692 (1927). 
4W. B. Cartmel, Phys. Rev. 49, 647, 649 (1936). 


22. A Cloud-Chamber Study of the Alpha-Rays of 
Actino-Uranium. T. R. WiLkins AND D. P. CRAWForD, 
University of Rochester —W. M. Rayton and T. R. Wilkins! 
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(1935) reported an analysis of the alpha-rays from a thin}, 
evaporated film of metallic uranium in which complete 
resolution of the UI and UII isotopes was secured, J, 
addition, there seemed to be evidence of a small group of 
alpha-rays, presumably due to actino-uranium, betwee 

the main peaks due to UI and U II. In the present wor 

minor improvements have been made in the apparati 

and a much larger number of tracks studied. The me; 

ranges check those reported (by R and W) for UI ana 
U IT (2.625+0.0020 and 3.179+0.0024). The main peaks 
are so well fitted by Gaussian curves that the residuum 
is unquestionably present. The mean range of actino- 
uranium appears to be nearer that of U II. The relative 
fewness of these tracks prevents high accuracy of range 
determination but it is felt that a mean air-range (15°C) 
of 2.99 cm will prove to be very close to the correct value, 
The determination of the relative activity is correspond- 
ingly uncertain. The number of tracks measured was 2140, 


1 Phys. Rev. 51, 818-25 (1937). 
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